Jans H. Alzate-Morales

Computational study on the interaction of N1 substituted pyrazole derivatives with B-raf kinase: an unusual water wire hydrogen-bond network and novel interactions at the entrance of the active site.

Docking and molecular dynamics (MD) simulations of N1 substituted pyrazole derivatives complexed with B-Raf kinase were performed to gain insight into the structural and energetic preferences of these inhibitors. First, a comparative study of fully automated docking programs AutoDock, ICM, GLIDE, and Surflex-Dock in closely approximating the X-ray crystal structure of the inhibitor (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime was performed. Afterward, the dynamics of the above-mentioned compound and the less active analogous compounds with 1-methyl-4-piperidinyl and tetrahydro-2H-pyran-4-yl groups at position N1 of pyrazole ring inside the B-Raf active site were analyzed by MD simulations. We found that the most active compound has stable interactions with residues Ile463 and His539 at the entrance of the B-Raf active site. Those interactions were in very good agreement with more reliable quantum mechanics/molecular mechanics calculations performed on the torsional angle phi between the pyrazole ring and the substituents at position N1. In addition, we identified a water wire connecting N2 of the pyrazole ring, Cys532, and Ser536, which is composed of three water molecules for the most active compound. We found some differences in the water wire hydrogen-bond network formed by less active compounds. We suggest that the differences between these structural features are responsible for the differences in activity among the studied compounds.

Docking and quantitative structure-activity relationship studies for 3-fluoro-4-(pyrrolo(2,1-f)(1,2,4)triazin-4-yloxy)aniline, 3-fluoro-4-(1H-pyrrolo(2,3-b)pyridin-4-yloxy)aniline, and 4-(4- amino-2-fluorophenoxy)-2-pyridinylamine derivatives as c-Met kinase inhibitors

We have performed docking of 3-fluoro-4-(pyrrolo[2,1-f][1,2,4]triazin-4-yloxy)aniline (FPTA), 3-fluoro-4-(1H-pyrrolo[2,3-b]pyridin-4-yloxy)aniline (FPPA), and 4-(4-amino-2-fluorophenoxy)-2-pyridinylamine (AFPP) derivatives complexed with c-Met kinase to study the orientations and preferred active conformations of these inhibitors. The study was conducted on a selected set of 103 compounds with variations both in structure and activity. Docking helped to analyze the molecular features which contribute to a high inhibitory activity for the studied compounds. In addition, the predicted biological activities of the c-Met kinase inhibitors, measured as IC50 values were obtained by using quantitative structure–activity relationship (QSAR) methods: Comparative molecular similarity analysis (CoMSIA) and multiple linear regression (MLR) with topological vectors. The best CoMSIA model included steric, electrostatic, hydrophobic, and hydrogen bond-donor fields; furthermore, we found a predictive model containing 2D-autocorrelation descriptors, GETAWAY descriptors (GETAWAY: Geometry, Topology and Atom-Weight AssemblY), fragment-based polar surface area (PSA), and MlogP. The statistical parameters: cross-validate correlation coefficient and the fitted correlation coefficient, validated the quality of the obtained predictive models for 76 compounds. Additionally, these models predicted adequately 25 compounds that were not included in the training set.

Synthesis, in silico, in vitro, and in vivo investigation of 5-[11C]methoxy-substituted sunitinib, a tyrosine kinase inhibitor of VEGFR-2

Sunitinib (SU11248) is a highly potent tyrosine kinase inhibitor targeting vascular endothelial growth factor receptor (VEGFR). Radiolabeled inhibitors of receptor tyrosine kinases (RTKs) might be useful tools for monitoring RTKs levels in tumor tissue giving valuable information for anti-angiogenic therapy. Herein we report the synthesis of 5-methoxy-sunitinib 5 and its (11)C-radiolabeled analog [(11)C]-5. The non-radioactive reference compound 5 was prepared by Knoevenagel condensation of 5-methoxy-2-oxindole with the corresponding substituted 5-formyl-1H-pyrrole. A binding constant (K(d)) of 20 nM for 5 was determined by competition binding assay against VEGFR-2. In addition, the binding mode of sunitinib and its 5-methoxy substituted derivative was studied by flexible docking simulations. These studies revealed that the substitution of the fluorine at position 5 of the oxindole scaffold by a methoxy group did not affect the inhibitor orientation, but affected the electrostatic and van der Waals interactions of the ligand with residues near the DFG motif of VEGFR-2. 5-[(11)C]methoxy-sunitinib ([(11)C]-5) was synthesized by reaction of the desmethyl precursor with [(11)C]CH(3)I in the presence of DMF and NaOH in 17 ± 3% decay-corrected radiochemical yield at a specific activity of 162-205 GBq/μmol (EOS). In vivo stability studies of [(11)C]-5 in rat blood showed that more than 70% of the injected compound was in blood stream, 60 min after administration.

Association of nicotinic acid with a poly(amidoamine) dendrimer studied by molecular dynamics simulations.

The interaction of poly(amidoamine)-G3 (PAMAM-G3) dendrimer with nicotinic acid (NA) was investigated by using molecular dynamics (MD) simulations. First, sample free energy profiles of NA crossing PAMAM-G3 at pH 6 and 3 were computed using the adaptive biasing force (ABF) method. We found that PAMAM-G3 provides a more appropriate environment for NA inclusion when internal tertiary amine groups are unprotonated (at pH 6). However, when internal tertiary amine groups are protonated (at pH 3), the PAMAM cavities are less hydrophobic; therefore the drug-dendrimer interactions become similar to drug-solvent interactions. Traditional MD simulations were also performed to investigate the structural stability of the PAMAM-NA complexes near the free energy minima at pH 6. We found that association of NA and PAMAM adopts a preferred binding mode around the surface of PAMAM, where hydrogen bond (HB) interactions with the amino and amide NH groups of the nearby monomers are established. These interactions are very stable whether additional van der Waals interactions between pyridine ring of NA and methylene groups of the more external monomers of PAMAM are established.

Molecular Dynamics of Protein Kinase-Inhibitor Complexes: A Valid Structural Information

Protein kinases (PKs) are key components of protein phosphorylation based signaling networks in eukaryotic cells. They have been identified as being implicated in many diseases. High-resolution X-ray crystallographic data exist for many PKs and, in many cases, these structures are co-complexed with inhibitors. Although this valuable information confirms the precise structure of PKs and their complexes, it ignores the dynamic movements of the structures which are relevant to explain the affinities and selectivity of the ligands, to characterize the thermodynamics of the solvated complexes, and to derive predictive models. Atomistic molecular dynamics (MD) simulations present a convenient way to study PK-inhibitor complexes and have been increasingly used in recent years in structure-based drug design. MD is a very useful computational method and a great counterpart for experimentalists, which helps them to derive important additional molecular information. That enables them to follow and understand structure and dynamics of protein-ligand systems with extreme molecular detail on scales where motion of individual atoms can be tracked. MD can be used to sample dynamic molecular processes, and can be complemented with more advanced computational methods (e.g., free energy calculations, structure-activity relationship analysis). This review focuses on the most commonly applications to study PK-inhibitor complexes using MD simulations. Our aim is that researchers working in the design of PK inhibitors be aware of the benefits of this powerful tool in the design of potent and selective PK inhibitors.

A Novel Class of Selective Acetylcholinesterase Inhibitors: Synthesis and Evaluation of (E)-2-(Benzo[d]thiazol-2-yl)-3-heteroarylacrylonitriles

(E)-2-(benzo[d]thiazol-2-yl)-3-heteroarylacrylonitriles are described as a new class of selective inhibitors of acetylcholinesterase (AChE). The most potent compound in the series exhibited good AChE inhibitory activity (IC50 = 64 µM). Compound 7f was found to be more selective than galanthamine in inhibiting AChE and it showed a moderate selectivity index. Kinetic studies on AChE indicated that a competitive type of inhibition pattern exist for these acrylonitrile derivates. Molecular docking models of the ligand-AChE complexes suggest that compound 7g is located on the periphery of the AChE active site.

Computational Study of the Interactions between Guanine Derivatives and Cyclin-Dependent Kinase 2 (CDK2) by CoMFA and QM/MM

Comparative molecular field analysis (CoMFA) and QM/MM hybrid calculations were performed on 9H-purine derivatives as CDK2 inhibitors. CoMFA was carried out to describe the activities of 78 analogues. The models were applied to a training set including 64 compounds. The best CoMFA model included steric and electrostatic fields, had a good Q(2) value of 0.845, and adequately predicted the compounds contained in the test set. Furthermore, plots of the steric CoMFA field allowed conclusions to be drawn for the choice of suitable inhibitors. In addition, the dynamical behavior of compounds with 4-(aminosulfonyl)phenyl, 4-[(methylamino)sulfonyl]phenyl, 4-[(dimethylamino)sulfonyl]phenyl, and [3-methoxy-4-(aminosulfonyl)]phenyl groups at position 2 of the 9H-purine scaffold inside the CDK2 active site were analyzed by QM/MM calculations. The interactions of these compounds with residues Lys89, Asp86, and Ile10 were characterized.

Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds.

N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.

The pH sensor of the plant K+-uptake channel KAT1 is built from a sensory cloud rather than from single key amino acids.

The uptake of potassium ions (K+) accompanied by an acidification of the apoplasm is a prerequisite for stomatal opening. The acidification (approximately 2-2.5 pH units) is perceived by voltage-gated inward potassium channels (K(in)) that then can open their pores with lower energy cost. The sensory units for extracellular pH in stomatal K(in) channels are proposed to be histidines exposed to the apoplasm. However, in the Arabidopsis thaliana stomatal K(in) channel KAT1, mutations in the unique histidine exposed to the solvent (His267) do not affect the pH dependency. We demonstrate in the present study that His267 of the KAT1 channel cannot sense pH changes since the neighbouring residue Phe266 shifts its pKa to undetectable values through a cation-π interaction. Instead, we show that Glu240 placed in the extracellular loop between transmembrane segments S5 and S6 is involved in the extracellular acid activation mechanism. Based on structural models we propose that this region may serve as a molecular link between the pH- and the voltage-sensor. Like Glu240, several other titratable residues could contribute to the pH-sensor of KAT1, interact with each other and even connect such residues far away from the voltage-sensor with the gating machinery of the channel.

Similarities between the Binding Sites of SB-206553 at Serotonin Type 2 and Alpha7 Acetylcholine Nicotinic Receptors: Rationale for Its Polypharmacological Profile

Evidence from systems biology indicates that promiscuous drugs, i.e. those that act simultaneously at various protein targets, are clinically better in terms of efficacy, than those that act in a more selective fashion. This has generated a new trend in drug development called polypharmacology. However, the rational design of promiscuous compounds is a difficult task, particularly when the drugs are aimed to act at receptors with diverse structure, function and endogenous ligand. In the present work, using docking and molecular dynamics methodologies, we established the most probable binding sites of SB-206553, a drug originally described as a competitive antagonist of serotonin type 2B/2C metabotropic receptors (5-HT2B/2CRs) and more recently as a positive allosteric modulator of the ionotropic α7 nicotinic acetylcholine receptor (nAChR). To this end, we employed the crystal structures of the 5-HT2BR and acetylcholine binding protein as templates to build homology models of the 5-HT2CR and α7 nAChR, respectively. Then, using a statistical algorithm, the similarity between these binding sites was determined. Our analysis showed that the most plausible binding sites for SB-206553 at 5-HT2Rs and α7 nAChR are remarkably similar, both in size and chemical nature of the amino acid residues lining these pockets, thus providing a rationale to explain its affinity towards both receptor types. Finally, using a computational tool for multiple binding site alignment, we determined a consensus binding site, which should be useful for the rational design of novel compounds acting simultaneously at these two types of highly different protein targets.