Mechanistic insights into the phosphoryl transfer reaction in cyclin-dependent kinase 2: A QM/MM study.

Cyclin-dependent kinase 2 (CDK2) is an important member of the CDK family exerting its most important function in the regulation of the cell cycle. It catalyzes the transfer of the gamma phosphate group from an ATP (adenosine triphosphate) molecule to a Serine/Threonine residue of a peptide substrate. Due to the importance of this enzyme, and protein kinases in general, a detailed understanding of the reaction mechanism is desired. Thus, in this work the phosphoryl transfer reaction catalyzed by CDK2 was revisited and studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our results suggest that the base-assisted mechanism is preferred over the substrate-assisted pathway when one Mg2+ is present in the active site, in agreement with a previous theoretical study. The base-assisted mechanism resulted to be dissociative, with a potential energy barrier of 14.3 kcal/mol, very close to the experimental derived value. An interesting feature of the mechanism is the proton transfer from Lys129 to the phosphoryl group at the second transition state, event that could be helping in neutralizing the charge on the phosphoryl group upon the absence of a second Mg2+ ion. Furthermore, important insights into the mechanisms in terms of bond order and charge analysis were provided. These descriptors helped to characterize the synchronicity of bond forming and breaking events, and to characterize charge transfer effects. Local interactions at the active site are key to modulate the charge distribution on the phosphoryl group and therefore alter its reactivity.

Experimental and Theoretical Approaches in the Study of Phenanthroline-Tetrahydroquinolines for Alzheimer's Disease.

The imino-Diels-Alder reaction is one of the most common strategies in organic chemistry and is an important tool for providing a broad spectrum of biologically active heterocyclic systems. A combined theoretical and experimental study of the imino-Diels-Alder reaction is described. The new phenanthroline-tetrahydroquinolines were evaluated as cholinesterase inhibitors. Their cytotoxicity in human neuroblastoma SH-SY5Y cells was also evaluated. The theoretical results suggest that compounds formation in stages can be explained by endo cycloadducts under the established reaction conditions, thereby confirming experimental results obtained for percentage yield. These results allowed us to establish that pyridine substituent remarkably influences activation energy and reaction yield, as well as in acetylcholinesterase (AChE) activity. Among these derivatives, compounds with 4-pyridyl and 4-nitrophenyl showed favorable AChE activity and proved to be non-cytotoxic.

Multiscale Simulations of Clavulanate Inhibition Identify the Reactive Complex in Class A ß-Lactamases and Predict the Efficiency of Inhibition.

Clavulanate is used as an effective drug in combination with ß-lactam antibiotics to treat infections of some antibiotic resistant bacteria. Here, we perform combined quantum mechanics/molecular mechanics simulations of several covalent complexes of clavulanate with class A ß-lactamases KPC-2 and TEM-1. Simulations of the deacylation reactions identify the decarboxylated trans-enamine complex as being responsible for inhibition. Further, the obtained free energy barriers discriminate clinically relevant inhibition (TEM-1) from less effective inhibition (KPC-2).

Design, facile synthesis, and evaluation of novel spiro- and pyrazolo[1,5-c]quinazolines as cholinesterase inhibitors: Molecular docking and MM/GBSA studies.

Given the wide spectrum of biological uses of pyrazolo[1,5-c]quinazoline and spiro-quinazoline derivatives as anticancer, anti-inflammatory analgesic agents, and their therapeutic applications in neurodegenerative disorders, it is compulsory to find easy, efficient, and simple methods to obtain and chemically diversify these families of compounds, thereby improving their biological applications. In this paper, we report the design and eco-friendly two-step synthesis of novel, fused spiro-pyrazolo[1,5-c]quinazoline derivatives as cholinesterase inhibitors. In addition, we studied their protein-ligand interactions via molecular docking and MM/GBSA calculations for a further rational design of more potent inhibitors. In first step, 2-(1H-pyrazol-5-yl)anilines were obtained through microwave (MW) assisted solvent-free/catalyst-free conditions and the second step involved the synthesis of the spiro-pyrazolo[1,5-c]quinazolines by a cyclocondensation reaction between 2-(1H-pyrazol-5-yl)anilines and cyclic ketones, or acetophenones, using stirring at room temperature. The compounds were obtained in high purity, good yields (50-97%), and at varying reaction times. The spiro-compounds were evaluated as acetylcholinesterase and butyrylcholinesterase inhibitors (AChEIs/BuChEIs) respectively, and the most potent compound exhibited a moderate AChE inhibitory activity (5f: IC50 = 84 µM). Molecular docking studies indicated that the binding mode of the compound 5f share common characteristics with the galantamine/donepezil-AChE complexes. Moreover, free binding energy (ΔG) calculations showed a good agreement with the experimental biological activity values. Our theoretical results indicated that halogen bond interactions could be involved with differential potency of these compounds and provide a new starting point to design novel pyrazolo[1,5-c]quinazolines as new anti-Alzheimer agents.

Combined molecular modelling and 3D-QSAR study for understanding the inhibition of NQO1 by heterocyclic quinone derivatives.

A combination of three-dimensional quantitative structure-activity relationship (3D-QSAR), and molecular modelling methods were used to understand the potent inhibitory NAD(P)H:quinone oxidoreductase 1 (NQO1) activity of a set of 52 heterocyclic quinones. Molecular docking results indicated that some favourable interactions of key amino acid residues at the binding site of NQO1 with these quinones would be responsible for an improvement of the NQO1 activity of these compounds. The main interactions involved are hydrogen bond of the amino group of residue Tyr128, π-stacking interactions with Phe106 and Phe178, and electrostatic interactions with flavin adenine dinucleotide (FADH) cofactor. Three models were prepared by 3D-QSAR analysis. The models derived from Model I and Model III, shown leave-one-out cross-validation correlation coefficients (q2LOO ) of .75 and .73 as well as conventional correlation coefficients (R2 ) of .93 and .95, respectively. In addition, the external predictive abilities of these models were evaluated using a test set, producing the predicted correlation coefficients (r2pred ) of .76 and .74, respectively. The good concordance between the docking results and 3D-QSAR contour maps provides helpful information about a rational modification of new molecules based in quinone scaffold, in order to design more potent NQO1 inhibitors, which would exhibit highly potent antitumor activity.

Structural and Affinity Determinants in the Interaction between Alcohol Acyltransferase from F. x ananassa and Several Alcohol Substrates: A Computational Study.

Aroma and flavor are important factors of fruit quality and consumer preference. The specific pattern of aroma is generated during ripening by the accumulation of volatiles compounds, which are mainly esters. Alcohol acyltransferase (AAT) (EC 2.3.1.84) catalyzes the esterification reaction of aliphatic and aromatic alcohols and acyl-CoA into esters in fruits and flowers. In Fragaria x ananassa, there are different volatiles compounds that are obtained from different alcohol precursors, where octanol and hexanol are the most abundant during fruit ripening. At present, there is not structural evidence about the mechanism used by the AAT to synthesize esters. Experimental data attribute the kinetic role of this enzyme to 2 amino acidic residues in a highly conserved motif (HXXXD) that is located in the middle of the protein. With the aim to understand the molecular and energetic aspects of volatiles compound production from F. x ananassa, we first studied the binding modes of a series of alcohols, and also different acyl-CoA substrates, in a molecular model of alcohol acyltransferase from Fragaria x ananassa (SAAT) using molecular docking. Afterwards, the dynamical behavior of both substrates, docked within the SAAT binding site, was studied using routine molecular dynamics (MD) simulations. In addition, in order to correlate the experimental and theoretical data obtained in our laboratories, binding free energy calculations were performed; which previous results suggested that octanol, followed by hexanol, presented the best affinity for SAAT. Finally, and concerning the SAAT molecular reaction mechanism, it is suggested from molecular dynamics simulations that the reaction mechanism may proceed through the formation of a ternary complex, in where the Histidine residue at the HXXXD motif deprotonates the alcohol substrates. Then, a nucleophilic attack occurs from alcohol charged oxygen atom to the carbon atom at carbonyl group of the acyl CoA. This mechanism is in agreement with previous results, obtained in our group, in alcohol acyltransferase from Vasconcellea pubescens (VpAAT1).

Performance of the MM/GBSA scoring using a binding site hydrogen bond network-based frame selection: the protein kinase case.

A conformational selection method, based on hydrogen bond (Hbond) network analysis, has been designed in order to rationalize the configurations sampled using molecular dynamics (MD), which are commonly used in the estimation of the relative binding free energy of ligands to macromolecules through the MM/GBSA or MM/PBSA method. This approach makes use of protein-ligand complexes obtained from X-ray crystallographic data, as well as from molecular docking calculations. The combination of several computational approaches, like long MD simulations on protein-ligand complexes, Hbond network-based selection by scripting techniques and finally MM/GBSA, provides better statistical correlations against experimental binding data than previous similar reported studies. This approach has been successfully applied in the ranking of several protein kinase inhibitors (CDK2, Aurora A and p38), which present both diverse and related chemical structures.

Medical care and organisation at the 2012 Roskilde Music Festival: a prospective observational study.

BACKGROUND: Music festivals, with a mix of music, alcohol and camping at the festival site combined with low hygiene, have become an integral part of society and attracts large numbers of guests. Our study was performed in order to increase the very limited knowledge on health-care issues and organisation of major music festivals. METHOD: Pre-defined illness and injury categories were constructed based on categories previously reported from music festivals. We prospectively recorded patient presentations to the Medical Health Care Organisation (MHCO) at the 2012 Roskilde Festival. RESULTS: During 10 days, more than 130,000 guests and volunteers attended the 2012 Roskilde Festival. Ten thousand six hundred thirty patient presentations were registered between the 30th of June and 9th of July 2012 by the MHCO. The majority of patient presentations, 6919, could be handled by first-aid volunteers with different levels of training. The remaining 3473 patient presentations were assessed to require further health-care competencies or additional resources such as prescriptions, medication or suturing. Two hundred thirty-eight patient presentations were triaged to a designated observation area. Two hundred sixty patients were referred to a local hospital, a general practitioner or a dentist. The overall patient presentation rate was 72/1000 attendees, and the transport-to-hospital rate was 1.8/1000 attendees. CONCLUSION: Our study demonstrates that illnesses and injuries are frequent, although mostly minor, in this normally low-risk population consisting of primarily young and healthy guests. However, comparison with other recent events was difficult as only limited data have been published from other music festivals. Future festivals should publish similar data.

Ultrasound-assisted phase-transfer catalysis method in an aqueous medium to promote the Knoevenagel reaction: advantages over the conventional and microwave-assisted solvent-free/catalyst-free method.

Given the broad spectrum of uses of acrylonitrile derivatives as fluorescent probes, AChE inhibitors, and others, it is necessary to find easy, efficient and simple methods to synthesize and diversify these compounds. We report the results of a comparative study of the effects of three techniques on the reactions between heterocyclic aldehydes and 2-(benzo[d]thiazol-2-yl)acetonitrile: stirring; ultrasound coupled to PTC conditions (US-PTC); and MW irradiation (MWI) under solvent and catalyst-free conditions. The effects of conditions on reaction parameters were evaluated and compared in terms of reaction time, yield, purity and outcomes. The US-PTC method is more efficient than the MWI and conventional methods. The reaction times were considerably shorter, with high yields (>90%) and good levels of purity. In addition, X-ray diffraction analysis and quantum mechanical calculations, at the level of density functional theory (DFT), ratify obtaining acrylonitrile isomers with E configurations. The crystal structure of 3c is stabilized by weak C-Ho⋯N intermolecular interactions (Ho⋯NC=2.45 Å, Co⋯NC=3.348(3) Å, Ho⋯NC=162°), forming centrosymmetric ring R2(2) (20) along the crystallographic a axis.

Analysis of intraprostatic therapeutic effects in prostate cancer patients using [(11)C]-choline pet/ct after external-beam radiation therapy.

PURPOSE: The objective of the present study was to analyze, with relatively high sensitivity and specificity, uptake properties of [(11)C]-choline in prostate cancer patients by means of positron-emission tomography (pet)/computed tomography (ct) imaging using objectively defined pet parameters to test for statistically significant changes before, during, and after external-beam radiation therapy (ebrt) and to identify the time points at which the changes occur. METHODS: The study enrolled 11 patients with intermediate-risk prostate cancer treated with ebrt, who were followed for up to 12 months after ebrt. The [(11)C]-choline pet scans were performed before treatment (baseline); at weeks 4 and 8 of ebrt; and at 1, 2, 3, 6, and 12 months after ebrt. RESULTS: Analysis of [(11)C]-choline uptake in prostate tissue before treatment resulted in a maximum standardized uptake value (suvmax) of 4.0 ± 0.4 (n = 11) at 40 minutes after injection. During week 8 of ebrt, the suvmax declined to 2.9 ± 0.1 (n = 10, p < 0.05). At 2 and 12 months after ebrt, suvmax values were 2.3 ± 0.3 (n = 10, p < 0.01) and 2.2 ± 0.2 (n = 11, p < 0.001) respectively, indicating that, after ebrt, maximum radiotracer uptake in the prostate was significantly reduced. Similar effects were observed when analyzing the tumour:muscle ratio (tmr). The tmr declined from 7.4 ± 0.6 (n = 11) before ebrt to 6.1 ± 0.4 (n = 11, nonsignificant) during week 8 of ebrt, to 5.6 ± 0.03 (n = 11, p < 0.05) at 2 months after ebrt, and to 4.4 ± 0.4 (n = 11, p < 0.001) at 12 months after ebrt. CONCLUSIONS: Our study demonstrated that intraprostatic [(11)C]-choline uptake in the 11 analyzed prostate cancer patients significantly declined during and after ebrt. The pet parameters SUVmax and tmr also declined significantly. These effects can be detected during radiation therapy and up to 1 year after therapy. The prognostic value of these early and statistically significant changes in intraprostatic [(11)C]-choline pet avidity during and after ebrt are not yet established. Future studies are indicated to correlate changes in [(11)C]-choline uptake parameters with long-term biochemical recurrence to further evaluate [(11)C]-choline pet changes as a possible, but currently unproven, biomarker of response.

Impact of a physician-staffed helicopter on a regional trauma system: a prospective, controlled, observational study.

INTRODUCTION: This study aims to compare the trauma system before and after implementing a physician-staffed helicopter emergency medical service (PS-HEMS). Our hypothesis was that PS-HEMS would reduce time from injury to definitive care for severely injured patients. METHODS: This was a prospective, controlled, observational study, involving seven local hospitals and one level I trauma centre using a before and after design. All patients treated by a trauma team within a 5-month period (1 December 2009-30 April 2010) prior to and a 12-month period (1 May 2010-30 April 2011) after implementing a PS-HEMS were included. We compared time from dispatch of the first ground ambulance to arrival in the trauma centre for patients with Injury Severity Score (ISS) > 15. Secondary end points were the proportion of secondary transfers and 30-day mortality. RESULTS: We included 1788 patients, of which 204 had an ISS > 15. The PS-HEMS transported 44 severely injured directly to the trauma centre resulting in a reduction of secondary transfers from 50% before to 34% after implementation (P = 0.04). Median delay for definitive care for severely injured patients was 218 min before and 90 min after implementation (P < 0.01). The 30-day mortality was reduced from 29% (16/56) before to 14% (21/147) after PS-HEMS (P = 0.02). Logistic regression showed PS-HEMS had an odds ratio (OR) for survival of 6.9 compared with ground transport. CONCLUSIONS: Implementation of a PS-HEMS was associated with significant reduction in time to the trauma centre for severely injured patients. We also observed significantly reduced proportions of secondary transfers and 30-day mortality.

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.

Properties of noise in positron emission tomography images reconstructed with filtered-backprojection and row-action maximum likelihood algorithm.

Noise levels observed in positron emission tomography (PET) images complicate their geometric interpretation. Post-processing techniques aimed at noise reduction may be employed to overcome this problem. The detailed characteristics of the noise affecting PET images are, however, often not well known. Typically, it is assumed that overall the noise may be characterized as Gaussian. Other PET-imaging-related studies have been specifically aimed at the reduction of noise represented by a Poisson or mixed Poisson + Gaussian model. The effectiveness of any approach to noise reduction greatly depends on a proper quantification of the characteristics of the noise present. This work examines the statistical properties of noise in PET images acquired with a GEMINI PET/CT scanner. Noise measurements have been performed with a cylindrical phantom injected with (11)C and well mixed to provide a uniform activity distribution. Images were acquired using standard clinical protocols and reconstructed with filtered-backprojection (FBP) and row-action maximum likelihood algorithm (RAMLA). Statistical properties of the acquired data were evaluated and compared to five noise models (Poisson, normal, negative binomial, log-normal, and gamma). Histograms of the experimental data were used to calculate cumulative distribution functions and produce maximum likelihood estimates for the parameters of the model distributions. Results obtained confirm the poor representation of both RAMLA- and FBP-reconstructed PET data by the Poisson distribution. We demonstrate that the noise in RAMLA-reconstructed PET images is very well characterized by gamma distribution followed closely by normal distribution, while FBP produces comparable conformity with both normal and gamma statistics.

Synthesis, in silico, in vitro, and in vivo investigation of 5-[¹¹C]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.

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 7 g is located on the periphery of the AChE active site.

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.

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.

Study of differences in the VEGFR2 inhibitory activities between semaxanib and SU5205 using 3D-QSAR, docking, and molecular dynamics simulations.

Semaxanib (SU5416) and 3-[4'-fluorobenzylidene]indolin-2-one (SU5205) are structurally similar drugs that are able to inhibit vascular endothelial growth factor receptor-2 (VEGFR2), but the former is 87 times more effective than the latter. Previously, SU5205 was used as a radiolabelled inhibitor (as surrogate for SU5416) and a radiotracer for positron emission tomography (PET) imaging, but the compound exhibited poor stability and only a moderate IC(50) toward VEGFR2. In the current work, the relationship between the structure and activity of these drugs as VEGFR2 inhibitors was studied using 3D-QSAR, docking and molecular dynamics (MD) simulations. First, comparative molecular field analysis (CoMFA) was performed using 48 2-indolinone derivatives and their VEGFR2 inhibitory activities. The best CoMFA model was carried out over a training set including 40 compounds, and it included steric and electrostatic fields. In addition, this model gave satisfactory cross-validation results and adequately predicted 8 compounds contained in the test set. The plots of the CoMFA fields could explain the structural differences between semaxanib and SU5205. Docking and molecular dynamics simulations showed that both molecules have the same orientation and dynamics inside the VEGFR2 active site. However, the hydrophobic pocket of VEGFR2 was more exposed to the solvent media when it was complexed with SU5205. An energetic analysis, including Embrace and MM-GBSA calculations, revealed that the potency of ligand binding is governed by van der Waals contacts.

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.

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 IC(50) 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.