Quantum chemical modeling of the GTP hydrolysis by the RAS–GAP protein complex

We present results of the modeling for the hydrolysis reaction of guanosine triphosphate (GTP) in the RAS–GAP protein complex using essentially ab initio quantum chemistry methods. One of the approaches considers a supermolecular cluster composed of 150 atoms at a consistent quantum level. Another is a hybrid QM/MM method based on the effective fragment potential technique, which describes interactions between quantum and molecular mechanical subsystems at the ab initio level of the theory. Our results show that the GTP hydrolysis in the RAS–GAP protein complex can be modeled by a substrate-assisted catalytic mechanism. We can locate a configuration on the top of the barrier corresponding to the transition state of the hydrolysis reaction such that the straightforward descents from this point lead either to reactants GTP+H₂O or to products guanosine diphosphate (GDP)+H₂PO₄^?. However, in all calculations such a single-step process is characterized by an activation barrier that is too high. Another possibility is a two-step reaction consistent with formation of an intermediate. Here the Pγ-O(Pβ) bond is already broken, but the lytic water molecule is still in the pre-reactive state. We present arguments favoring the assumption that the first step of the GTP hydrolysis reaction in the RAS–GAP protein complex may be assigned to the breaking of the Pγ-O(Pβ) bond prior to the creation of the inorganic phosphate.     

Activated bleomycin. A transient complex of drug, iron, and oxygen that degrades DNA

"Activated bleomycin" is an oxygenated iron drug complex which embodies the drug's DNA-cleaving activity. This activity is exercised on DNA, if present, but if DNA is absent, the drug itself is inactivated. Hyperfine interactions in the EPR spectra of activated bleomycin prepared with 57Fe(II) and 17O2 demonstrate the presence of iron as Fe(III) and of bound oxygen originating in dioxygen. Bleomycin can also be activated with Fe(III) and either H2O2 or ethyl hydroperoxide. These latter reactions do not produce a ferrous intermediate nor do they require O2. But O2 is required for the reaction of activated bleomycin with DNA to yield the malondialdehyde-like chromogens used to monitor DNA degradation. The attack on DNA is quantitatively concurrent with the decay of activated bleomycin, however generated.     

Cooperativity effect involving drug–DNA/RNA intermolecular interaction: A B3LYP-D3 and MP2 theoretical investigation on ketoprofen⋯cytosine⋯H2O system

In order to examine the origin of the drug action and design new DNA/RNA-targeted drugs, the cooperativity effect involving drug–DNA/RNA intermolecular interaction in ketoprofen?cytosine?H₂O ternary system were investigated by the B3LYP, B3LYP-D3, and MP2 methods with the 6-311++G(2d,p) basis set. The thermodynamic cooperativity was also evaluated at 310.15 K. The N–H?O, O–H?O, O–H?N, C–H?N, and C–H?O H bonds coexist in ternary complexes. The intermolecular interactions obtained by B3LYP-D3 are close to those calculated by MP2. The steric effects and van der Waals interactions have little influence on the cooperativity effects. The anti-cooperativity effect in ket?cyt?H₂O is far more notable than the cooperativity effect, and the stability of the cyclic structure with anti-cooperativity effect is higher than that of the linear structure with cooperativity effect, as is confirmed by the AIM (atoms in molecules) and RDG (reduced density gradient) analysis. Thus, it can be inferred that, in the presence of H₂O, the anti-cooperativity effect plays a dominant role in the drug–DNA/RNA interaction, and the nature of the hydration in the binding of drugs to DNA/RNA bases is the H-bonding anti-cooperativity effect. Furthermore, the drug always links simultaneously with DNA/RNA base and H₂O, and only in this way can the biological activity of drugs play a role. In most cases, the enthalpy change is the major factor driving the cooperativity, as is different from most of biomacromolecule complexes.     

Drug screening of α-amylase inhibitors as candidates for treating diabetes

In the present study, the identification of potential α-amylase inhibitors is explored as a potential strategy for treating type-2 diabetes mellitus. A computationally driven approach using molecular docking was employed to search for new α-amylase inhibitors. The interactions of potential drugs with the enzyme's active site were investigated and compared with the contacts established by acarbose (a reference drug for α-amylase inhibition) in the crystallographic structure 1B2Y. For this active site characterization, both molecular docking and molecular dynamics simulations were performed, and the residues involved in the α-amylase–acarbose complex were considered to analyse the potential drug's interaction with the enzyme. Two potential α-amylase inhibitors (AN-153I105594 and AN-153I104845) have been selected following this computational strategy. Both compounds established a large number of interactions with key binding site α-amylase amino acids and obtained a comparable docking score concerning the reference drug (acarbose). Aiming to further analyse candidates' properties, their ADME (absorption, distribution, metabolism, excretion) parameters, druglikeness, organ toxicity, toxicological endpoints and median lethal dose (LD₅₀) were estimated. Overall estimations are promising for both candidates, and in silico toxicity predictions suggest that a low toxicity should be expected.     

Biophysical and molecular docking insight into interaction mechanism and thermal stability of human serum albumin isoforms with a semi-synthetic water-soluble camptothecin analog irinotecan hydrochloride

In the present work, we have examined the binding parameters, thermodynamics, and stability of human serum albumin (HSA) isoforms at pH 7.4 and 9.0, using spectroscopic, calorimetric, and molecular docking methods in the presence of water-soluble camptothecin analog irinotecan hydrochloride (CPT-11). We observed that CPT-11 binds to HSA through a static quenching procedure of ground-state complex formation with N-isoform and B-isoform. Hydrogen bond and hydrophobic interactions are the major governing forces that participating in the formation of protein–drug complex. To determine the binding site of CPT-11 within HSA molecules, we also have performed molecular docking experiments. We explored the CPT-11-mediated stability and modulation of HSA by performing dynamic light scattering (DLS) and differential scanning calorimetry (DSC) experiments. DLS and DSC techniques are used to determine the size and the melting point (T_m) of HSA, which was decreased in the presence of CPT-11. Therefore, CPT-11 plays an important role in HSA stability and protein–ligand interactions. The present study provides valuable information in the field of pharmacokinetics, pharmaco-dynamics, and drug discovery.     

Analysis of proton wires in the enzyme active site suggests a mechanism of c‐di‐GMP hydrolysis by the EAL domain phosphodiesterases

We report for the first time a hydrolysis mechanism of the cyclic dimeric guanosine monophosphate (c‐di‐GMP) by the EAL domain phosphodiesterases as revealed by molecular simulations. A model system for the enzyme‐substrate complex was prepared on the base of the crystal structure of the EAL domain from the BlrP1 protein complexed with c‐di‐GMP. The nucleophilic hydroxide generated from the bridging water molecule appeared in a favorable position for attack on the phosphorus atom of c‐di‐GMP. The most difficult task was to find a pathway for a proton transfer to the O3' atom of c‐di‐GMP to promote the O3'? P bond cleavage. We show that the hydrogen bond network extended over the chain of water molecules in the enzyme active site and the Glu359 and Asp303 side chains provides the relevant proton wires. The suggested mechanism is consistent with the structural, mutagenesis, and kinetic experimental studies on the EAL domain phosphodiesterases. Proteins 2016; 84:1670–1680. © 2016 Wiley Periodicals, Inc.     

Molecular modeling studies of the artemisinin (qinghaosu)-hemin interaction: Docking between the antimalarial agent and its putative receptor

Artemisinin (qinghaosu, QHS) is a promising new antimalarial agent that is effective against drug-resistant strains of malaria. The antimalarial activity of this drug appears to be mediated by an interaction of the drug's endoperoxide bridge with intraparasitic hemin. We have carried out a computer-assisted docking of QHS with hemin from various starting configurations and found that, in the most stable docked configuration, the endoperoxide bridge is in close proximity to the hemin iron. In contrast, an inactive analog, deoxyartemisinin (DQHS), docks in a different manner. Further computer analysis of the drug-hemin interaction might aid in the design of new QHS congeners.     

Structural impact of the leukemia drug 1-beta-D-arabinofuranosylcytosine (Ara-C) on the covalent human topoisomerase I-DNA complex

1-beta-d-Arabinofuranosylcytosine (Ara-C) is a potent antineoplastic drug used in the treatment of acute leukemia. Previous biochemical studies indicated the incorporation of Ara-C into DNA reduced the catalytic activity of human topoisomerase I by decreasing the rate of single DNA strand religation by the enzyme by 2-3-fold. We present the 3.1 A crystal structure of human topoisomerase I in covalent complex with an oligonucleotide containing Ara-C at the +1 position of the non-scissile DNA strand. The structure reveals that a hydrogen bond formed between the 2'-hydroxyl of Ara-C and the O4' of the adjacent -1 base 5' to the damage site stabilizes a C3'-endo pucker in the Ara-C arabinose ring. The structural distortions at the site of damage are translated across the DNA double helix to the active site of human topoisomerase I. The free sulfhydryl at the 5'-end of the nicked DNA strand in this trapped covalent complex is shifted out of alignment with the 3'-phosphotyrosine linkage at the catalytic tyrosine 723 residue, producing a geometry not optimal for religation. The subtle structural changes caused by the presence of Ara-C in the DNA duplex may contribute to the cytotoxicity of this leukemia drug by prolonging the lifetime of the covalent human topoisomerase I-DNA complex.     

Structural Basis for Human DNA Polymerase Kappa to Bypass Cisplatin Intrastrand Cross-Link (Pt-GG) Lesion as an Efficient and Accurate Extender

Cisplatin (cis-diamminedichloroplatinum) is a common chemotherapeutic drug that reacts with the N7 atoms of adjacent guanines in DNA to form the Pt-1,2-d(GpG) intrastrand cross-link (Pt-GG), a major product to block DNA replication. Translesion DNA synthesis has been implicated in chemoresistance during cisplatin treatment of cancer due to Pt-GG lesion bypass. Gene knockdown studies in human cells have indicated a role for polκ during translesion synthesis of the Pt-GG lesion. However, the bypass activity of polκ with cisplatin lesions has not been well characterized. In this study, we investigated polκ's ability to bypass Pt-GG lesion in vitro and determined two crystal structures of polκ in complex with Pt-GG DNA. The ternary complex structures represent two consecutive stages of lesion bypass: nucleotide insertion opposite the 5′G (Pt-GG2) and primer extension immediately after the lesion (Pt-GG3). Our biochemical data showed that polκ is very efficient and accurate in extending DNA primers after the first G of the Pt-GG lesion. The structures demonstrate that the efficiency and accuracy is achieved by stably accommodating the bases with the cisplatin adduct in the active site for proper Watson–Crick base pairing with the incoming nucleotide in both the second insertion and post-insertion complexes. Our studies suggest that polκ works as an extender for efficient replication of the Pt-GG lesion in cells. This work holds promise for considering polκ, along with polη, as potential targets for drug design, which together could improve the efficacy of cisplatin treatment for cancer therapy.     

Relationship of inhibition of prostaglandin biosynthesis by analgesics to asthma attacks in aspirin-sensitive patients.

Eleven patients with asthma and aspirin hypersensitivity have been challenged with eight non-steroidal anti-inflammatory drugs. Each drug was given by mouth in at least three different doses and the patients'' symptoms and peak expiratory flow (PEF) rates were observed over a three-hour period. Indomethacin 5 mg caused bronchoconstriction in all patients. Therapeutic doses of mefenamic or flufenamic acid caused bronchoconstriction in most patients. Phenylbutazone 200-400 mg induced a moderate fall in PEF. There were no reactions to therapeutic doses of salicylamide, paracetamol, benzydamine, and chloroquine. Microsomal prostaglandin synthetase, activity was inhibited by aspirin, indomethacin, mefenamic acid, flufenamic acid, and phenylbutazone. The other four drugs had no inhibitory effect. We suggest that precipitation of attacks in asthmatic patients hypersensitive to certain anti-inflammatory drugs is related to drug''s ability to inhibit prostaglandin biosynthesis.     

Nonorthogonal tight-binding model with H–C–N–O parameterisation

A parametric nonorthogonal tight-binding model (NTBM1) with the set of parameters for H–C–N–O systems is presented. This model compares well with widely used semi-empirical AM1 and PM3/PM7 models but contains less fitting parameters per atom. All NTBM1 parameters are derived based on a criterion of the best agreement between the calculated and experimental values of bond lengths, valence angles and binding energies for various H–C–N–O molecules. Results for more than 200 chemical compounds are reported. Parameters are currently available for hydrogen, carbon, nitrogen, oxygen atoms and corresponding interatomic interactions. The model has a good transferability and can be used for both relaxation of large molecular systems (e.g., high-molecular compounds or covalent cluster complexes) and long-timescale molecular dynamics simulation (e.g., modelling of thermal decomposition processes). The program package based on this model is available for download at no cost from http://ntbm.info.     

Involvement of a low molecular weight component(s) in the mechanism of action of the glucocorticoid receptor.

[³H]Dexamethasone-receptor complexes from rat liver cytosol preincubated at 0° bind poorly to DNA-cellulose. However, if the steroid-receptor complex is subjected to gel filtration at 0–4° separating it from the low molecular weight components of cytosol, the steroid-receptor complex becomes “activated” enabling its binding to DNA-cellulose. This activation can be prevented if the gel filtration column is first equilibrated with the low molecular weight components of cytosol. In addition, if adrenalectomized rat liver cytosol, in the absence of exogeneous steroid, is subjected to gel filtration the macromolecular fractions separated from the “small molecules” of that cytosol have much reduced binding activity towards [³H]dexamethasone. These results suggest that rat liver cytosol contains a low molecular weight component(s) which maintains the glucocorticoid receptor in a conformational state that allows the binding of dexamethasone. Furthermore, this component must be removed from the steroid-receptor complex before binding to DNA can occur.     

Improving anterior deltoid activity in a musculoskeletal shoulder model – an analysis of the torque-feasible space at the sternoclavicular joint

Modelling the shoulder's musculature is challenging given its mechanical and geometric complexity. The use of the ideal fibre model to represent a muscle's line of action cannot always faithfully represent the mechanical effect of each muscle, leading to considerable differences between model-estimated and in vivo measured muscle activity. While the musculo–tendon force coordination problem has been extensively analysed in terms of the cost function, only few works have investigated the existence and sensitivity of solutions to fibre topology. The goal of this paper is to present an analysis of the solution set using the concepts of torque-feasible space (TFS) and wrench-feasible space (WFS) from cable-driven robotics. A shoulder model is presented and a simple musculo–tendon force coordination problem is defined. The ideal fibre model for representing muscles is reviewed and the TFS and WFS are defined, leading to the necessary and sufficient conditions for the existence of a solution. The shoulder model's TFS is analysed to explain the lack of anterior deltoid (DLTa) activity. Based on the analysis, a modification of the model's muscle fibre geometry is proposed. The performance with and without the modification is assessed by solving the musculo–tendon force coordination problem for quasi-static abduction in the scapular plane. After the proposed modification, the DLTa reaches 20% of activation.     

Interaction between Rad9–Hus1–Rad1 and TopBP1 activates ATR–ATRIP and promotes TopBP1 recruitment to sites of UV-damage

The checkpoint clamp Rad9–Hus1–Rad1 (9–1–1) interacts with TopBP1 via two casein kinase 2 (CK2)-phosphorylation sites, Ser-341 and Ser-387 in Rad9. While this interaction is known to be important for the activation of ATR-Chk1 pathway, how the interaction contributes to their accumulation at sites of DNA damage remains controversial. Here, we have studied the contribution of the 9–1–1/TopBP1 interaction to the assembly and activation of checkpoint proteins at damaged DNA. UV-irradiation enhanced association of Rad9 with chromatin and its localization to sites of DNA damage without a direct interaction with TopBP1. TopBP1, as well as RPA and Rad17 facilitated Rad9 recruitment to DNA damage sites. Similar to Rad9, TopBP1 also localized to sites of UV-induced DNA damage. The DNA damage-induced TopBP1 redistribution was delayed in cells expressing a TopBP1 binding-deficient Rad9 mutant. Pharmacological inhibition of ATR recapitulated the delayed accumulation of TopBP1 in the cells, suggesting that ATR activation will induce more efficient accumulation of TopBP1. Taken together, TopBP1 and Rad9 can be independently recruited to damaged DNA. Once recruited, a direct interaction of 9–1–1/TopBP1 occurs and induces ATR activation leading to further TopBP1 accumulation and amplification of the checkpoint signal. Thus, we propose a new positive feedback mechanism that is necessary for successful formation of the damage-sensing complex and DNA damage checkpoint signaling in human cells.     

Systematical investigation of binding interaction between novel ruthenium(II) arene complex with curcumin analogs and ctDNA

In this study, the interaction between a novel ruthenium(II) arene complex with curcumin analogs and calf thymus DNA (ctDNA) was investigated systematically by viscosity measurement, the DNA melting approach, multispectroscopic techniques and electrochemical methods. The absorption spectra of the ctDNA–drug complex showed a slight red shift and a weak hypochromic effect. The relative viscosity and melting temperature of ctDNA increased on addition of the drug. The evidence obtained from fluorescence competitive experiments indicated that the binding mode of the drug with ctDNA was intercalative. Using acridine orange (AO) as a fluorescence probe, the drug statically quenched the fluorescence of the ctDNA–AO complex, and hydrogen bonding and van der Waals interactions played vital roles in the binding interaction between the drug and ctDNA. The influences of ionic strength, chemical denaturants and pH on the binding interaction were also investigated. Circular dichroism and Fourier transform infrared spectra suggested that this drug might bond with the G–C base pairs of ctDNA and the right‐handed B‐form helicity of ctDNA remained after drug binding. The intercalative binding between the drug and ctDNA was further investigated using electrochemical techniques. All these results suggested that the biological activity of ctDNA was affected by ruthenium(II) arene complex with curcumin analogs. Copyright © 2016 John Wiley & Sons, Ltd.     

Base specific interaction of reductively activated nitroimidazoles with DNA

To exert biological activity, nitroimidazole drugs require reductive activation in vivo. Nucleic acids are susceptible to the activated drug in vitro and are presumably the major target in vivo. We carried out electrolytical reduction of several 5-nitroimidazoles at a controlled potential either in the presence or prior to the addition of DNA. Using a nucleotide sequence specific test to analyse cleavage products, specific interaction of the reduced nitroimidazole intermediate(s) towards the guanine residues is prominent. Since the strand scission depends on subsequent piperidine treatment, it can be concluded that the primary interaction between the activated drug and guanine is a covalent modification weakening the glycosidic bond.     

Quantum chemical modeling of the GTP hydrolysis by the RAS-GAP protein complex

We present results of the modeling for the hydrolysis reaction of guanosine triphosphate (GTP) in the RAS-GAP protein complex using essentially ab initio quantum chemistry methods. One of the approaches considers a supermolecular cluster composed of 150 atoms at a consistent quantum level. Another is a hybrid QM/MM method based on the effective fragment potential technique, which describes interactions between quantum and molecular mechanical subsystems at the ab initio level of the theory. Our results show that the GTP hydrolysis in the RAS-GAP protein complex can be modeled by a substrate-assisted catalytic mechanism. We can locate a configuration on the top of the barrier corresponding to the transition state of the hydrolysis reaction such that the straightforward descents from this point lead either to reactants GTP+H(2)O or to products guanosine diphosphate (GDP)+H(2)PO(4)(-). However, in all calculations such a single-step process is characterized by an activation barrier that is too high. Another possibility is a two-step reaction consistent with formation of an intermediate. Here the Pgamma-O(Pbeta) bond is already broken, but the lytic water molecule is still in the pre-reactive state. We present arguments favoring the assumption that the first step of the GTP hydrolysis reaction in the RAS-GAP protein complex may be assigned to the breaking of the Pgamma-O(Pbeta) bond prior to the creation of the inorganic phosphate.     

Spectroscopic study on the interaction of bovine serum albumin with zinc(II) phthalocyanine

The interaction between the photosensitive antitumour drug, 2(3),9(10),16(17),23(24)‐tetra‐(((2‐aminoethylamino)methyl)phenoxy)phthalocyaninato‐zinc(II) (ZnPc) and bovine serum albumin (BSA) has been investigated using various spectroscopic methods. This work may provide some useful information for understanding the interaction mechanism of anticancer drug–albumin binding and gain insight into the biological activity and metabolism of the drug in blood. Based on analysis of the fluorescence spectra, ZnPc could quench the intrinsic fluorescence of BSA and the quenching mechanism was static by forming a ground state complex. Meanwhile, the Stern–Volmer quenching constant (K_SV), binding constant (K_b), number of binding sites (n) and thermodynamic parameters were obtained. Results showed that the interaction of ZnPc with BSA occurred spontaneously via hydrogen bond and van der Waal's force. According to Foster's non‐radioactive energy transfer theory, the energy transfer from BSA to ZnPc occurred with high possibility. Synchronous fluorescence and circular dichroism (CD) spectra also demonstrated that ZnPc induced the secondary structure of and conformation changes in BSA, especially α helix. Copyright © 2015 John Wiley & Sons, Ltd.     

O-O bond splitting mechanism in cytochrome oxidase

Hybrid density functional theory (DFT) calculations have been used to investigate different mechanisms for O–O bond splitting in cytochrome oxidase. It is shown that the requirement for a low activation barrier for the O–O bond splitting is that two protons, apart from the tyrosine hydroxyl proton, are available at the binuclear center. A mechanism is suggested for the transformation from a species with a molecularly coordinated O₂, to an O–O cleaved species with an oxo-ferryl group. The mechanism has a calculated activation barrier in reasonable agreement with experimental estimates, and the overall reaction is close to thermoneutral, in line with the requirement that the energy wasted as heat should be minimized. The rate limiting step in the mechanism occurs at the initial Fe–O₂ intermediate, consistent with experimental observations that the decay of the oxy intermediate parallels the increase of the oxo product. The formation of a radical at the cross-linked tyrosine–histidine structure is a possible source for one of the electrons required in the bond cleavage process. Possible sources for the two protons are discussed, including a suggested key role for the hydroxyl group on the farnesyl side chain of heme a₃.