Strain energy in enzyme–substrate binding: An energetic insight into the flexibility versus rigidity of enzyme active site

Sunsetting Binding MOAD with its last data update and the addition of 3D-ligand polypharmacology tools

Binding MOAD is a database of protein-ligand complexes and their affinities with many structured relationships across the dataset. The project has been in development for over 20 years, but now, the time has come to bring it to a close. Currently, the database contains 41,409 structures with affinity coverage for 15,223 (37%) complexes. The website BindingMOAD.org provides numerous tools for polypharmacology exploration. Current relationships include links for structures with sequence similarity, 2D ligand similarity, and binding-site similarity. In this last update, we have added 3D ligand similarity using ROCS to identify ligands which may not necessarily be similar in two dimensions but can occupy the same three-dimensional space. For the 20,387 different ligands present in the database, a total of 1,320,511 3D-shape matches between the ligands were added. Examples of the utility of 3D-shape matching in polypharmacology are presented. Finally, plans for future access to the project data are outlined.

Computational methods for calculation of protein-ligand binding affinities in structure-based drug design

Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.

Molecular conversion of MIG6 hotspot-3 peptide from the nonbinder to a moderate binder of HER2 by rational design of an orthogonal interaction system at the HER2-peptide interface

Human epidermal growth factor receptor 2 (HER2) has been established as an approved druggable target for the treatment of patients with diverse gynecological tumors such as ovarian, cervical and breast cancers. The mitogen-inducible gene 6 (MIG6) protein is a negative regulator of HER2 signaling by using its Seg1 segment to disrupt the allosteric dimerization of HER2 kinase domain. Previous studies found that the Seg1 adopts three separated hotspots to interact with the HER2 dimerization interface, in which the third hotspot (H3) is located at the core region of the interface but its derived H3 peptide (356PKYVS360) and Tyr358Phe mutant (356PKFVS360) cannot bind effectively to the interface in an independent manner. In this study, we demonstrate that the H3 peptide can be converted from nonbinder to a moderate binder of HER2 by just adding an orthogonal noncovalent interaction system (X⋯O┄H) between a halogen bond (X⋯O) and a hydrogen bond (H┄O) involving peptide Phe358 residue and HER2 Val948/Trp951 residues. High-level calculations are utilized to rigorously characterize and rationally design the X⋯O┄H system, which is then optimized with different halogen atoms and at different substituting positions. It is revealed that there is a synergistic effect between the X⋯O and H┄O of the orthogonal interaction system; formation of the halogen bond can enhance the interaction strength of the hydrogen bond. In silico analysis and in vitro assay reach a consistence that Br-substitution at the m-position of peptide Phe358 phenyl moiety is the best choice that can render strong interaction for the X⋯O┄H system, which also makes the peptide 'bindable' to HER2 kinase domain, while F/Cl/I-substitution at the same position can only improve the peptide affinity moderately or modestly. In contrast, the Br-substitution at the o- and p-positions of peptide Phe358 phenyl moiety cannot define effective X⋯O┄H interaction and thus does not confer additional affinity to the HER2-peptide complex.

Rational Design and Intramolecular Cyclization of Hotspot Peptide Segments at YAP–TEAD4 Complex Interface

Background:

The Yes-Associated Protein (YAP) is a central regulator of Hippo pathway involved in carcinogenesis, which functions through interaction with TEA Domain (TEAD) transcription factors. Pharmacological disruption of YAP–TEAD4 complexes has been recognized as a potential therapeutic strategy against diverse cancers by suppressing the oncogenic activity of YAP.

Objective:

Two peptides, termed PS-1 and PS-2 are split from the interfacial context of YAP protein. Dynamics simulations, energetics analyses and fluorescence polarizations are employed to characterize the intrinsic disorder as well as binding energy/affinity of the two YAP peptides to TEAD4 protein.

Methods:

Two peptides, termed PS-1 and PS-2 are split from the interfacial context of YAP protein. Dynamics simulations, energetics analyses and fluorescence polarizations are employed to characterize the intrinsic disorder as well as binding energy/affinity of the two YAP peptides to TEAD4 protein.

Result:

The native conformation of PS-2 peptide is a cyclic loop, which is supposed to be constrained by adding a disulfide bond across the spatially vicinal residue pair Arg87-Phe96 or Met86- Phe95 at the peptide’s two ends, consequently resulting in two intramolecular cyclized counterparts of linear PS-2 peptide, namely PS-2(cyc87,96) and PS-2(cyc86,95). The linear PS-2 peptide is determined as a weak binder of TEAD4 (Kd = 190 μM), while the two cyclic PS-2(cyc87,96) and PS-2(cyc86,95) peptides are measured to have moderate or high affinity towards TEAD4 (Kd = 21 and 45 μM, respectively).

Conclusion:

PS-1 and PS-2 peptides are highly flexible and cannot maintain in native active conformation when splitting from the interfacial context, and thus would incur a considerable entropy penalty upon rebinding to the interface. Cyclization does not influence the direct interaction between PS-2 peptide and TEAD4 protein, but can largely reduce the intrinsic disorder of PS-2 peptide in free state and considerably minimize indirect entropy effect upon the peptide binding.

Effects of the Hydrophobicity of Key Residues on the Characteristics and Stability of Glucose Oxidase on a Graphene Surface

Glucose oxidase (GOx) is one of the most widely investigated enzymes in the field of bioelectrochemistry. It is mainly used for the detection of glucose in solutions and enzyme-based biofuel cells. On the basis of the combination of GOx with graphene, novel nanodevices exceeding conventional limits can be developed. To develop a hybrid enzyme-graphene nanodevice with a good performance, it is important that GOx is deposited well on the graphene surface while maintaining its structure and not impeding the oxidation activity of the GOx. In this study, we propose a method to improve the stability of GOx and secure its immobility on the graphene sheet and its glucose-binding affinity by single-point mutation of GOx using molecular dynamics simulations. We confirm that the structural stability, immobility, and substrate binding affinity of GOx can be modified by changing the hydrophobicity of a key residue. We demonstrate that biosensors or biofuel cells can be redesigned and their properties can be improved by using molecular dynamics simulation.

Second-generation 4,5,6,7-tetrahydrobenzo[d]thiazoles as novel DNA gyrase inhibitors

Aim: DNA gyrase and topoisomerase IV are essential bacterial enzymes, and in the fight against bacterial resistance, they are important targets for the development of novel antibacterial drugs. Results: Building from our first generation of 4,5,6,7-tetrahydrobenzo[d]thiazole-based DNA gyrase inhibitors, we designed and prepared an optimized series of analogs that show improved inhibition of DNA gyrase and topoisomerase IV from Staphylococcus aureus and Escherichia coli, with IC₅₀ values in the nanomolar range. Importantly, these inhibitors also show improved antibacterial activity against Gram-positive strains. Conclusion: The most promising inhibitor, 29, is active against Enterococcus faecalis, Enterococcus faecium and S. aureus wild-type and resistant strains, with minimum inhibitory concentrations between 4 and 8 μg/ml, which represents good starting point for development of novel antibacterials.

Universality of fold-encoded localized vibrations in enzymes

Enzymes speed up biochemical reactions at the core of life by as much as 15 orders of magnitude. Yet, despite considerable advances, the fine dynamical determinants at the microscopic level of their catalytic proficiency are still elusive. In this work, we use a powerful mathematical approach to show that rate-promoting vibrations in the picosecond range, specifically encoded in the 3D protein structure, are localized vibrations optimally coupled to the chemical reaction coordinates at the active site. Remarkably, our theory also exposes an hithertho unknown deep connection between the unique localization fingerprint and a distinct partition of the 3D fold into independent, foldspanning subdomains that govern long-range communication. The universality of these features is demonstrated on a pool of more than 900 enzyme structures, comprising a total of more than 10,000 experimentally annotated catalytic sites. Our theory provides a unified microscopic rationale for the subtle structure-dynamics-function link in proteins.

Systematic profiling of substrate binding response to multidrug-resistant mutations in HIV-1 protease: Implication for combating drug resistance

Human immunodeficiency virus 1 (HIV-1) protease (PR) represents one of the primary targets for developing antiviral agents for the treatment of HIV-infected patients. However, a number of multidrug-resistant mutations in the enzyme have been observed over the past decades, largely limiting the application of PR inhibitors in antiviral therapy. A systematic investigation of the intermolecular interaction between the multidrug-resistant mutants of HIV-1 PR and its substrates would help to establish a complete profile of substrate response to PR mutations and to design new antiviral agents combating drug resistance. Here, we describe an integrative method to profile 6 clinical multidrug-resistant PR mutants against a panel of 16 substrate octapeptides that flank 12 distinct PR cleavage sites in viral precursor polyproteins. It is found that most multidrug-resistant mutations have only a modest or moderate effect on substrate peptide binding, although these mutations would cause a large free energy loss in PR inhibitor binding. Structural analysis reveals that the substrate peptides are loosely bound within PR active pocket to form a wide contact interface between them, and thus mutation of just single or few residues seems not to influence PR-substrate binding considerably. In addition, peptides derived from variable cleavage sites are generally more sensitive to the mutations as compared to those derived from conserved sites, supporting the co-evaluation mechanism of HIV-1 PR and its substrates under drug suppression. We also identify 12 functionally conserved key residues around the enzyme's active site, which play crucial role in substrate recognition. In vitro fluorescence anisotropy assays confirm that wild-type PR can bind substrate peptides ARVL/AEAM and NLAF/PQGE with a moderately high affinity (K_D=2 and 16μM, respectively). In contrast, the key residue mutations N25D/D29N can completely eliminate (K_D=n.d.) or largely reduce (K_D=32 and 120μM, respectively) the binding capability of the two peptides, suggesting that these PR residues could be the potential target sites for developing resistance-free anti-HIV agents.

Probing the influence of carboxyalkyl groups on the molecular flexibility and the charge density of apigenin derivatives

Apigenin is an important flavonoids due to its antidiabetic bioactivity. It was reported experimentally that the 7-substituent derivative of apigenin has higher biological activity than 4'- and 5-substituted derivatives while introducing sole carboxyalkyl group -(CH₂)₇COOH into the parent structure. Molecular docking studies indicated that the other two derivatives have lower binding affinities than the 7-substituent derivative (-7.52 kcal mol^-1), which is considered to be a better inhibitor than the parent molecule. Almost all of the carbon atoms and oxygen atoms are coplaner for all three molecules in solution phase, however, all carboxyalkyl groups bend inside into the parent molecules in the active site, and the jagged geometries of the carbon chains are destroyed correspondingly. In addition, most of the electron densities of the chemical bonds for all molecules are decreased, especially the 7-substituent derivative. In contrast, most of the Laplacian values for three molecules are increased in the active site, which suggests that the charge densities at the bond critical point (bcp) are much more depleted than the solution phase. Dipole moments of derivatives are all increased in the active site, suggesting strong intermolecular interactions. After interacting with the S. cerevisiae α-glucosidase, only the 7-substituent derivative has the lowest energy gap ΔE _HOMO-LUMO, which indicates the lowest stability and the highest inhibition activity. Graphical abstract Probing the influence of carboxyalkyl groups on the molecular flexibility and the charge density of apigenin derivatives.

Computational analysis and enzyme assay of inhibitor response to disease single nucleotide polymorphisms (SNPs) in lipoprotein lipase

Lipoprotein lipase (LPL) is the rate-limiting enzyme for the hydrolysis of the triglyceride (TG) core of circulating TG-rich lipoproteins, chylomicrons, and very low-density lipoproteins. The enzyme has been established as an efficacious and safe therapeutic target for the management of obesity. Here, a systematic profile of the lipase inhibitor response of three anti-obesity agents (Orlistat, Lipstatin, and Cetilistat) to clinical LPL missense mutations arising from disease single nucleotide polymorphisms (SNPs) was established by integrating complex structure modeling, virtual mutagenesis, molecular dynamics (MD) simulations, binding energy analysis, and radiolabeled TG hydrolysis assays. The profile was then used to characterize the resistance and sensitivity of systematic mutation–inhibitor pairs. It is suggested that the Orlistat and Lipstatin have a similar response profile to the investigated mutations due to their homologous chemical structures, but exhibit a distinct profile to that of Cetilistat. Most mutations were predicted to have a modest or moderate effect on inhibitor binding; they are located far away from the enzyme active site and thus can only influence the binding limitedly. A number of mutations were found to sensitize or cause resistance for lipase inhibitors by directly interacting with the inhibitor ligands or by indirectly addressing allosteric effect on enzyme active site. Long-term MD simulations revealed a different noncovalent interaction network at the complex interfaces of Orlistat with wild-type LPL as well as its sensitized mutant H163R and resistant mutant I221T.

Structure-Based Design of a Br Halogen Bond at the Complex Interface of the Human Placental HtrA1 PDZ Domain with Its Heptapeptide Ligand

The shock-induced serine protease HtrA1 is a potential regulator of human placenta development during pregnancy. The protein contains a functional PDZ domain that has been solved in complex with a phage display-derived heptapeptide: Asp-6 Ser-5 Arg-4 Ile-3 Trp-2 Trp-1 Val0 . In this study, a rationally designed halogen bond was introduced to the domain-peptide complex based on its NMR structure in solution. We computationally compared the stabilization energies and hindrance effects due to the presence of different halogens X (X = F, Cl, Br, or I), using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach, and found that the Br atom could considerably promote the peptide binding free energy (ΔΔG = -5.2 kcal/mol). Fluorescence assays confirmed that the peptide affinity to the HtrA1 PDZ domain was improved by approximately sevenfold upon bromination. Structural analysis identified a geometrically perfect halogen bond between the Br atom of the peptide Trp-1 residue and the carbonyl O atom of the HtrA1 Ile385 residue, with a bond length and an interaction energy of d = 3.20 Å and ΔE = -3.7 kcal/mol, respectively.

Targeting druggable enzymome by exploiting natural medicines: An in silico-in vitro integrated approach to combating multidrug resistance in bacterial infection

CONTEXT

Antibiotic resistance is a major clinical and public health problem. Development of new therapeutic approaches to prevent bacterial multidrug resistance during antimicrobial chemotherapy has thus been becoming a primary consideration in the medicinal chemistry community.

OBJECTIVE

We described a new strategy that combats multidrug resistance by using natural medicines to target the druggable enzymome (i.e., enzymatic proteome) of Staphylococcus aureus.

MATERIALS AND METHODS

A pipeline of integrating in silico analysis and in vitro assay was purposed to identify antibacterial agents from a large library of natural products with diverse structures, high drug-likeness, and relatively low flexibility, with which a systematic interactome of 826 natural product candidates with 125 functionally essential S. aureus enzymes was constructed via a high-throughput cross-docking approach. The obtained docking score matrix was then converted into an array of synthetic scores; each corresponds to a natural product candidate. By systematically examining the docking results, a number of highly promising candidates with potent antibacterial activity were suggested.

RESULTS

Three natural products, i.e., radicicol, jorumycin, and amygdalin, have been determined to possess strong broad-spectrum potency combating both the drug-resistant and drug-sensitive strains (MIC value <10 μg/ml). In addition, some natural products such as tetrandrine, bilobalide, and arbutin exhibited selective inhibition on different strains.

DISCUSSION AND CONCLUSION

Combined quantum mechanics/molecular mechanics analysis revealed diverse non-bonded interactions across the complex interfaces of newly identified antibacterial agents with their putative targets GyrB ATPase and tyrosyl-tRNA synthetase.

Cutoff lensing: predicting catalytic sites in enzymes

Predicting function-related amino acids in proteins with unknown function or unknown allosteric binding sites in drug-targeted proteins is a task of paramount importance in molecular biomedicine. In this paper we introduce a simple, light and computationally inexpensive structure-based method to identify catalytic sites in enzymes. Our method, termed cutoff lensing, is a general procedure consisting in letting the cutoff used to build an elastic network model increase to large values. A validation of our method against a large database of annotated enzymes shows that optimal values of the cutoff exist such that three different structure-based indicators allow one to recover a maximum of the known catalytic sites. Interestingly, we find that the larger the structures the greater the predictive power afforded by our method. Possible ways to combine the three indicators into a single figure of merit and into a specific sequential analysis are suggested and discussed with reference to the classic case of HIV-protease. Our method could be used as a complement to other sequence- and/or structure-based methods to narrow the results of large-scale screenings.

Identification of new members of the MAPK gene family in plants shows diverse conserved domains and novel activation loop variants

Background

Mitogen Activated Protein Kinase (MAPK) signaling is of critical importance in plants and other eukaryotic organisms. The MAPK cascade plays an indispensible role in the growth and development of plants, as well as in biotic and abiotic stress responses. The MAPKs are constitute the most downstream module of the three tier MAPK cascade and are phosphorylated by upstream MAP kinase kinases (MAPKK), which are in turn are phosphorylated by MAP kinase kinase kinase (MAPKKK). The MAPKs play pivotal roles in regulation of many cytoplasmic and nuclear substrates, thus regulating several biological processes.

Results

A total of 589 MAPKs genes were identified from the genome wide analysis of 40 species. The sequence analysis has revealed the presence of several N- and C-terminal conserved domains. The MAPKs were previously believed to be characterized by the presence of TEY/TDY activation loop motifs. The present study showed that, in addition to presence of activation loop TEY/TDY motifs, MAPKs are also contain MEY, TEM, TQM, TRM, TVY, TSY, TEC and TQY activation loop motifs. Phylogenetic analysis of all predicted MAPKs were clustered into six different groups (group A, B, C, D, E and F), and all predicted MAPKs were assigned with specific names based on their orthology based evolutionary relationships with Arabidopsis or Oryza MAPKs.

Conclusion

We conducted global analysis of the MAPK gene family of plants from lower eukaryotes to higher eukaryotes and analyzed their genomic and evolutionary aspects. Our study showed the presence of several new activation loop motifs and diverse conserved domains in MAPKs. Advance study of newly identified activation loop motifs can provide further information regarding the downstream signaling cascade activated in response to a wide array of stress conditions, as well as plant growth and development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1244-7) contains supplementary material, which is available to authorized users.