Marawan Ahmed

Discovery of new HER2/EGFR dual kinase inhibitors based on the anilinoquinazoline scaffold as potential anti-cancer agents

Abstract Herein, we designed and synthesized certain anilinoquinazoline derivatives bearing bulky arylpyridinyl, arylpropenoyl and arylpyrazolyl moieties at the 4′ position of the anilinoquinazoline, as potential dual HER2/EGFR kinase inhibitors. A detailed molecular modeling study was performed by docking the synthesized compounds in the active site of the epidermal growth factor receptor (EGFR). The synthesized compounds were further tested for their inhibitory activity on EGFR and HER2 tyrosine kinases. The aryl 2-imino-1,2-dihydropyridine derivatives 5d and 5e displayed the most potent inhibitory activity on EGFR with IC50 equal to 2.09 and 1.94 μM, respectively, and with IC50 equal to 3.98 and 1.04 μM on HER2, respectively. Furthermore, the anti-proliferative activity of these most active compounds on MDA-MB-231 breast cancer cell lines, known to overexpress EGFR, showed an IC50 range of 2.4 and 2.5 μM, respectively.

Assessment of new anti-HER2 ligands using combined docking, QM/MM scoring and MD simulation

In the development of new anti-cancer drugs to tackle the problem of resistance to current chemotherapeutic agents, a new series of anti-HER2 (human epidermal growth factor receptors 2) agents has been synthesized and investigated using different computational methods. Although non-selective, the most active inhibitor in the new series shows higher activity toward HER2 than EGFR. The induced fit docking protocol (IFD) is performed to find possible binding poses of the new inhibitors in the active site of the HER2 receptor. Molecular dynamic simulations of the inhibitor-protein complexes for the two most active compounds from the new series are carried out. Simulations stability is checked using different stability parameters. Different scoring functions are employed.

In silico design: Extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2

Based on the hit structures that have been identified in our previous studies against EGFR and HER2, new potential inhibitors that share the same scaffold of the hit structures are designed and screened in silico. Insights into understanding the potential inhibitory effect of the new inhibitors against both EGFR and HER2 receptors is obtained using extended molecular dynamics (MD) simulations and different scoring techniques. The binding mechanisms and dynamics are detailed with respect to two approved inhibitors against EGFR (lapatinib) and HER2 (SYR127063). The best scoring inhibitor (T9) is chosen for additional in silico investigation against both the wild-type and T790M mutant strain of EGFR and the wild-type HER2. The results reveal that certain substitution patterns increase the stability and assure stronger binding and higher H-bond occupancy of the conserved water molecule that is commonly observed with kinase crystal structures. Furthermore, the new inhibitor (T9) forms stable interactions with the mutant strain as a direct consequence of the enhanced ability to form additional hydrogen bonding interactions with binding site residues.

Photoelectron spectra of some antibiotic building blocks: 2-azetidinone and thiazolidine-carboxylic acid.

X-ray photoelectron spectra of the core and valence levels of the fundamental building blocks of β-lactam antibiotics have been investigated and compared with theoretical calculations. The spectra of the compounds 2-azetidinone and the 2- and 4-isomers of thiazolidine-carboxylic acid are interpreted in the light of theoretical calculations. The spectra of the two isomers of thiazolidine-carboxylic acid are rather similar, as expected, but show clear effects due to isomerization. Both isomers are analogues of proline, which is well-known to populate several low energy conformers in the gas phase. We have investigated the low energy conformers of thiazolidine-4-carboxylic acid theoretically in more detail and find some spectroscopic evidence that multiple conformers may be present. The measured valence levels are assigned for all three compounds, and the character of the frontier orbitals is identified and analyzed.

Adsorption of cytosine and AZA derivatives of cytidine on Au single crystal surfaces

The adsorption of cytosine on the Au(111) and Au(110) surfaces has been studied using both aqueous deposition and evaporation in vacuum to prepare the samples. Soft X-ray photoelectron spectroscopy (XPS) and near edge X-ray absorption fine structure spectroscopy (NEXAFS) were used to determine the electronic structure and orientation of the adsorbates. In addition, three derivatives of cytosine, 6-azacytosine, 6-azacytidine and 5- azacytidine, were studied. Monolayer films of the latter three samples were adsorbed on Au(111) from aqueous solution, and the nature of bonding was determined. Spectra have been interpreted in the light of published calculations of free cytosine molecules and new ab initio calculations of the other compounds. Surface core level shifts of Au 4f imply that all of these compounds are chemisorbed. Cytosine adsorbs as a single tautomer, but in two chemical states with different surface-molecule bonding. For deposition in vacuum, a flat-lying molecular state bonded through the N(3) atom of the pyrimidine ring dominates, but a second state is also present. For deposition from solution, the second state dominates, with the molecular plane no longer parallel to the surface. This state also bonds through the N(3) atom, but in addition interacts with the surface via the amino group. Two tautomers of 6-azacytosine were observed, and they and 6-azacytidine adsorb with similar geometries, chemically bonding via the azacytosine ring. The ribose ring does not appear to perturb the adsorption of azacytidine compared with azacytosine. The azacytosine ring is nearly but not perfectly parallel to the surface, like 5-azacytidine, which adsorbs as an imino tautomer. ...

Modeling the human Na v 1.5 sodium channel: structural and mechanistic insights of ion permeation and drug blockade

Abnormalities in the human Nav1.5 (hNav1.5) voltage-gated sodium ion channel (VGSC) are associated with a wide range of cardiac problems and diseases in humans. Current structural models of hNav1.5 are still far from complete and, consequently, their ability to study atomistic interactions of this channel is very limited. Here, we report a comprehensive atomistic model of the hNav1.5 ion channel, constructed using homology modeling technique and refined through long molecular dynamics simulations (680 ns) in the lipid membrane bilayer. Our model was comprehensively validated by using reported mutagenesis data, comparisons with previous models, and binding to a panel of known hNav1.5 blockers. The relatively long classical MD simulation was sufficient to observe a natural sodium permeation event across the channel’s selectivity filters to reach the channel’s central cavity, together with the identification of a unique role of the lysine residue. Electrostatic potential calculations revealed the existence of two potential binding sites for the sodium ion at the outer selectivity filters. To obtain further mechanistic insight into the permeation event from the central cavity to the intracellular region of the channel, we further employed “state-of-the-art” steered molecular dynamics (SMD) simulations. Our SMD simulations revealed two different pathways through which a sodium ion can be expelled from the channel. Further, the SMD simulations identified the key residues that are likely to control these processes. Finally, we discuss the potential binding modes of a panel of known hNav1.5 blockers to our structural model of hNav1.5. We believe that the data presented here will enhance our understanding of the structure–property relationships of the hNav1.5 ion channel and the underlying molecular mechanisms in sodium ion permeation and drug interactions. The results presented here could be useful for designing safer drugs that do not block the hNav1.5 channel.

Photoelectron Spectra and Electronic Structures of the Radiosensitizer Nimorazole and Related Compounds

Soft X-ray photoelectron spectroscopy has been used to investigate the radiosensitizer nimorazole and related model compounds. We report the valence and C, N, and O 1s photoemission spectra and K-edge NEXAFS spectra of gas-phase nimorazole, 1-methyl-5-nitroimidazole, and 4(5)-nitroimidazole in combination with theoretical calculations. The valence band and core level spectra are in agreement with theory. We determine the equilibrium populations of the two tautomers in 4(5)-nitroimidazole and find a ratio of 1:0.7 at 390 K. The NEXAFS spectra of the studied nitroimidazoles show excellent agreement with spectra of compounds available in the literature that exhibit a similar chemical environment. By comparing 1-methyl-5-nitroimidazole (single tautomer) with 4(5)-nitroimidazole, we are able to disentangle the photoemission and photoabsorption spectra and identify features due to each single tautomer.

Structures of cycloserine and 2-oxazolidinone probed by X-ray photoelectron spectroscopy: theory and experiment.

The electronic structures and properties of 2-oxazolidinone and the related compound cycloserine (CS) have been investigated using theoretical calculations and core and valence photoelectron spectroscopy. Isomerization of the central oxazolidine heterocycle and the addition of an amino group yield cycloserine. Theory correctly predicts the C, N, and O 1s core spectra, and additionally, we report theoretical natural bond orbital (NBO) charges. The valence ionization energies are also in agreement with theory and previous measurements. Although the lowest binding energy part of the spectra of the two compounds shows superficial similarities, further analysis of the charge densities of the frontier orbitals indicates substantial reorganization of the wave functions as a result of isomerization. The highest occupied molecular orbital (HOMO) of CS shows leading carbonyl π character with contributions from other heavy (non-H) atoms in the molecule, while the HOMO of 2-oxazolidinone (OX2) has leading nitrogen, carbon, and oxygen pπ characters. The present study further theoretically predicts bond resonance effects of the compounds, evidence for which is provided by our experimental measurements and published crystallographic data.

Sitting above the maze: recent model discoveries in molecular science

The 2013 Nobel Prize in Chemistry has convinced the world that how important the role that computational sciences play in chemical and materials sciences. In this review, computational methods and rational molecule design, including quantum mechanics and molecular mechanics methods, have been applied to study electronic structures and the interactions in a number of important applications at molecular level. The applications which include bioactive compounds, drug candidates and photoactive molecules at Swinburne University in the past several years are discussed. The research is in close collaboration with world class experimental groups from spectroscopy, organic and medicinal synthesis laboratories and most recently to γ-ray spectroscopy as well as other theory groups in the world. Ionisation spectra of biomolecules and bioactive compounds including amino acids, DNA bases, cyclic dipeptides, drug candidates, complexes and photoactive molecules are discussed. Most recent projects such as infrared spectral studies of ferrocene, rational design of organic dyes in solar cell applications, and recent development in γ-ray spectra of positron annihilation in molecules are highlighted.

Effects of protein-protein interactions and ligand binding on the ion permeation in KCNQ1 potassium channel

The voltage-gated KCNQ1 potassium ion channel interacts with the type I transmembrane protein minK (KCNE1) to generate the slow delayed rectifier (IKs) current in the heart. Mutations in these transmembrane proteins have been linked with several heart-related issues, including long QT syndromes (LQTS), congenital atrial fibrillation, and short QT syndrome. Off-target interactions of several drugs with that of KCNQ1/KCNE1 ion channel complex have been known to cause fatal cardiac irregularities. Thus, KCNQ1/KCNE1 remains an important avenue for drug-design and discovery research. In this work, we present the structural and mechanistic details of potassium ion permeation through an open KCNQ1 structural model using the combined molecular dynamics and steered molecular dynamics simulations. We discuss the processes and key residues involved in the permeation of a potassium ion through the KCNQ1 ion channel, and how the ion permeation is affected by (i) the KCNQ1-KCNE1 interactions and (ii) the binding of chromanol 293B ligand and its derivatives into the complex. The results reveal that interactions between KCNQ1 with KCNE1 causes a pore constriction in the former, which in-turn forms small energetic barriers in the ion-permeation pathway. These findings correlate with the previous experimental reports that interactions of KCNE1 dramatically slows the activation of KCNQ1. Upon ligand-binding onto the complex, the energy-barriers along ion permeation path are more pronounced, as expected, therefore, requiring higher force in our steered-MD simulations. Nevertheless, pulling the ion when a weak blocker is bound to the channel does not necessitate high force in SMD. This indicates that our SMD simulations have been able to discern between strong and week blockers and reveal their influence on potassium ion permeation. The findings presented here will have some implications in understanding the potential off-target interactions of the drugs with the KCNQ1/KCNE1 channel that lead to cardiotoxic effects.