Expression of XBP1 in lymphocytes of the small intestine in rats under chronic social stress and modulation of intestinal microflora composition.

The present study was conducted to investigate of the influence of chronic social stress and modulation of the composition of intestinal microflora on the distribution of Xbp(1+)-lymphocytes in the gut-associated lymphoid tissue of ileum of the rats. Structure of population of Xbp(1+)-cells has been studied by the analysis of serial histological sections using the method of indirect immunofluorescence with monoclonal antibodies to Xbp1 of rat. Chronic social stress development is accompanied with the reduction of total number of Xbp(1+)-lymphocytes in lymphoid structures of ileum (31% -3 fold reduction, p < 0.05), mostly expressed in lymphoidfollicles, and changes the concentration of Xbp1 protein in immunopositive cells. Modulation of the composition of intestinal microflora by antibiotics and probiotics under chronic social stress results in the increase of total number of Xbp1+ lymphocytes in gut-associated lymphoid tissue, the degree of it depends on the kind ofstress. The discovered alterations of Xbp1 expression under stress may be one of the triggers for development of autoimmune and inflammatory bowel diseases. Thus, increased understanding of the molecular actions and transcriptional networks regulated by XBP1 in immune cells may aid in the development of potential therapeutics targeting immune disorders.

Mechanism of the chemical step for the guanosine triphosphate (GTP) hydrolysis catalyzed by elongation factor Tu.

Elongation factor Tu (EF-Tu), the protein responsible for delivering aminoacyl-tRNAs (aa-tRNAs) to ribosomal A site during translation, belongs to the group of guanosine-nucleotide (GTP/GDP) binding proteins. Its active 'on'-state corresponds to the GTP-bound form, while the inactive 'off'-state corresponds to the GDP-bound form. In this work we focus on the chemical step, GTP+H(2)O-->GDP+Pi, of the hydrolysis mechanism. We apply molecular modeling tools including molecular dynamics simulations and the combined quantum mechanical-molecular mechanical calculations for estimates of reaction energy profiles for two possible arrangements of switch II regions of EF-Tu. In the first case we presumably mimic binding of the ternary complex EF-Tu.GTP.aa-tRNA to the ribosome and allow the histidine (His85) side chain of the protein to approach the reaction active site. In the second case, corresponding to the GTP hydrolysis by EF-Tu alone, the side chain of His85 stays away from the active site, and the chemical reaction GTP+H(2)O-->GDP+Pi proceeds without participation of the histidine but through water molecules. In agreement with the experimental observations which distinguish rate constants for the fast chemical reaction in EF-Tu.GTP.aa-tRNA.ribosome and the slow spontaneous GTP hydrolysis in EF-Tu, we show that the activation energy barrier for the first scenario is considerably lower compared to that of the second case.

Conformation dependence of pKa's of the chromophores from the purple asFP595 and yellow zFP538 fluorescent proteins.

Two members of the green fluorescent protein family, the purple asFP595 and yellow zFP538 proteins, are perspective fluorescent markers for use in multicolor imaging and resonance energy-transfer applications. We report the results of quantum based calculations of the solution pKa values for selected protonation sites of the denatured asFP595 and zFP538 chromophores in the trans- and cis-conformations in order to add in the interpretation of photophysical properties of these proteins. The pKa values were determined from the theromodynamic cycle based on B3LYP/6-311++G(2df,2p) calculations of the gas phase free energies of the molecules and the B3LYP/6-311++G(d,p) calculations of solvation energies. The results show that the pKa's of the protonation sites of the chromophore from asFP595 noticeably depend on the isomer conformation (cis- or trans-), while those of zFP538 are much less sensitive to isomerization.

Flexible effective fragment QM/MM method: validation through the challenging tests.

A new version of the QM/MM method, which is based on the effective fragment potential (EFP) methodology [Gordon, M. et al., J Phys Chem A 2001, 105, 293] but allows flexible fragments, is verified through calculations of model molecular systems suggested by different authors as challenging tests for QM/MM approaches. For each example, the results of QM/MM calculations for a partitioned system are compared to the results of an all-electron ab initio quantum chemical study of the entire system. In each case we were able to achieve approximately similar or better accuracy of the QM/MM results compared to those described in original publications. In all calculations we kept the same set of parameters of our QM/MM scheme. A new test example is considered when calculating the potential of internal rotation in the histidine dipeptide around the C(alpha)bond;C(beta) side chain bond.

alpha- and 3(10)-helix interconversion: a quantum-chemical study on polyalanine systems in the gas phase and in aqueous solvent.

Helices are among the predominant secondary structures in globular proteins. About 90% of the residues in them are found to be in the alpha-helical conformation, and another 10% in the 3(10) conformation. There is a standing controversy between experimental and some theoretical results, and controversy among theoretical results concerning the predominance of each conformation, in particular, helices. We address this controversy by ab initio Hartree-Fock and density functional theory studies of helices with different lengths in a vacuum and in the aqueous phase. Our results show that (1) in a vacuum, all oligo(Ala) helices of 4-10 residues adopt the 3(10) - conformation; (2) in aqueous solution, the 6-10 residue peptides adopt the alpha-helical conformation; (3) there might be two intermediates between these helical conformers allowing for their interconversion. The relevance of these results to the structure and folding of proteins is discussed.

Experimental determination and calculations of redox potential descriptors of compounds directed against retroviral zinc fingers: Implications for rational drug design.

A diverse set of electrophilic compounds that react with cysteine thiolates in retroviral nucleocapsid (NC) proteins and abolish virus infectivity has been identified. Although different in chemical composition, these compounds are all oxidizing agents that lead to the ejection of Zn(II) ions bound to conserved structural motifs (zinc fingers) present in retroviral NC proteins. The reactivity of a congeneric series of aromatic disulfides toward the NC protein of the human immunodeficiency virus type 1 (HIV-1), NCp7, has been characterized by HPLC separation of starting reagents from reaction products. We calculated the absolute redox potentials of these compounds in the gas phase and in aqueous solvent, using a density functional theory method and a continuum solvation model. Pulsed polarography experiments were performed and showed a direct correlation between calculated and experimentally determined redox propensities. A dependence between protein reactivity and redox potential for a specific compound was shown: Reaction with NCp7 did not take place below a threshold value of redox potential. This relationship permits the distinction between active and nonactive compounds targeted against NCp7, and provides a theoretical basis for a scale of reactivity with retroviral zinc fingers. Our results indicate that electrophilic agents with adequate thiophilicity to react with retroviral NC fingers can now be designed using known or calculated electrochemical properties. This may assist in the design of antiretroviral compounds with greater specificity for NC protein. Such electrophilic agents can be used in retrovirus inactivation with the intent of preparing a whole-killed virus vaccine formulation that exhibits unaffected surface antigenic properties.

Theoretical calculations of glycine and alanine gas-phase acidities.

The gas-phase acidities of glycine and alanine were determined by using a variety of high level theoretical methods to establish which of these would give the best results with accessible computational efforts. MP2, MP4, QCISD, G2 ab initio procedures, hybrid Becke3-LYP (B3LYP) and gradient corrected Becke-Perdew (BP) and Perdew-Wang and Perdew (PWP) nonlocal density functionals were used for the calculations. A maximum deviation of approximately 13 and 18 kJ/mol from experimental data was observed for the computed delta Hacid and delta Gacid values, respectively. The best result was obtained at G2 level, but comparable reliability was reached when the considerably less time consuming B3LYP, BP, and PWP density functional approaches were employed.

Cyclopropane-derived peptidomimetics. Design, synthesis, evaluation, and structure of novel HIV-1 protease inhibitors.

Toward establishing the general efficacy of using trisubstituted cyclopropanes as peptide mimics to stabilize extended peptide structures, the cyclopropanes 20a-d were incorporated as replacements into 9-13, which are analogues of the known HIV-1 protease inhibitors 14 and 15. The syntheses of 20a-d commenced with the Rh2[5(S)-MEPY]4-catalyzed cyclization of the allylic diazoesters 16a-d to give the cyclopropyl lactones 17a-d in high enantiomeric excess. Opening of the lactone moiety using the Weinreb protocol and straightforward refunctionalization of the intermediate amides 18a-d gave 20a-d. A similar sequence of reactions was used to prepare the N-methyl-2-pyridyl analogue 28. Coupling of 20a-d and 28 with the known diamino diol 22 delivered 9-13. Pseudopeptides 9-12 were found to be competitive inhibitors of wild-type HIV-1 protease in biological assays having Kis of 0.31-0.35 nM for 9, 0.16-0.21 nM for 10, 0.47 nM for 11, and 0.17 nM for 12; these inhibitors were thus approximately equipotent to the known inhibitor 14(IC50 = 0.22 nM) from which they were derived. On the other hand 13 (Ki = 80 nM) was a weaker inhibitor than its analogue 15 (Ki = 0.11 nM). The solution structures of 9 and 10 were analyzed by NMR spectroscopy and simulated annealing procedures that included restraints derived from homo- and heteronuclear coupling constants and NOEs; because of the molecular symmetry of9 and 10, a special protocol to treat the NOE data was used. The final structure was checked by restrained and free molecular dynamic calculations using an explicit DMSO solvent box. The preferred solution conformations of 9 and 10 are extended structures that closely resemble the three-dimensional structure of 10 bound to HIV-1 protease as determined by X-ray crystallographic analysis of the complex. This work convincingly demonstrates that extended structures of peptides may be stabilized by the presence of substituted cyclopropanes that serve as peptide replacements. Moreover, the linear structure enforced in solution by the two cyclopropane rings in the pseudopeptides 9-12 appears to correspond closely to the biologically active conformation of the more flexible inhibitors 14 and 15. The present work, which is a combination of medicinal, structural, and quantum chemistry, thus clearly establishes that cyclopropanes may be used as structural constraints to reduce the flexibility of linear pseudopeptides and to help enforce the biologically active conformation of such ligands in solution.

Molecular mechanisms of resistance: free energy calculations of mutation effects on inhibitor binding to HIV-1 protease.

The changes in the inhibitor binding constants due to the mutation of isoleucine to valine at position 84 of HIV-1 protease are calculated using molecular dynamics simulations. The calculations are done for three potent inhibitors--KNI-272, L-735,524 (indinavir or MK-639), and Ro 31-8959 (saquinavir). The calculations agree with the experimental data both in terms of an overall trend and in the magnitude of the resulting free energy change. HIV-1 protease is a homodimer, so each mutation causes two changes in the enzyme. The decrease in the binding free energy from each mutated side chain differs among the three inhibitors and correlates well with the size of the cavities induced in the protein interior near the mutated residue. The cavities are created as a result of a mutation to a smaller side chain, but the cavities are less than would be predicted from the wild-type structures, indicating that there is significant relaxation to partially fill the cavities.

All-atom models for the non-nucleoside binding site of HIV-1 reverse transcriptase complexed with inhibitors: a 3D QSAR approach.

Several molecular modeling techniques were used to generate an all-atom molecular model of a receptor binding site starting only from Ca atom coordinates. The model consists of 48 noncontiguous residues of the non-nucleoside binding site of HIV-1 reverse transcriptase and was generated using a congeneric series of nevirapine analogs as structural probes. On the basis of the receptor-ligand atom contacts, the program HINT was used to develop a 3D quantitative structure activity relationship that predicted the rank order of binding affinities for the series of inhibitors. Electronic profiles of the ligands in their docked conformations were characterized using electrostatic potential maps and frontier orbital calculations. These results led to the development of a 3D stereoelectronic pharmacophore which was used to construct 3D queries for database searches. A search of the National Cancer Institute's open database identified a lead compound that exhibited moderate antiviral activity.

Structure of HIV-1 protease with KNI-272, a tight-binding transition-state analog containing allophenylnorstatine.

BACKGROUND: HIV-1 protease (HIV PR), an aspartic protease, cleaves Phe-Pro bonds in the Gag and Gag-Pol viral polyproteins. Substrate-based peptide mimics constitute a major class of inhibitors of HIV PR presently being developed for AIDS treatment. One such compound, KNI-272, which incorporates allophenylnorstatine (Apns)-thioproline (Thp) in place of Phe-Pro, has potent antiviral activity and is undergoing clinical trials. The structure of the enzyme-inhibitor complex should lead to an understanding of the structural basis for its tight binding properties and provide a framework for interpreting the emerging resistance to this drug. RESULTS: The three-dimensional crystal structure of KNI-272 bound to HIV PR has been determined to 2.0 A resolution and used to analyze structure-activity data and drug resistance for the Arg8-->Gln and ILe84-->Val mutations in HIV PR. The conformationally constrained Apns-Thp linkage is favorably recognized in its low energy trans conformation, which results in a symmetric mode of binding to the active-site aspartic acids and also explains the unusual preference of HIV PR for the S, or syn, hydroxyl group of the Apns residue. The inhibitor recognizes the enzyme via hydrogen bonds to three bridging water molecules, including one that is coordinated directly to the catalytic Asp125 residue. CONCLUSIONS: The structure of the HIV PR/KNI-272 complex illustrates the importance of limiting the conformational degrees of freedom and of using protein-bound water molecules for building potent inhibitors. The binding mode of HIV PR inhibitors can be predicted from the stereochemical relationship between adjacent hydroxyl-bearing and side chain bearing carbon atoms of the P1 substituent. Our structure also provides a framework for designing analogs targeted to drug-resistant mutant enzymes.

Quantitative evaluation of hydration thermodynamics with a continuum model.

We attempt to analyze whether experimental entropies, enthalpies and free energies of hydration of small uncharged molecules can be quantitatively rationalized with a continuum model including a classical reaction field formalism. We find that a simple proportionality to accessible surface with five different atom types allows satisfactory (within 1-1.5 kcal/mol) reproduction of hydration entropies (T delta S) of over 40 solutes. The agreement with experiment can possibly be improved if proximity effects and configurational contributions to transfer entropies are taken into account. In calculations of hydration enthalpies a reasonable agreement with experimental data can be obtained only when solute polarizability is taken into account. Electrostatic contributions to calculated hydration enthalpies exhibit strong dependencies on both the magnitude and the direction of molecular dipole moments. We demonstrate that for 20 molecules with experimentally measured vacuum dipole moments density functional calculations with DZVPD basis set including diffuse functions on d-orbitals allows prediction of experimental dipole moments within 0.1 D. At a fixed direction of the molecular dipole moment, mu, the electrostatic component of hydration enthalpy varies as mu 2. Thus an uncertainty of 0.1 D corresponds to uncertainties of 0.5-0.7 kcal/mol in hydration enthalpies of most small dipolar solutes. A 30 degree change in the direction of the molecular dipole together with the corresponding change in the quadrupole moment can result in a change of hydration enthalpy of 3 kcal/mol. Changes in the quadrupole moment alone can result in hydration enthalpy changes of over 1 kcal/mol. Representations of multipole expansions by point charges on nuclei fitted to molecular electrostatic potentials cannot accurately reproduce all these factors. Use of such point charges in calculations of hydration enthalpies predictably leads to discrepancies with experiment of approximately 3 kcal/mol for some solutes. However, errors in hydration enthalpies and hydration entropies are usually compensating leading in most cases to agreement between calculated and experimental free energies of hydration within 1.5 kcal/mol.