Theoretical density functional theory studies on interactions of small biologically active molecules with isolated heme group

We present ab-initio density functional theory studies on the interactions of small biologically active molecules, namely NO, CO, O(2), H(2)O, and NO(2) (-) with the full-size heme group. Our results show that the small molecule-iron bond is the strongest in carbonyl and the weakest in nitrite system. Trans influence induced by NO binding to the five-coordinate heme complex is shown. Nitric oxide in the resulting complex might be described as NO(-). The differences among the small ligands of XO type (CO, NO, O(2)), and their distant chemical behavior from H(2)O and NO(2) (-) ligands in binding to the Fe(II) ion, are shown. Moreover, the role of the heme ring as a reservoir of electrons in the studied complexes is invoked. The analysis of the parameters defining the iron-histidine bond indicates that this bond is longer and weaker in nitrosyl and carbonyl complexes than in the other systems. Our findings support the proposed mechanism of soluble guanylate cyclase (sGC) activation and suggest that the first step of sGC activation by CO may be the same as during the activation by NO. Obtained results are then compared with the data concerning smaller model of the heme, the porphyrin complexes, available in the literature.     

Theoretical treatment of multiple site reactivity in large molecules

The method of electrostatic potentials is generalized via a double perturbation theory to describe the simultaneous attack of a large molecule by two or more charged species. The interaction term is found to reproduce accurately the full “super-molecule” calculation. The scheme is illustrated by treating the diprotonation of adenine and N7-methyladenine. For adenine, the most stable diprotonated species is found to be the N1, N7 combination, while the N3, N9 couple is predicted for the N7 substituted derivative if N3 is the first protonation site.     

Physicochemical simulation and experimental study of interactions in the Si-Al-Ca-Mg-Fe-Na-Cr-Cl-H-O system

The interactions in the SiO₂-Al₂O₃-CaO-MgO-Fe₂O₃-FeO-Na₂O-Cr₂O₃-CaCl₂-HCl-H₂O system were studied by physicochemical simulation methods. Conditions for the most complete extraction of chromium into hydrochloric acid solutions were estimated. The equilibrium compositions of reaction products were calculated. The effect of calcium chloride on the solubility of chromium-containing minerals in hydrochloric acid was determined.     

Physical nature of interactions within the active site of cytosine-5-methyltransferase

The physical nature of interactions within the active site of cytosine-5-methyltransferase (CMT) was studied using a variation-perturbation energy decomposition scheme defining a sequence of approximate intermolecular interaction energy models. These models have been used to analyze the catalytic activity of residues constituting cytosine-5-methyltransferase active site as well their role in the binding group of de novo designed inhibitors. Our results indicate that Glu119, Arg163, and Arg165 appear to play the dominant role in stabilizing the protonated transition state structure and their influence can be qualitatively approximated by electrostatic interactions alone. The stabilization of neutral structures of the alternative reaction pathway is small, which might suggest the protonated pathway as preferred by the enzyme. Exchange and delocalization terms are negligible in most cases, or they cancel each other to some extent. Interactions of inhibitors with the CMT active site are dominated by electrostatic multipole contributions in analogy with previously studied transition state analogue inhibitors of leucyl aminopeptidase.     

Theoretical and experimental investigations of the electronic Rydberg states of diazomethane: assignments and state interactions

The electronic states of diazomethane in the region 3.00-8.00 eV have been characterized by ab initio calculations, and electronic transitions in the region 6.32-7.30 eV have been examined experimentally using a combination of 2 + 1 REMPI spectroscopy and photoelectron imaging in a molecular beam. In the examined region, three Rydberg states of 3p character contribute to the transitions, 2(1)A2(3p(y) <-- pi), 2(1)B1(3p(z) <-- pi), and 3(1)A1(3p(x) <-- pi). The former two states are of mostly pure Rydberg character and exhibit a resolved K structure, whereas the 3(1)A1(3p(x) <-- pi) state is mixed with the valence 2(1)A1(pi* <-- pi) state, which is unbound and is strongly predissociative. Analyses of photoelectron kinetic energy distributions indicate that the ground vibrational level of the 2(1)B1(3p(z)) state is mixed with the 2(1)A2(3p(y)) nu(9) level, which is of B1 vibronic symmetry. The other 2(1)A2(3p(y)) vibronic states exhibit pure Rydberg character, generating ions in single vibrational levels. The photoelectron spectra of the 3(1)A1(3p(x) <-- pi) state, on the other hand, give rise to many states of the ion as a result of strong mixing with the valence state, as evidenced also in the ab initio calculations. The equilibrium geometries of the electronic states of neutral diazomethane were calculated by CCSD(T), using the cc-pVTZ basis, and by B3LYP, using the 6-311G(2df,p) basis. Geometry and frequencies of the ground state of the cation were calculated by CCSD(T)/cc-pVTZ, using the unrestricted (UHF) reference. Vertical excitation energies were calculated using EOM-CCSD/6-311(3+,+)G* at the B3LYP optimized geometry. The theoretical results show that the 2(1)A2(3p(y) <-- pi) and 2(1)B1(3p(z) <-- pi) states have geometries similar to the ion, which has C(2v) symmetry, with slight differences due to the interactions of the electron in the 3p orbital with the nuclei charge distributions. The geometry of the 3(1)A1(3p(x) <-- pi) state is quite different and has Cs symmetry. The experimental and theoretical results agree very well, both in regard to excitation energies and to vibrational modes of the ion.     

Theoretical studies on substrate binding to the active site of carbonic anhydrase

The role of some amino acids and metal ions in the catalytic activity of carbonic anhydrase (carbonate dehydratase EC 4.2.1.1) has been investigated. The additional stabilization of the transition state complex was used as the qualitative measure of the effect of molecular surrounding on CO₂ hydration reaction calculated within the approximate CNDO /2 approach. The effect of the molecular environment has been simulated by inclusion into the SCF LCAO MO Hamiltonian, a term representing the interaction with a set of point charges and atomic dipoles centered on experimentally determined position of Zn²⁺ ion and all atoms of histidine 94, 96, 119 and threonine 199, 200 (or histidine 200 in the case of carbonic anhydrase B). The possible molecular mechanism of CO₂ hydration inside the active site has been also discussed.     

Theoretical and experimental studies of the structure and intermolecular interactions in 2-biphenylmethanole

The crystalline structure of 2-biphenylmethanol has been investigated by the radiographic method at room temperature, and its IR transmission spectra were measured in the 400–4000 cm^?1 range. Structural-dynamic models of free molecules of 2-biphenylmethanol, biphenyl and methanol, hydrogen-bonded complexes of 2-biphenylmethanol with methanol and a tetramer of methanol molecules have been built by the density functional method (B3LYP/6-31G*). The structure, the energy, mechanical and electro-optical parameters, the frequency and intensity of normal vibrations in IR spectra of each of the mentioned molecular systems were calculated. The crystalline structure of 2-biphenylmethanol, its features determined by the hydrogen bond formation in the crystal has been found, its structure-forming role and energy were assessed, and IR spectra interpreted on the basis of single crystal X-ray diffraction data, the analysis of measured IR spectra, and results of molecular modeling.     

细胞死亡蛋白的三维结构及活性位点的理论研究

利用同源模建和分子动力学模拟方法,构建了细胞死亡蛋白（CED-10）的三维结构,并预测了其可能的活性位点,为细胞死亡蛋白家族催化机理的深入研究提供了、重要的参考信息．     

Theoretical study of the unusual protonation properties of the active site cysteines in thioredoxin

The thiol/disulfide oxidoreductases of the thioredoxin family have, in the active site, two cysteines that can be in a reduced or an oxidized form. One of the cysteines in the reduced state is deprotonated, and it is called nucleophilic cysteine. The pK(a) of this cysteine is different from that of a normal cysteine and varies widely among the different enzymes of this family. However, the factors responsible for the different degrees of stabilization of nucleophilic cysteine thiolate are not fully understood. Here, we have studied the well-known hypothesis of proton sharing between the active site thiols by performing a linear transit scan for the transfer of the proton between the active site cysteines. We used a two-layered (DFT/MM) ONIOM formalism, with the active site region treated at the B3LYP/6-31+G(d) level and the remains of the protein treated with the Amber Parm94 force field. The solvation free energy was accounted for with a continuum solvent model, by solving the Poisson-Boltzmann equation using the program Delphi. We have obtained excellent agreement with the experimental data available in the literature. Besides refuting the proton sharing hypothesis, our results include a value of 14.0 for the pK(a) of the buried cysteine, a quantity that has not been possible to obtain experimentally but which has been proven to be higher than 11. Additionally, this study also provides detailed information on the very interesting and so far unknown fact that the contribution of the enzymatic structure (8.3 kcal/mol) prevails in relation to that of the solvent (0.60 kcal/mol) concerning the differential stabilization of the active site thiolates.     

Theoretical and experimental investigation of stereoelectronic interactions in H-complexes of sulfenamide derivatives

The photoelectron and IR spectra of a number of sulfenamide derivatives and their H-complexes have been investigated, A correlation between an increase in the vertical ionization potential of the lone electron pair of the nitrogen atom and a decrease in the frequency shift of the stretching OH-vibrations in the H-complexes of compounds R₃N, R₂NCH₂OR, R₂NSR, and R₂NSOR was found. The electronic and geometric structures of the starting bases and their H-complexes were calculated by theab initio and MNDO methods. Anomeric interactions were found to decrease the energy of the n(N) orbital and to hinder the formation of H-complexes. The calculations of the sulfenamides and their H-complexes in unstable conformations, characterized by increased energies of H-complexation and proton affinity, were also carried out.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2894–2897, December, 1996.     

Automatic superposition of drug molecules based on their common receptor site

Summary We have prevously developed a new rational method for superposing molecules in terms of submolecular physical and chemical properties, but not in terms of atom positions or chemical structures as has been done in the conventional methods. The program was originally developed for interactive use on a three-dimensional graphic display, providing goodness-of-fit indices on molecular shape, hydrogen bonds, electrostatic interactions and others.Here, we report a new unbiased searching method for the best superposition of molecules, covering all the superposing modes and conformational freedom, as an additional function of the program. The function is based on a novel least-squares method which superposes the expected positions and orientations of hydrogen bonding partners in the receptor that are deduced from both molecules. The method not only gives reliability and reproducibility to the result of the superposition, but also allows us to save labor and time. It is demonstrated that this method is very efficient for finding the correct superposing mode in such systems where hydrogen bonds play important roles.     

Theoretical studies of [FeFe]-hydrogenase: infrared fingerprints of the dithiol-bridging ligand in the active site

An unresolved structural issue for [FeFe]-hydrogenases is the nature of the dithiol-bridging ligand in the diiron subcluster of the active site. The two most probable candidates are 1,3-dithiopropane (propane dithiol, PDT) and di-(thiomethyl)-amine (DTN). In the latter case, the dithiol-bridging ligand is assumed to play a major role in the reaction cycle. We report density-functional theory studies of the differing roles of these dithiol-bridging ligands in the infrared spectra of synthetic models and of computational representations of the diiron cluster of the active site. Our analysis shows distinct spectral features associated with the dithiol-bridging NH mode for compounds having a DTN bridge, which, however, would have been obscured by the H2O vibrations in existing measurements. However, if indeed nitrogen is present in the dithiol-bridging ligand, a combination of selective deuteration and chemical inactivation with CO would create a unique signature in an accessible region of the infrared spectrum, whose position and intensity are predicted.     

Characterization of an RNA active site: interactions between a Diels-Alderase ribozyme and its substrates and products

Ribozymes have recently been shown to catalyze the stereoselective formation of carbon-carbon bonds between small organic molecules. The interactions of these Diels-Alderase ribozymes with their substrates and products have now been elucidated by chemical substitution analysis by using 44 different, systematically varied analogues. RNA-diene interaction is governed by stacking interactions, while hydrogen bonding and metal ion coordination appear to be less important. The diene has to be an anthracene derivative, and substituents at defined positions are permitted, thereby shedding light on the geometry of the binding site. The dienophile must be a five-membered maleimidyl ring with an unsubstituted reactive double bond, and a hydrophobic side chain makes a major contribution to RNA binding. The ribozyme distinguishes between different enantiomers of chiral substrates and accelerates cycloadditions with both enantio- and diastereoselectivity. The stereochemistry of the reaction is controlled by RNA-diene interactions. The RNA interacts strongly and stereoselectively with the cycloaddition products, requiring several structural features to be present. Taken together, the results highlight the intricacy of ribozyme active sites which can control chemical reaction pathways based on minute differences in substrate stereochemistry and substitution pattern.     

Influence of protonation on substrate and inhibitor interactions at the active site of human monoamine oxidase-A

Although substrate conversion mediated by human monoaminooxidase (hMAO) has been associated with the deprotonated state of their amine moiety, data regarding the influence of protonation on substrate binding at the active site are scarce. Thus, in order to assess protonation influence, steered molecular dynamics (SMD) runs were carried out. These simulations revealed that the protonated form of the substrate serotonin (5-HT) exhibited stronger interactions at the protein surface compared to the neutral form. The latter displayed stronger interactions in the active site cavity. These observations support the possible role of the deprotonated form in substrate conversion. Multigrid docking studies carried out to rationalize the role of 5-HT protonation in other sites besides the active site indicated two energetically favored docking sites for the protonated form of 5-HT on the enzyme surface. These sites seem to be interconnected with the substrate/inhibitor cavity, as revealed by the tunnels observed by means of CAVER program. pK(a) calculations in the surface loci pointed to Glu32?, Asp32?, His???, and Asp132 as candidates for a possible in situ deprotonation step. Docking analysis of a group of inhibitors (structurally related to substrates) showed further interactions with the same two docking access sites. Interestingly, the protonated/deprotonated amine moiety of almost all compounds attained different docking poses in the active site, none of them oriented to the flavin moiety, thus producing a more variable and less productive orientations to act as substrates. Our results highlight the role of deprotonation in facilitating substrate conversion and also might reflect the necessity of inhibitor molecules to adopt specific orientations to achieve enzyme inhibition.     

The measurement of molecular diversity by receptor site interaction simulation

The assembly of large compound libraries for the purpose of screening against various receptor targets to identify chemical leads for drug discovery programs has created a need for methods to measure the molecular diversity of such libraries. The method described here, for which we propose the acronym RESIS (for Receptor Site Interaction Simulation), relates directly to this use. A database is built of three-dimensional representations of the compounds in the library and a set of three-point three-dimensional theoretical receptor sites is generated based on putative hydrophobic and polar interactions. A series of flexible, three-dimensional searches is then performed over the database, using each of the theoretical sites as the basis for one such search. The resulting pattern of hits across the grid of theoretical receptor sites provides a measure of the molecular diversity of the compound library. This can be conveniently displayed as a density map which provides a readily comprehensible visual impression of the library diversity characteristics. A library of 7500 drug compounds derived from the CIPSLINEPC databases was characterized with respect to molecular diversity using the RESIS method. Some specific uses for the information obtained from application of the method are discussed. A comparison was made of the results from the RESIS method with those from a recently published two-dimensional approach for assessing molecular diversity using sets of compounds from the Maybridge database (MAY).     

The tethered agonist approach to mapping ion channel proteins--toward a structural model for the agonist binding site of the nicotinic acetylcholine receptor

BACKGROUND: The integral membrane proteins of neurons and other excitable cells are generally resistant to high resolution structural tools. Structure-function studies, especially those enhanced by the nonsense suppression methodology for unnatural amino acid incorporation, constitute one of the most powerful probes of ion channels and related structures. The nonsense suppression methodology can also be used to incorporate functional side chains designed to deliver novel structural probes to membrane proteins. In this vein, we sought to generalize a potentially powerful tool - the tethered agonist approach - for mapping the agonist binding site of ligand-gated ion channels. RESULTS: Using the in vivo nonsense suppression method for unnatural amino acid incorporation, a series of tethered quaternary ammonium derivatives of tyrosine have been incorporated into the nicotinic acetylcholine receptor. At three sites a constitutively active receptor results, but the pattern of activation as a function of chain length is different. At position alpha149, there is a clear preference for a three-carbon tether, while at position alpha93 tethers of 2-5 carbons are comparably effective. At position gamma55/delta57 all tethers except the shortest one can activate the receptor. Based on these and other data, a model for the receptor binding site can be developed by analogy to the acetylcholine esterase crystal structure. CONCLUSION: Through the use of nonsense suppression techniques, the tethered agonist approach has been made into a general tool for probing receptor structures. When applied to the nicotinic receptor, the method places new restrictions on developing models for the agonist binding site.     

Molecular dynamics simulation of the human adenosine A3 receptor: agonist induced conformational changes of Trp243

The adenosine A(3) receptor together with rhodopsin belongs to Class A of the G-protein coupled receptors. As the crystal structure of bovine rhodopsin represents the dark (inactive) state of the receptor, the details of GPCR activation are still unknown. In this molecular dynamics study we investigate how the homology model of the human adenosine A(3) receptor responds to ligand exposure. To this end we placed the homology model in a POPC membrane model. After equilibrating for 13 ns an agonist (Cl-IB-MECA) and an inverse agonist (PSB-10) were placed inside the putative binding pocket. In the following 10 ns molecular dynamics simulation we observed a different behaviour of the side-chain torsions of Trp243(6.48), depending on the presence or absence of the agonist or inverse agonist. This conformational change of Trp243 correlates with the assumed influence of ligands on receptor activation. Other predicted conformational changes of the receptor could not be observed yet. So Trp243 may represent the first switch in receptor activation.