In silico modeling of the specific inhibitory potential of thiophene-2,3-dihydro-1,5-benzothiazepine against BChE in the formation of β-amyloid plaques associated with Alzheimer's disease

<p>Abstract</p> <p>Background</p> <p>Alzheimer's disease, known to be associated with the gradual loss of memory, is characterized by low concentration of acetylcholine in the hippocampus and cortex part of the brain. Inhibition of acetylcholinesterase has successfully been used as a drug target to treat Alzheimer's disease but drug resistance shown by butyrylcholinesterase remains a matter of concern in treating Alzheimer's disease. Apart from the many other reasons for Alzheimer's disease, its association with the genesis of fibrils by β-amyloid plaques is closely related to the increased activity of butyrylcholinesterase. Although few data are available on the inhibition of butyrylcholinesterase, studies have shown that that butyrylcholinesterase is a genetically validated drug target and its selective inhibition reduces the formation of β-amyloid plaques.</p> <p>Rationale</p> <p>We previously reported the inhibition of cholinesterases by 2,3-dihydro-1, 5-benzothiazepines, and considered this class of compounds as promising inhibitors for the cure of Alzheimer's disease. One compound from the same series, when substituted with a hydroxy group at C-3 in ring A and 2-thienyl moiety as ring B, showed greater activity against butyrylcholinesterase than to acetylcholinesterase. To provide insight into the binding mode of this compound (Compound A), molecular docking in combination with molecular dynamics simulation of 5000 ps in an explicit solvent system was carried out for both cholinesterases.</p> <p>Conclusion</p> <p>Molecular docking studies revealed that the potential of Compound A to inhibit cholinesterases was attributable to the cumulative effects of strong hydrogen bonds, cationic-π, π-π interactions and hydrophobic interactions. A comparison of the docking results of Compound A against both cholinesterases showed that amino acid residues in different sub-sites were engaged to stabilize the docked complex. The relatively high affinity of Compound A for butyrylcholinesterase was due to the additional hydrophobic interaction between the 2-thiophene moiety of Compound A and Ile69. The involvement of one catalytic triad residue (His438) of butyrylcholinesterase with the 3'-hydroxy group on ring A increases the selectivity of Compound A. C-C bond rotation around ring A also stabilizes and enhances the interaction of Compound A with butyrylcholinesterase. Furthermore, the classical network of hydrogen bonding interactions as formed by the catalytic triad of butyrylcholinesterase is disturbed by Compound A. This study may open a new avenue for structure-based drug design for Alzheimer's disease by considering the 3D-pharmacophoric features of the complex responsible for discriminating these two closely-related cholinesterases.</p

Dynamic changes in the secondary structure of ECE-1 and XCE account for their different substrate specificities

<p>Abstract</p> <p>Background</p> <p>X-converting enzyme (XCE) involved in nervous control of respiration, is a member of the M13 family of zinc peptidases, for which no natural substrate has been identified yet. In contrast, it’s well characterized homologue endothelin-converting enzyme-1 (ECE-1) showed broad substrate specificity and acts as endopeptidase as well as dipeptidase. To explore the structural differences between XCE and ECE-1, homology model of XCE was built using the complex structure of ECE-1 with phosphoramidon (pdb-id: 3DWB) as template. Phosphoramidon was docked into the binding site of XCE whereas phosphate oxygen of the inhibitor was used as water molecule to design the apo forms of both enzymes. Molecular dynamics simulation of both enzymes was performed to analyze the dynamic nature of their active site residues in the absence and presence of the inhibitor.</p> <p>Results</p> <p>Homology model of XCE explained the role of non-conserved residues of its S2’ subsite. Molecular dynamics (MD) simulations identified the flexible transitions of F149/I150, N566/N571, W714/W719, and R145/R723 residues of ECE-1/XCE for the strong binding of the inhibitor. Secondary structure calculations using DSSP method reveals the folding of R145/R723 residue of ECE-1/XCE into β-sheet structure while unfolding of the S2’ subsite residues in aECE-1 and sustained compact folding of that of aXCE. The results evaluated are in good agreement with available experimental data, thus providing detailed molecular models which can explain the structural and specificities differences between both zinc peptidases.</p> <p>Conclusions</p> <p>Secondary structure changes of both enzymes during the simulation time revealed the importance of β-sheet structure of R145/R723 for its binding with the terminal carboxylate group of the inhibitor. Unfolding of the α-helix comprising the S2’ subsite residues in aECE-1 correlate well with its endopeptidase activity while their compact folding in aXCE may account for the inactivity of the enzyme towards large C-terminal containing substrates.</p