A comparison of structural and dynamic properties of different simulation methods applied to SH3.

The dynamic and static properties of molecular dynamics simulations using various methods for treating solvent were compared. The SH3 protein domain was chosen as a test case because of its small size and high surface-to-volume ratio. The simulations were analyzed in structural terms by examining crystal packing, distribution of polar residues, and conservation of secondary structure. In addition, the "essential dynamics" method was applied to compare each of the molecular dynamics trajectories with a full solvent simulation. This method proved to be a powerful tool for the comparison of large concerted atomic motions in SH3. It identified methods of simulation that yielded significantly different dynamic properties compared to the full solvent simulation. Simulating SH3 using the stochastic dynamics algorithm with a vacuum (reduced charge) force field produced properties close to those of the full solvent simulation. The application of a recently described solvation term did not improve the dynamic properties. The large concerted atomic motions in the full solvent simulation as revealed by the essential dynamics method were analyzed for possible biological implications. Two loops, which have been shown to be involved in ligand binding, were seen to move in concert to open and close the ligand-binding site.     

Virtual screening to enrich a compound collection with CDK2 inhibitors using docking, scoring, and composite scoring models

Docking programs can generate subsets of a compound collection with an increased percentage of actives against a target (enrichment) by predicting their binding mode (pose) and affinity (score), and retrieving those with the highest scores. Using the QXP and GOLD programs, we compared the ability of six single scoring functions (PLP, Ligscore, Ludi, Jain, ChemScore, PMF) and four composite scoring models (Mean Rank: MR, Rank-by-Vote: Vt, Bayesian Statistics: BS and PLS Discriminant Analysis: DA) to separate compounds that are active against CDK2 from inactives. We determined the enrichment for the entire set of actives (IC50 < 10 microM) and for three activity subsets. In all cases, the enrichment for each subset was lower than for the entire set of actives. QXP outperformed GOLD at pose prediction, but yielded only moderately better enrichments. Five to six scoring functions yielded good enrichments with GOLD poses, while typically only two worked well with QXP poses. For each program, two scoring functions generally performed better than the others (Ligscore2 and Ludi for GOLD; QXP and Jain for QXP). Composite scoring functions yielded better results than single scoring functions. The consensus approaches MR and Vt worked best when separating micromolar inhibitors from inactives. The statistical approaches BS and DA, which require training data, performed best when distinguishing between low and high nanomolar inhibitors. The key observation that all hit rate profiles for all four activity intervals for all scoring schemes for both programs are significantly better than random, is evidence that docking can be successfully applied to enrich compound collections.     

Inhibition of protein-protein interactions: the discovery of druglike beta-catenin inhibitors by combining virtual and biophysical screening

The interaction between beta-catenin and Tcf family members is crucial for the Wnt signal transduction pathway, which is commonly mutated in cancer. This interaction extends over a very large surface area (4800 A(2)), and inhibiting such interactions using low molecular weight inhibitors is a challenge. However, protein surfaces frequently contain "hot spots," small patches that are the main mediators of binding affinity. By making tight interactions with a hot spot, a small molecule can compete with a protein. The Tcf3/Tcf4-binding surface on beta-catenin contains a well-defined hot spot around residues K435 and R469. A 17,700 compounds subset of the Pharmacia corporate collection was docked to this hot spot with the QXP program; 22 of the best scoring compounds were put into a biophysical (NMR and ITC) screening funnel, where specific binding to beta-catenin, competition with Tcf4 and finally binding constants were determined. This process led to the discovery of three druglike, low molecular weight Tcf4-competitive compounds with the tightest binder having a K(D) of 450 nM. Our approach can be used in several situations (e.g., when selecting compounds from external collections, when no biochemical functional assay is available, or when no HTS is envisioned), and it may be generally applicable to the identification of inhibitors of protein-protein interactions.     

O2 versus N2O respiration in a continuous microbial enrichment

Despite its ecological importance, essential aspects of microbial N₂O reduction—such as the effect of O₂ availability on the N₂O sink capacity of a community—remain unclear. We studied N₂O vs. aerobic respiration in a chemostat culture to explore (i) the extent to which simultaneous respiration of N₂O and O₂ can occur, (ii) the mechanism governing the competition for N₂O and O₂, and (iii) how the N₂O-reducing capacity of a community is affected by dynamic oxic/anoxic shifts such as those that may occur during nitrogen removal in wastewater treatment systems. Despite its prolonged growth and enrichment with N₂O as the sole electron acceptor, the culture readily switched to aerobic respiration upon exposure to O₂. When supplied simultaneously, N₂O reduction to N₂ was only detected when the O₂ concentration was limiting the respiration rate. The biomass yields per electron accepted during growth on N₂O are in agreement with our current knowledge of electron transport chain biochemistry in model denitrifiers like Paracoccus denitrificans. The culture’s affinity constant (K_S) for O₂ was found to be two orders of magnitude lower than the value for N₂O, explaining the preferential use of O₂ over N₂O under most environmentally relevant conditions.

     

Molecular Dynamics Simulations of some Small Organic Molecules

Abstract

Free energy differences between different conformers of D-ribofuranose, L-malic acid and meso-tartaric acid in solution were calculated using Molecular Dynamics simulations. In case of ribose the α → β transition was studied. For the acids attention was focussed on the transitions between the three possible staggered conformers with respect to the central C-C bond. In all cases a thermodynamic integration method was employed to evaluate the free energy difference. The use of an alternative technique, umbrella sampling, for ribose did not give promising results.

It was shown that one needs a fairly accurate picture of the accessible conformational space in case of flexible molecules like the ones considered here before one can determine meaningful free energy differences. Large hysteresis effects between forward and reverse simulated transitions were observed, but contrary to the general belief they are no direct measure of the accuracy of the calculated ΔG values. In all cases the ΔG values resulting from the simulations and from NMR experiments agree within the, considerable, error limits and for the different forms of D-ribose, L-malic acid and L-tartaric acid the relative order of their populations is also correctly reproduced.     

New triple-helical model for the shaft of the adenovirus fibre.

The adenovirus fibre is a trimeric protein with a globular head on a long thin shaft that projects from the twelve fivefold vertices of the virion. The shaft region of the fibre primary sequence has a unique pseudo-repeating motif of 15 residues. Using constraints derived from sequence analysis, the trimeric nature of the fibre, the experimental determination of the shaft length and general knowledge about protein structure, an atomic model of the fibre shaft has been constructed by computer modelling techniques. In the final model the three monomers form a left-handed triple-helical structure with threefold symmetry and with successive 15-residue repeats on the same chain related by an axial rise of 13.1 A and a left-handed azimuthal rotation of close to 300 degrees. Three threefold related beta-sheets with short strands are formed by inter-monomer main-chain hydrogen bonds giving rise to superhelical ribbons covering the surface of the shaft. The model satisfies criteria of extensive hydrogen bonding, reasonable backbone torsion angles, burial of most hydrophobic residues and good packing of the hydrophobic core. Furthermore, the model is consistent with the observed shaft length of about 290 A and its calculated X-ray fibre diffraction patterns shows the characteristic features found in the diffraction pattern of crystals of fibre, notably layer lines with a spacing of about 1/26 A-1 and strong meridional intensity at 1/4.4 A-1.     

Chronic Allopurinol Treatment during the Last Trimester of Pregnancy in Sows: Effects on Low and Normal Birth Weight Offspring

Low-birth-weight (LBW) children are born with several risk factors for disease, morbidity and neonatal mortality, even if carried to term. Placental insufficiency leading to hypoxemia and reduced nutritional supply is the main cause for LBW. Brain damage and poor neurological outcome can be the consequence. LBW after being carried to term gives better chances for survival, but these children are still at risk for poor health and the development of cognitive impairments. Preventive therapies are not yet available. We studied the risk/efficacy of chronic prenatal treatment with the anti-oxidative drug allopurinol, as putative preventive treatment in piglets. LBW piglets served as a natural model for LBW. A cognitive holeboard test was applied to study the learning and memory abilities of these allopurinol treated piglets after weaning. Preliminary analysis of the plasma concentrations in sows and their piglets suggested that a daily dose of 15 mg.kg⁻¹ resulted in effective plasma concentration of allopurinol in piglets. No adverse effects of chronic allopurinol treatment were found on farrowing, birth weight, open field behavior, learning abilities, relative brain, hippocampus and spleen weights. LBW piglets showed increased anxiety levels in an open field test, but cognitive performance was not affected by allopurinol treatment. LBW animals treated with allopurinol showed the largest postnatal compensatory body weight gain. In contrast to a previous study, no differences in learning abilities were found between LBW and normal-birth-weight piglets. This discrepancy might be attributable to experimental differences. Our results indicate that chronic prenatal allopurinol treatment during the third trimester of pregnancy is safe, as no adverse side effects were observed. Compensatory weight gain of treated piglets is a positive indication for the chronic prenatal use of allopurinol in these animals. Further studies are needed to assess the possible preventive effects of allopurinol on brain functions in LBW piglets.     

Interchain cysteine bridges control entry of progesterone to the central cavity of the uteroglobin dimer.

The progesterone-binding protein uteroglobin has been expressed in Escherichia coli in an unfused, soluble form. Like mature uteroglobin from rabbit endometrium (UG), the E.coli produced uteroglobin (UG1) dimerizes in vitro, forms an antiparallel dimer with Cys3-Cys69' and Cys69-Cys3' disulfide bonds and binds progesterone under reducing conditions. In order to analyze the dimerization and the reduction dependence of progesterone binding in more detail, we separately replaced cysteine 3 and cysteine 69 by serines. Under reducing conditions, both uteroglobin variants (UG1-3Ser and UG1-69Ser) bind progesterone with the same affinity as the wild-type suggesting that both cysteine residues are not directly involved in progesterone binding. In contrast to the wild-type protein, both cysteine variants also bind progesterone with high affinity in the absence of reducing agents. In addition, UG1-3Ser and UG1-69Ser both form covalently linked homodimers. Thus, unnatural Cys69-69' and Cys3-3' disulfide bonds exist in UG1-3Ser and UG1-69Ser, respectively. These data together with computer models based on X-ray diffraction data strongly support the idea that progesterone reaches its binding site located in an internal hydrophobic cavity via a hydrophobic tunnel along helices 1 and 4. Under non-reducing conditions the tunnel is closed by two disulfide bridges (Cys3-Cys69' and Cys69-Cys3') that lie in the most flexible region of the dimer. Reduction or replacement of a cysteine residue enables conformational changes that open the channel allowing progesterone to enter.     

Monoclonal antibodies for immunohistochemical labeling of immunocompetent cells in frozen sections of rhesus monkey tissues.

Twenty-eight anti-human and two rhesus specific monoclonal antibodies (MAbs) were evaluated for immunohistochemical peroxidase staining reactivity in rhesus monkeys lymph nodes, thymus, liver, and skin allografts. Reactivity with the following antigens was assessed: MHC class I, II-DR, -DQ, and -DP antigens; leukocyte markers CD1, CD2, CD3, CD4, CD8, CD14, CD16, CD25, CD57; a proliferation associated nuclear antigen; interferon-gamma and tumor necrosis factor-alpha. Twenty-three MAbs proved to be suitable for immunohistochemical staining on frozen sections.     

(N-acyl-N-alkyl)glycyl borolysine analogs: A new class of potent thrombin inhibitors

In this report the structure-activity relationships of a series of novel (N-acyl-N-alkyl)glycyl borolysine thrombin inhibitors are described. This work culminates in the discovery of (N-3-phenylpropanoyl-N-phenethyl)glycyl borolysine (12j), a potent, orally active inhibitor with a binding conformation in which the N-phenethyl group occupies the aryl binding pocket of thrombin.