Pharmacological studies from our group [Lima et al. Pharmacol Biochem Behav 92:508, (2009)] revealed that geissospermine (GSP), the major alkaloid of the bark extract of Brazilian Geissospermum vellosii, inhibits acetylcholinesterases (AChEs) in the brains of rats and electric eels (Electrophorus electricus). However, the binding mode (i.e., conformation and orientation) of this indole-indoline alkaloid into the AChE active site is unknown. Therefore, in order to propose a plausible binding mode between GSP and AChE, which might explain the observed experimental inhibitory activity, we performed comparative automatic molecular docking simulations using the AutoDock and Molegro Virtual Docker (MVD) programs. A sample of ten crystal structures of the Pacific electric ray (Torpedo californica) TcAChE, in complex with ten diverse active site ligands, was selected as a robust re-docking validation test, and also for GSP docking. The MVD results indicate a preferential binding mode between GSP and AChE, in which GSP functional groups may perform specific interactions with residues in the enzyme active site, according to the ligand-protein contacts detected by the LPC/CSU server. Four hydrogen bonds were detected between GSP and Tyr121, Ser122, Ser200, and His440, in which the last two residues belong to the catalytic triad (Ser200···His440···Glu327). Hydrophobic and π-π stacking interactions were also detected between GSP and Phe330 and Trp84, respectively; these are involved in substrate stabilization at the active site. This study provides the basis to propose structural changes to the GSP structure, such as molecular simplification and isosteric replacement, in order to aid the design of new potential AChE inhibitors that are relevant to the treatment of Alzheimer's disease.
A series of 74 dihydroalkoxybenzyloxopyrimidines (DABOs), a class of highly potent non-nucleoside reverse transcriptase inhibitors (NNRTIs), was retrieved from the literature and studied by receptor-dependent (RD) three-dimensional quantitative structure-activity relationship (3D-QSAR) analysis to derive RD-3D-QSAR models. The descriptors in this new method are the steric and electrostatic interaction energies of the protein-ligand complexes (per residue) simulated by molecular dynamics, an approach named
Residue-Ligand Interaction Energy
(ReLIE). This study was performed using a training set of 59 compounds and the MKC-442/RT complex structure as reference. The ReLIE-3D-QSAR models were constructed and evaluated by genetic algorithm (GA) and partial least squares (PLS). In the best equations, at least one term is related to one of the amino acid residues of the p51 subunit:
Asn136
,
Asn137
,
Glu138
, and
Thr139
. This fact implies the importance of interchain interaction (p66-p51) in the equations that best describe the structure-activity relationship for this class of compounds. The best equation shows q
2
= 0.660, SE
cv
= 0.500, r
2
= 0.930, and SEE = 0.226. The external predictive ability of this best model was evaluated using a test set of 15 compounds. In order to design more potent DABO analogues as anti-HIV/AIDS agents, substituents capable of interactions with residues like
Ile94
,
Lys101
,
Tyr181
, and
Tyr188
should be selected. Also, given the importance of the conserved
Asn136
, this residue could become an attractive target for the design of novel NNRTIs with improved potency and increased ability to avoid the development of drug-resistant viruses.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.