Search

A series of tetrahydroisoquinolines functionalized with carbamates is reported here as highly selective ligands on the dopamine D2 receptor. These compounds were selected by means of a molecular modeling study. The studies were carried out in three stages: first an exploratory study was carried out using combined docking techniques and molecular dynamics simulations. According to these results, the bioassays were performed; these experimental studies corroborated the results obtained by molecular modeling. In the last stage of our study, a QTAIM analysis was performed in order to determine the main molecular interactions that stabilize the different ligand-receptor complexes. Our results show that the adequate use of combined simple techniques is a very useful tool to predict the potential affinity of new ligands at dopamine D1 and D2 receptors. In turn the QTAIM studies show that they are very useful to evaluate in detail the molecular interactions that stabilize the different ligand-receptor complexes; such information is crucial for the design of new ligands.

The inclusion complexes formed by chalcone and 2',4'-dihydroxychalcone with β-cyclodextrin have been studied combining experimental (phase solubility diagrams, Fourier transform infrared spectroscopy) and molecular modeling (molecular dynamics, quantum mechanics/molecular mechanics calculations) techniques. The formation constants of the complexes were determined at different temperatures, and the thermodynamic parameters of the process were obtained. The inclusion of chalcone in β-cyclodextrin is an exothermic process, while the inclusion of 2',4'-dihydroxychalcone is endothermic. Free energy profiles, derived from umbrella sampling using molecular dynamics simulations, were constructed to analyze the binding affinity and the complexation reaction at a molecular level. Hybrid QM/MM calculations were also employed to obtain a better description of the energetic and structural aspects of the complexes. The intermolecular interactions that stabilize both inclusion complexes were characterized by means of quantum atoms in molecules theory and reduce density gradient method. The calculated interactions were experimentally observed using FTIR.

Antimicrobial peptides form part of the immune system as protection against the action of external pathogens. The differences that exist between mammalian and microbial cell membrane architectures are key aspects of the ability of these peptides to discriminate between pathogens and host cells. Given that the pathogen membrane is the non-specific target of these cationic peptides, different molecular mechanisms have been suggested to describe the rules that permit them to distinguish between pathogens and mammalian cells. In this context, and setting aside the old fashion idea that cationic peptides act through one mechanism alone, this work will provide insight into the molecular action mechanism of small antimicrobial peptides, based on molecular dynamics simulations of phospholipid bilayers that mimic different cell membrane architectures. After measuring different properties of these lipid bilayers, in the absence and presence of peptides, a four-step action mechanism was suggested on the basis of the formation of phospholipid rafts induced by the presence of these cationic peptides. Thus, this work shows how differences in the bending modulus (k(b)) of these lipid rafts and differences in the free energy profiles (ΔG(z)) associated with the insertion of these peptides into these lipid rafts are key aspects for explaining the action mechanism of these cationic peptides at the molecular level.

Dopamine-mediated neurotransmission plays an important role in relevant psychiatric and neurological disorders. Nowadays, there is an enormous interest in the development of new dopamine receptors (DR) acting drugs as potential new targets for the treatment of schizophrenia or Parkinson's disease. Previous studies have revealed that isoquinoline compounds such as tetrahydroisoquinolines (THIQs) and tetrahydroprotoberberines (THPBs) can behave as selective D2 dopaminergic alkaloids since they share structural similarities with dopamine. In the present study we have synthesized eleven 2,3,9- and 2,3,11-trisubstituted THPB compounds (six of them are described for the first time) and evaluated their potential dopaminergic activity. Binding studies on rat striatal membranes were used to evaluate their affinity and selectivity towards D1 and D2 DR and establish the structure-activity relationship (SAR) as dopaminergic agents. In general, all the tested THPBs with protected phenolic hydroxyls showed a lower affinity for D1 and D2 DR than their corresponding homologues with free hydroxyl groups. In previous studies in which dopaminergic affinity of 1-benzyl-THIQs (BTHIQs) was evaluated, the presence of a Cl into the A-ring resulted in increased affinity and selectivity towards D2 DR. This is in contrast with the current study since the existence of a chlorine atom into the A-ring of the THPBs caused increased affinity for D1 DR but dramatically reduced the selectivity for D2 DR. An OH group in position 9 of the THPB (9f) resulted in a higher affinity for DR than its homologue with an OH group in position 11 (9e) (250 fold for D2 DR). None of the compounds showed any cytotoxicity in freshly isolated human neutrophils. A molecular modelling study of three representative THPBs was carried out. The combination of MD simulations with DFT calculations provided a clear picture of the ligand binding interactions from a structural and energetic point of view. Therefore, it is likely that compound 9d (2,3,9-trihydroxy-THPB) behave as D2 DR agonist since serine residues cluster are crucial for agonist binding and receptor activation., (Copyright © 2013 Elsevier Masson SAS. All rights reserved.)

A molecular modeling study on dihydrofolate reductase (DHFR) inhibitors was carried out. By combining molecular dynamics simulations with semiempirical (PM6), ab initio, and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of DHFR inhibitors interacting with the human enzyme is reported here, providing a clear picture of the binding interactions of these ligands from both structural and energetic viewpoints. A reduced model for the binding pocket was used. This approach allows us to perform more accurate quantum mechanical calculations as well as to obtain a detailed electronic analysis using the quantum theory of atoms in molecules (QTAIM) technique. Thus, molecular aspects of the binding interactions between inhibitors and the DHFR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental IC₅₀ values was obtained, predicting with an acceptable qualitative accuracy the potential inhibitor effect of nonsynthesized compounds. Such correlation was experimentally corroborated synthesizing and testing two new inhibitors reported in this paper.

A molecular modeling study on 16 1-benzyl tetrahydroisoquinolines (BTHIQs) acting as dopaminergic ligands was carried out. By combining molecular dynamics simulations with ab initio and density functional theory (DFT) calculations, a simple and generally applicable procedure to evaluate the binding energies of BTHIQs interacting with the human dopamine D2 receptor (D2 DR) is reported here, providing a clear picture of the binding interactions of BTHIQs from both structural and energetic viewpoints. Molecular aspects of the binding interactions between BTHIQs and the D2 DR are discussed in detail. A significant correlation between binding energies obtained from DFT calculations and experimental pKi values was obtained, predicting the potential dopaminergic effect of non-synthesized BTHIQs.

We report here an exhaustive and complete conformational study on the conformational potential energy hypersurface (PEHS) of dopamine (DA) interacting with the dopamine D2 receptor (D2-DR). A reduced 3D model for the binding pocket of the human D2-DR was constructed on the basis of the theoretical model structure of bacteriorhodopsin. In our reduced model system, only 13 amino acids were included to perform the quantum mechanics calculations. To obtain the different complexes of DA/D2-DR, we combined semiempirical (PM6), DFT (B3LYP/6-31G(d)), and QTAIM calculations. The molecular flexibility of DA interacting with the D2-DR was evaluated from potential energy surfaces and potential energy curves. A comparative study between the molecular flexibility of DA in the gas phase and at D2-DR was carried out. In addition, several molecular dynamics simulations were carried out to evaluate the molecular flexibility of the different complexes obtained. Our results allow us to postulate the complexes of type A as the "biologically relevant conformations" of DA. In addition, the theoretical calculations reported here suggested that a mechanistic stepwise process takes place for DA in which the protonated nitrogen group (in any conformation) acts as the anchoring portion, and this process is followed by a rapid rearrangement of the conformation allowing the interaction of the catecholic OH groups.