Human farnesyl pyrophosphate synthase inhibition by nitrogen bisphosphonates: a 3D-QSAR study

Role of Metabolism in Bone Development and Homeostasis

Carbohydrates, fats, and proteins are the underlying energy sources for animals and are catabolized through specific biochemical cascades involving numerous enzymes. The catabolites and metabolites in these metabolic pathways are crucial for many cellular functions; therefore, an imbalance and/or dysregulation of these pathways causes cellular dysfunction, resulting in various metabolic diseases. Bone, a highly mineralized organ that serves as a skeleton of the body, undergoes continuous active turnover, which is required for the maintenance of healthy bony components through the deposition and resorption of bone matrix and minerals. This highly coordinated event is regulated throughout life by bone cells such as osteoblasts, osteoclasts, and osteocytes, and requires synchronized activities from different metabolic pathways. Here, we aim to provide a comprehensive review of the cellular metabolism involved in bone development and homeostasis, as revealed by mouse genetic studies.

New insights into human farnesyl pyrophosphate synthase inhibition by second-generation bisphosphonate drugs

Pamidronate, alendronate, APHBP and neridronate are a group of drugs, known as second-generation bisphosphonates (2G-BPs), commonly used in the treatment of bone-resorption disorders, and recently their use has been related to some collateral side effects. The therapeutic activity of 2G-BPs is related to the inhibition of the human Farnesyl Pyrophosphate Synthase (hFPPS). Available inhibitory activity values show that 2G-BPs act time-dependently, showing big differences in their initial inhibitory activities but similar final IC₅₀ values. However, there is a lack of information explaining this similar final inhibitory potency. Although different residues have been identified in the stabilization of the R₂ side chain of 2G-BPs into the active site, similar free binding energies were obtained that highlighted a similar stability of the ternary complexes, which in turns justified the similar IC₅₀ values reported. Free binding energy calculations also demonstrated that the union of 2G-BPs to the active site were 38 to 54 kcal mol^-1 energetically more favourable than the union of the natural substrate, which is the basis of the inhibition potency of the hFPPS activity.

Physiologic Activity of Bisphosphonates – Recent Advances

Background:

Bisphosphonates are drugs commonly used for the medication and prevention of diseases caused by decreased mineral density. Despite such important medicinal use, they display a variety of physiologic activities, which make them promising anti-cancer, anti-protozoal, antibacterial and antiviral agents.

Objective:

To review physiological activity of bisphosphonates with special emphasis on their ongoing and potential applications in medicine and agriculture.

Method:

Critical review of recent literature data.

Results:

Comprehensive review of activities revealed by bisphosphonates.

Conclusion:

although bisphosphonates are mostly recognized by their profound effects on bone physiology their medicinal potential has not been fully evaluated yet. Literature data considering enzyme inhibition suggest possibilities of far more wide application of these compounds. These applications are, however, limited by their low bioavailability and therefore intensive search for new chemical entities overcoming this shortage are carried out.

Mechanistic insights into protonation state as a critical factor in hFPPS enzyme inhibition

Zoledronate and risedronate are the most powerful available nitrogen-containing bisphosphonates used in the treatment of bone-resorption disorders. Knowledge about inhibition mechanisms of these molecules is based on available crystallographic structures of human farnesyl pyrophosphate synthase (hFPPS). However, there is a lack of information explaining the inhibition potency of these two molecules compared to the natural substrate, dimethylallyl pyrophosphate. We carried out a molecular dynamics study that shown: (1) that NBPs potency is related to higher electrostatic interactions with the metallic cluster of the active site than to the natural substrate, and (2) the protonation of the R2 side chain is a critical factor to stabilize the NBPs into a closely irreversible ternary complex with the hFPPS.

Combined pharmacophore modeling, 3D-QSAR, homology modeling and docking studies on CYP11B1 inhibitors

The mitochondrial cytochrome P450 enzymes inhibitor steroid 11β-hydroxylase (CYP11B1) can decrease the production of cortisol. Therefore, these inhibitors have an effect in the treatment of Cushing’s syndrome. A pharmacophore model generated by Genetic Algorithm with Linear Assignment for Hypermolecular Alignment of Datasets (GALAHAD) was used to align the compounds and perform comparative molecular field analysis (CoMFA) with Q² = 0.658, R² = 0.959. The pharmacophore model contained six hydrophobic regions and one acceptor atom, and electropositive and bulky substituents would be tolerated at the A and B sites, respectively. A three-dimensional quantitative structure-activity relationship (3D-QSAR) study based on the alignment with the atom root mean square (RMS) was applied using comparative molecular field analysis (CoMFA) with Q² = 0.666, R² = 0.978, and comparative molecular similarity indices analysis (CoMSIA) with Q² = 0.721, R² = 0.972. These results proved that all the models have good predictability of the bioactivities of inhibitors. Furthermore, the QSAR models indicated that a hydrogen bond acceptor substituent would be disfavored at the A and B groups, while hydrophobic groups would be favored at the B site. The three-dimensional (3D) model of the CYP11B1 was generated based on the crystal structure of the CYP11B2 (PDB code 4DVQ). In order to probe the ligand-binding modes, Surflex-dock was employed to dock CYP11B1 inhibitory compounds into the active site of the receptor. The docking result showed that the imidazolidine ring of CYP11B1 inhibitors form H bonds with the amino group of residue Arg155 and Arg519, which suggested that an electronegative substituent at these positions could enhance the activities of compounds. All the models generated by GALAHAD QSAR and Docking methods provide guidance about how to design novel and potential drugs for Cushing’s syndrome treatment.