Behavior of uncharged oximes compared to HI6 and 2-PAM in the human AChE-tabun conjugate: a molecular modeling approach.

Tabun is one of the most dangerous nerve agents because it has deleterious effects like inhibition of the essential enzymes acetylcholinesterase (AChE) and butyrylcholinesterase. Some oximes such HI6 as 2-PAM are nucleophiles that are capable to reactivate inhibited human AChE under some conditions. Zwitterionic and cationic species have the best chance of productive action on inhibited AChE. However uncharged oximes can give important interaction information. In order to investigate the interaction and behavior of cationic and uncharged oximes, we performed molecular docking simulations and molecular dynamics and calculated binding energies of complexes of these compounds with human AChE. The uncharged oximes of larger structure were more susceptible to the influence of the substituents on the phosphorus atom and presented low binding energies. In contrast, HI 6 and 2-PAM showed high binding energy values with great contribution of the amino acid Asp74, demonstrating the importance of the quaternary nitrogen to the affinity and interaction of the oximes/AChE tabun-inhibited complexes.

Dynamical study, hydrogen bond analysis, and constant pH simulations of the beta carbonic anhydrase of Methanobacterium thermoautotrophicum.

Within the five classes (α, ß, γ, δ, and ζ) of carbonic anhydrases (CAs) the first two, containing mammal and plant representatives, are the most studied among all CAs. In this study, we have focused our investigation on the beta-class carbonic anhydrase of Methanobacterium thermoautotrophicum. We investigated both the importance of the Asp-Arg dyad near the catalytic zinc-bound water and the possible roles that water molecules within the active site and residues near the entrance of the catalytic cleft have on the first step of the enzyme's reaction mechanism. Hydrogen-bonding analysis of selected residues within the active site and constant pH replica exchange molecular dynamics constant pH replica exchange simulations were performed. The latter was done in order to evaluate the pKa values of possible proton acceptors. We found an intricate hydrogen-bonding network involving two acidic residues within the active site, Asp16 and Asp34, and the catalytic water molecule. We also observed a very strong interaction between the zinc-bound water and residues Asp34 and Arg36. This interaction was not significantly affected by the change in the protonation state of both the catalytic water and aspartate residue 34. The pKa analysis show that the effect of the R36A mutation affects not only the possible proton acceptors, but also the catalytic water itself.

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors.

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20-40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (ELJ) and electrostatic (EC) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pKi=7.772-0.1885×Gly96+0.0517×Tyr158 (R²=0.80; n=10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.

Hologram QSAR models of a series of 6-arylquinazolin-4-amine inhibitors of a new Alzheimer's disease target: dual specificity tyrosine-phosphorylation-regulated kinase-1A enzyme.

Dual specificity tyrosine-phosphorylation-regulated kinase-1A (DYRK1A) is an enzyme directly involved in Alzheimer's disease, since its increased expression leads to ß-amyloidosis, Tau protein aggregation, and subsequent formation of neurofibrillary tangles. Hologram quantitative structure-activity relationship (HQSAR, 2D fragment-based) models were developed for a series of 6-arylquinazolin-4-amine inhibitors (36 training, 10 test) of DYRK1A. The best HQSAR model (q2 = 0.757; SEcv = 0.493; R2 = 0.937; SE = 0.251; R2pred = 0.659) presents high goodness-of-fit (R2 > 0.9), as well as high internal (q2 > 0.7) and external (R2pred > 0.5) predictive power. The fragments that increase and decrease the biological activity values were addressed using the colored atomic contribution maps provided by the method. The HQSAR contribution map of the best model is an important tool to understand the activity profiles of new derivatives and may provide information for further design of novel DYRK1A inhibitors.

A versatile low temperature solid-state synthesis of vanadium nitride (VN) via a "guanidinium-route": experimental and theoretical studies from the key-intermediate to the final product.

An easy and efficient method ("guanidinium route") to synthesize vanadium nitride (VN) is evaluated in this paper. Initially, ammonium m-vanadate was mixed with guanidinium carbonate, producing an important intermediate, ammonium m-vanadate (GmV), through a solid-state reaction. GmV was decomposed as a function of the temperature and studied with TGA, DRX, FT-Infrared, Temperature-Programmed Surface Reaction (TPSR) and DFT. We show that GmV is transformed into bulk crystalline VN below 800 °C. XPS, periodic DFT calculations, and elemental analyses show that the surface is not carbon-free.

Application of 4D-QSAR studies to a series of raloxifene analogs and design of potential selective estrogen receptor modulators.

Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis was applied on a series of 54 2-arylbenzothiophene derivatives, synthesized by Grese and coworkers, based on raloxifene (an estrogen receptor-alpha antagonist), and evaluated as ERa ligands and as inhibitors of estrogen-stimulated proliferation of MCF-7 breast cancer cells. The conformations of each analogue, sampled from a molecular dynamics simulation, were placed in a grid cell lattice according to three trial alignments, considering two grid cell sizes (1.0 and 2.0 Å). The QSAR equations, generated by a combined scheme of genetic algorithms (GA) and partial least squares (PLS) regression, were evaluated by "leave-one-out" cross-validation, using a training set of 41 compounds. External validation was performed using a test set of 13 compounds. The obtained 4D-QSAR models are in agreement with the proposed mechanism of action for raloxifene. This study allowed a quantitative prediction of compounds' potency and supported the design of new raloxifene analogs.

Residue-ligand interaction energy (ReLIE) on a receptor-dependent 3D-QSAR analysis of S- and NH-DABOs as non-nucleoside reverse transcriptase inhibitors.

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² = 0.660, SE(cv) = 0.500, r² = 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.

Docking of the alkaloid geissospermine into acetylcholinesterase: a natural scaffold targeting the treatment of Alzheimer's disease.

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.

Receptor-dependent (RD) 3D-QSAR approach of a series of benzylpiperidine inhibitors of human acetylcholinesterase (HuAChE).

Acetylcholine inhibitors (AChEIs) are currently considered as potential drugs for treating Alzheimer disease. In this work, we developed a receptor-dependent 3D-QSAR (RD-3D-QSAR) models based on a series of 60 benzylpiperidine inhibitors of human acetylcholinesterase to support the design of new AChEIs. The best two models, A-F (N = 47, q(2) = 0.736, r(2) = 0.860) and C-F (N = 47, q(2) = 0.753, r(2) = =0.900) were developed and validated by a combined GA-PLS approach, available in WOLF. Residues of the aromatic gorge (Tyr341 and Trp439) and catalytic triad (His447) are related to both equations showing the consistency of these models with the SAR. Based on those models we have proposed four new benzylpiperidine derivatives and predicted the pIC(50) for each molecule. The good predicted potency of benzylpiperidine derivative, IIa, indicates that it is a potential candidate as a new HuAChE inhibitor.

Investigating the differential activation of vascular endothelial growth factor (VEGF) receptors.

The vascular endothelial growth factors are key mediators of angiogenesis and are also related to several physiological processes such as monocyte chemotaxis, dendritic cell development, hematopoietic stem cell survival, and many others. PlGF, VEGF, VEGFB, VEGFC and VEGFD were identified as members of the vascular endothelial growth factor family. They act by differential activation of three receptors: Flt-1, KDR and Flt-4. PlGF and VEGFB only activate Flt-1. VEGF activates both Flt-1 and KDR. VEGFC and VEGFD activate KDR and Flt-4. The available three dimensional structures of VEGF and PlGF, in complex with the domain-2 of Flt-1, show that both proteins bind in a very similar way to Flt-1 receptor. Here we construct the three dimensional model of the domain-2 of KDR receptor using the same domain of Flt-1 as template. We also construct the model complexes VEGF/KDR, VEGFB/Flt-1, VEGFB/KDR and PlGF/KDR. Molecular dynamics simulations with explicit solvent are carried out on eleven molecular systems: unbound VEGF, VEGF/Flt-1(D2), VEGF/KDR(D2), unbound PlGF, PlGF/Flt-1(D2), PlGF/KDR(D2), unbound VEGFB, VEGFB/Flt-1(D2), VEGFB/KDR(D2), unbound Flt-1(D2) and unbound KDR(D2). We analyze protein-protein interactions, shape complementarity, charge complementarity and hydrogen bonds. As a coarse estimation of the desolvation penalties, we assume a correlation to the number of hydrogen bonds with solvent molecules that are lost upon complex formation. The results herein are consistent with the experimental selectivity profile (VEGF being able to activate both Flt-1 and KDR receptors while VEGFB and PlGF being only able to activate Flt-1), and provide a collection of evidences sustaining the complementarity of polar interactions as the main responsible for protein recognition and selectivity.

Molecular docking of a series of peptidomimetics in the trypanothione binding site of T. cruzi trypanothione reductase.

Chagas' disease (CD) has been responsible for many deaths and disabilities mainly in South America. Currently, 40 million people are at risk of acquiring this disease and, existing therapies are still unsatisfactory, presenting harsh side effects. Therefore, the development of new chemical entities to reverse this state is critical. A series of peptidomimetics, developed by Mc Kie et al. (2001) [11], showed a reversible and competitive inhibition against Trypanosoma cruzi Trypanothione Reductase (TR). These inhibitors may be used as basis of lead compounds in the design of new drug candidates for the treatment of CD. In this work, we have docked this series of peptidomimetics into the TR binding site, using the FlexX algorithm as implemented in the Sybyl program, in order to access the binding mode of this class of compounds in the target enzyme.

The receptor-dependent QSAR paradigm: an overview of the current state of the art.

The original quantitative structure-activity relationship (QSAR) formulation was proposed by Hansch and Fujita in the 1960's and, since then QSAR analysis has evolved as a mature science, due mainly to the advances that occurred in the past two decades in the fields of molecular modeling, data analysis algorithms, chemoinformatics, and the application of graph theory in chemistry. Moreover, it is also worthy of note the exponential progress that have occurred in software and hardware development. In this context, a myriad of QSAR methods exist; from the considered "classical" approaches (known as two-dimensional (2D) QSAR), to three-dimensional (3D) and multidimensional (nD) QSAR ones. A distinct QSAR approach has been recently proposed, the receptor-dependent-QSAR, where explicit information regarding the receptor structure (usually a protein) is extensively used during modeling process. Indeed, a limited, but growing number of receptor-dependent QSAR methods are reported in the literature. With no intention to be comprehensive, an overview of receptor-dependent QSAR methods will be discussed along with an in-depth examination of their applications in drug design, virtual screen, and ADMET modeling in silico.

3D-QSAR CoMFA/CoMSIA models based on theoretical active conformers of HOE/BAY-793 analogs derived from HIV-1 protease inhibitor complexes.

The three-dimensional quantitative structure-activity relationships (3D-QSAR) of a series of HOE/BAY-793 analogs (C(2)-symmetric diol peptidomimetics), developed by Budt and co-workers [Bioorg. Med. Chem. 3 (1995) 559] as inhibitors of HIV-1 protease (HIV-PR), were studied using Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA). Theoretical active conformers for these peptidomimetics were generated, derived from modeled protease inhibitor complexes, in order to orient the compounds superposition and to afford a consistent alignment. The best CoMFA model (N=27, q(2)=0.637, R(2)=0.991) showed contributions of the steric (45.7%) and electrostatic (54.3%) fields to the activity, while the best CoMSIA model (N=27, q(2)=0.511, R(2)=0.987) showed contributions of the electrostatic (68.5%) and hydrogen bond donor (37.5%) fields. The models were also external validated using four compounds (test set) not included in the model generation process. The statistical parameters from both models indicate that the data are well fitted and have high predictive ability. Moreover, the resulting 3D CoMFA/CoMSIA contour maps provide useful guidance for designing highly active ligands. The CoMFA/CoMSIA models were also compared with previous 4D-QSAR models [E.F.F. da Cunha, M.G. Albuquerque, O.A.C. Antunes, R.B. de Alencastro, QSAR Comb. Sci. 24 (2005), 240-253.].

3D-QSAR CoMFA of a series of DABO derivatives as HIV-1 reverse transcriptase non-nucleoside inhibitors.

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 comparative molecular field analysis (CoMFA) in order to derive three-dimensional quantitative structure-activity relationship (3D-QSAR) models. The CoMFA study has been performed with a training set of 59 compounds, testing three alignments and four charge schemes (DFT, HF, AM1, and PM3) and using defaults probe atom (Csp (3), +1 charge), cutoffs (30 kcal.mol (-1) for both steric and electrostatic fields), and grid distance (2.0 A). The best model ( N = 59), derived from Alignment 1 and PM3 charges, shows q (2) = 0.691, SE cv = 0.475, optimum number of components = 6, r (2) = 0.930, SEE = 0.226, and F-value = 115.544. The steric and electrostatic contributions for the best model were 43.2% and 56.8%, respectively. The external predictive ability (r (2) pred = 0.918) of the resultant best model was evaluated using a test set of 15 compounds. In order to design more potent DABO analogues as anti-HIV/AIDS agents, attention should be taken in order to select a substituent for the 4-oxopyrimidine ring, since, as revealed by the best CoMFA model, there are a steric restriction at the C2-position, a electron-rich group restriction at the C6-position ( para-substituent of the 6-benzyl group), and a steric allowed region at the C5-position.

On the structure, interactions, and dynamics of bound VEGF.

The vascular endothelial growth factor (VEGF) is believed to be the most important protein in the regulation of the angiogenic cascade. Thus, exploring the structure and dynamical properties of this growth factor and the influence of receptor and inhibitor binding to these properties may reveal new insights on VEGF's biological process and inhibition opportunities. Here we describe an analysis of molecular dynamics simulations of VEGF bound to the Flt-1 receptor, VEGF bound to the v107 peptide inhibitor, and also VEGF bound to a mutant v107. We analyze the effects of binding to VEGF regarding three aspects: structure, interactions, and dynamics. We found that the structure of VEGF is not significantly perturbed upon binding. We analyze the individual contribution of the VEGF residues to the total interaction energy of binding to Flt-1 and v107. We also compare dynamical variables such as thermal fluctuations and correlations with those of the unbound form. We found that receptor binding is able to promote stronger perturbations on the VEGF dynamical behavior than VEGF inhibitor binding. VEGF motions in the receptor bound complex are shown to be less correlated than motions of unbound VEGF. The work addresses the changes on conformational flexibility of the isolated VEGF upon binding, as well as changes in structure and side-chain rearrangements.

Human Cytomegalovirus Protease: Why is the Dimer Required for Catalytic Activity?

Human cytomegalovirus (HCMV) is a pathogenic agent responsible for morbidity and mortality in immunocompromised and immunosuppressed individuals. HCMV encodes a serine protease that is essential for the production of infectious virions. In this work, we applied molecular dynamics (MD) simulations on HCMV protease models in order to investigate the experimentally observed (i) catalytic activity of the enzyme homodimer and (ii) induced-fit mechanism upon the binding of substrates and peptidyl inhibitors. Long and stable trajectories were obtained for models of the monomeric and dimeric states, free in solution and bound covalently and noncovalently to a peptidyl-activated carbonyl inhibitor, with very good agreement between theoretical and experimental results. The MD results suggest that HCMV protease indeed operates by an induced-fit mechanism. Also, our analysis indicates that the catalytic activity of the dimer is a result of more favorable interactions between the oxyanion in the covalently bound state and the backbone nitrogen of Arg165, resulting in a reaction that is 7.0 kcal/mol more exergonic and a more significant thermodynamic driving force. The incipient oxyanion in the transition state should also benefit from the stronger interactions with Arg165, reducing in this manner the intrinsic activation barrier for the reaction in the dimeric state.

Free-energy force-field three-dimensional quantitative structure-activity relationship analysis of a set of p38-mitogen activated protein kinase inhibitors.

The p38-mitogen-activated protein kinases (p38-MAPKs) belong to a family of serine-threonine kinases activated by pro-inflammatory or stressful stimuli that are known to be involved in several diseases. Their biological importance, related to the release of inflammatory pro-cytokines such as tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1), has generated many studies aiming at the development of selective inhibitors for the treatment of inflammatory diseases. In this work, we developed receptor-based three dimensional (3D) quantitative structure-activity relationship (QSAR) models for a series of 33 pyridinyl imidazole compounds [Liverton et al. (1999) 42:2180], using a methodology named free-energy force-field (FEFF) [Tokarski and Hopfinger (1997) 37:792], in which scaled intra- and intermolecular energy terms of the Assisted Model Building Energy Refinement (AMBER) force field combined with a hydration-shell solvation model are the independent variables used in the QSAR studies. Multiple temperature molecular-dynamics simulations (MDS) of ligand-protein complexes and genetic-function approximation (GFA) were employed using partial least squares (PLS) as the fitting functions to develop FEFF-3D-QSAR models for the binding process. The best model obtained in the FEFF-3D-QSAR receptor-dependent (RD) method shows the importance of the van der Waals energy change upon binding and the electrostatic energy in the interaction of ligands with the receptor. The QSAR equations described here show good predictability and may be regarded as representatives of the binding process of ligands to p38-MAPK. Additionally, we have compared the top FEFF-3D-QSAR model with receptor independent (RI) 4D-QSAR models developed in a recent study [Romeiro et al. (2005) 19:385].

Construction of 4D-QSAR models for use in the design of novel p38-MAPK inhibitors.

The p38-mitogen-activated protein kinase (p38-MAPK) plays a key role in lipopolysaccharide-induced tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1) release during the inflammatory process, emerging as an attractive target for new anti-inflammatory agents. Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis [Hopfinger et al., J. Am. Chem. Soc., 119 (1997) 10509] was applied to a series of 33 (a training set of 28 and a test set of 5) pyridinyl-imidazole and pyrimidinyl-imidazole inhibitors of p38-MAPK, with IC50 ranging from 0.11 to 2100 nM [Liverton et al., J. Med. Chem., 42 (1999) 2180]. Five thousand conformations of each analogue were sampled from a molecular dynamics simulation (MDS) during 50 ps at a constant temperature of 303 K. Each conformation was placed in a 2 angstroms grid cell lattice for each of three trial alignments. 4D-QSAR models were constructed by genetic algorithm (GA) optimization and partial least squares (PLS) fitting, and evaluated by leave-one-out cross-validation technique. In the best models, with three to six terms, the adjusted cross-validated squared correlation coefficients, Q2adj, ranged from 0.67 to 0.85. Model D (Q2adj = 0.84) was identified as the most robust model from alignment 1, and it is representative of the other best models. This model encompasses new molecular regions as containing pharmacophore sites, such as the amino-benzyl moiety of pyrimidine analogs and the N1-substituent in the imidazole ring. These regions of the ligands should be further explored to identify better anti-inflammatory inhibitors of p38-MAPK.

Theoretical evaluation of adiabatic and vertical electron affinity of some radiosensitizers in solution using FEP, ab initio and DFT methods.

The biological activity of radiosensitizers is associated to their electron affinity (EA), which can be divided in two main processes: vertical and adiabatic. In this work, we calculated the EAs of nitrofurans and nitroimidazoles (Fig. 2) using Hartree-Fock (HF) and density functional theory (DFT) methods and evaluated solvent effects (water and carbon tetrachloride) on EAs. For water, we combined the polarized continuum model (PCM) and free energy perturbation (FEP) (finite difference thermodynamic integration, FDTI) methods. For carbon tetrachloride, we used the FDTI method. The values of adiabatic EA obtained are in agreement with experimental data (deviations of 0.013 eV). The vertical EAs were calculated according to Cederbaum's outer valence Green function (OVGF) method. This methodology, which relies on theoretical aspects of free energy calculations on charged molecules in solution, was used to select potential selective radiosensitizers from recently reported compounds and could be helpful in the rational design of new and more selective bioreductive anticancer drugs.