Lucas J. Gutierrez

Molecular modeling study of dihydrofolate reductase inhibitors. Molecular dynamics simulations, quantum mechanical calculations, and experimental corroboration.

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.

Structural and Thermodynamic Characteristics of the Exosite Binding Pocket on the Human BACE1: A Molecular Modeling Approach

We report a molecular modeling study aimed to locate and provide the full structural characteristics of the exosite binding site of the BACE1. A three-step procedure was followed. In the first stage, we performed blind docking studies on the whole target surface. In a second stage, the mode of binding was further refined by molecular dynamics (MD) simulation. Finally, binding free energy calculations, through the MM-PBSA protocol, were carried out to gain insight into the stability and thermodynamics of the inhibitor located at the selected binding pockets. Twelve binding pockets were identified on the surface of BACE1 by blind docking studies. The calculations of binding free energies for the 12 complexes show that van der Waals interactions dominate the mode of binding of these complexes. The best ranked complex shows that residues Glu255-Pro258, Phe261, Gly264-Ala272, Asp311-Ala313, Ser315, and Asp317-Tyr320 are located within 6 Å from the INH located at the exosite. The hydrogen bonds formed between the INH peptide, residues Tyr1, Tyr3, and Leu7 with the BACE1 residues Leu267, Cys269, Trp270, Asp311, and Asp 317 can strengthen the binding of the BACE1−INH complex.

Pentameric models as alternative molecular targets for the design of new antiaggregant agents

The structure-based drug design has been an extremely useful technique used for searching and developing of new therapeutic agents in various biological systems. In the case of AD, this approach has been difficult to implement. Among other several causes, the main problem might be the lack of a specific stable and reliable molecular target. In this paper the results obtained using a pentameric amyloid beta (Aβ) model as a molecular target are discussed. Our MD simulations have shown that this system is relatively structured and stable, displaying a lightly conformational flexibility during 2.0 μs of simulation time. This study allowed us to distinguish characteristic structural features in specific regions of the pentamer which should be taken into account when choosing this model as a molecular target. This represents a clear advantage compared to the monomer or dimer models which are highly flexible structures with large numbers of possible conformers. Using this pentameric model we performed two types of studies usually carried out on a molecular target: a virtual screening and the design on structural basis of new mimetic peptides with antiaggregant properties. Our results indicate that this pentameric model might be a good molecular target for these particular studies of molecular modeling. Details about the predictive power of our virtual screening as well as about the molecular interactions that stabilize the mimetic peptide-pentamer Aβ complexes are discussed in this paper.

Structural and functional insights into the anti-BACE1 Fab fragment that recognizes the BACE1 exosite

A molecular modeling study giving structural, functional, and mutagenesis insights into the anti-BACE1 Fab fragment that recognizes the BACE1 exosite is reported. Our results allow extending experimental data resulting from X-ray diffraction experiments in order to examine unknown aspects for the Fab-BACE1 recognition and its binding mode. Thus, the study performed here allows extending the inherently static nature of crystallographic structures in order to gain a deeper understanding of the structural and dynamical basis at the atomic level. The characteristics and strength of the interatomic interactions involved in the immune complex formation are exhaustively analyzed. The results might explain how the anti-BACE1 Fab fragment and other BACE1 exosite binders are capable to produce an allosteric modulation of the BACE1 activity. Our site-directed mutagenesis study indicated that the functional anti-BACE1 paratope, residues Tyr32 (H1), Trp50 (H2), Arg98 (H3), Phe101 (H3), Trp104 (H3) and Tyr94 (L3), strongly dominates the binding energetics with the BACE1 exosite. The mutational studies described in this work might accelerate the development of new BACE1 exosite binders with interesting pharmacological activity.

New substituted aminopyrimidine derivatives as BACE1 inhibitors: in silico design, synthesis and biological assays

We report in this work new substituted aminopyrimidine derivatives acting as inhibitors of the catalytic site of BACE1. These compounds were obtained from a molecular modeling study. The theoretical and experimental study reported here was carried out in several steps: docking analysis, Molecular Dynamics (MD) simulations, Quantum Theory Atom in Molecules (QTAIM) calculations, synthesis and bioassays and has allowed us to propose some compounds of this series as new inhibitors of the catalytic site of BACE1. The QTAIM study has allowed us to obtain an excellent correlation between the electronic densities and the experimental data of IC50. Also, using combined techniques (MD simulations and QTAIM calculations) enabled us to describe in detail the molecular interactions that stabilize the different L-R complexes. In addition, our results allowed us to determine what portion of these compounds should be changed in order to increase their affinity with the BACE1. Another interesting result is that a sort of synergism was observed when the effects of these new catalytic site inhibitors were combined with Ac-Tyr5-Pro6-Tyr7-Asp8-Ile9-Pro10-Leu11-NH2, which we have recently reported as a modulator of BACE1 acting on its exosite.

Tetrahydroisoquinolines functionalized with carbamates as selective ligands of D2 dopamine receptor

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.

Collection of alkenylcoumarin derivatives as Taq DNA polymerase inhibitors: SAR and in silico simulations

Using a feasible method, we generated a small focused library of structurally related alkenylcoumarins. These compounds were evaluated as potential antitumoral agents against Taq DNA polymerase. 6-(pent-4-enyloxy)-coumarin (7) IC50 = 48.33 ± 2.85 μM was defined as a small molecule able to disturb DNA replication. Docking and Molecular Dynamic Simulations suggest an active-site binding. Structure/activity relationship was reasonably established. Compound 7 represents a potential structure for further studies in the development of new anti-cancer DNA/polymerase binding agents.

Theoretical models to predict the inhibitory effect of ligands of sphingosine kinase 1 using QTAIM calculations and hydrogen bond dynamic propensity analysis

We report here the results of two theoretical models to predict the inhibitory effect of inhibitors of sphingosine kinase 1 that stand on different computational basis. The active site of SphK1 is a complex system and the ligands under the study possess a significant conformational flexibility; therefore for our study we performed extended simulations and proper clusterization process. The two theoretical approaches used here, hydrogen bond dynamics propensity analysis and Quantum Theory of Atoms in Molecules (QTAIM) calculations, exhibit excellent correlations with the experimental data. In the case of the hydrogen bond dynamics propensity analysis, it is remarkable that a rather simple methodology with low computational requirements yields results in excellent accord with experimental data. In turn QTAIM calculations are much more computational demanding and are also more complex and tedious for data analysis than the hydrogen bond dynamic propensity analysis. However, this greater computational effort is justified because the QTAIM study, in addition to giving an excellent correlation with the experimental data, also gives us valuable information about which parts or functional groups of the different ligands are those that should be replaced in order to improve the interactions and thereby to increase the affinity for SphK1. Our results indicate that both approaches can be very useful in order to predict the inhibiting effect of new compounds before they are synthesized.

Essential Dynamics on Different Biological Systems: Fis Protein, tvMyb1 Transcriptional Factor and BACE1 Enzyme

Proteins and enzymes poses a non-covalent 3D structure and therefore their intrinsic flexibility allows the existence of an ensemble of different conformers which are separated by low-energy barriers. These ranges of available conformers for proteins in solution are due to the relative movements among the different domains. Domain motions are important for a variety of protein functions, including catalysis, regulation of activity, transport of metabolites, formation of protein assemblies, and cellular locomotion.