Antipneumococcal activity of neuraminidase inhibiting artocarpin.

Streptococcus (S.) pneumoniae is a major cause of secondary bacterial pneumonia during influenza epidemics. Neuraminidase (NA) is a virulence factor of both pneumococci and influenza viruses. Bacterial neuraminidases (NAs) are structurally related to viral NA and susceptible to oseltamivir, an inhibitor designed to target viral NA. This prompted us to evaluate the antipneumococcal potential of two NA inhibiting natural compounds, the diarylheptanoid katsumadain A and the isoprenylated flavone artocarpin. Chemiluminescence, fluorescence-, and hemagglutination-based enzyme assays were applied to determine the inhibitory efficiency (IC(50) value) of the tested compounds towards pneumococcal NAs. The mechanism of inhibition was studied via enzyme kinetics with recombinant NanA NA. Unlike oseltamivir, which competes with the natural substrate of NA, artocarpin exhibits a mixed-type inhibition with a Ki value of 9.70 µM. Remarkably, artocarpin was the only NA inhibitor (NAI) for which an inhibitory effect on pneumococcal growth (MIC: 0.99-5.75 µM) and biofilm formation (MBIC: 1.15-2.97 µM) was observable. In addition, we discovered that the bactericidal effect of artocarpin can reduce the viability of pneumococci by a factor of >1000, without obvious harm to lung epithelial cells. This renders artocarpin a promising natural product for further investigations.

Significance of ligand tails for interaction with the minor groove of B-DNA.

Minor groove binding ligands are of great interest due to their extraordinary importance as transcription controlling drugs. We performed three molecular dynamics simulations of the unbound d(CGCGAATTCGCG)(2) dodecamer and its complexes with Hoechst33258 and Netropsin. The structural behavior of the piperazine tail of Hoechst33258, which has already been shown to be a contributor in sequence-specific recognition, was analyzed. The simulations also reveal that the tails of the ligands are able to influence the width of the minor groove. The groove width is even sensitive for conformational transitions of these tails, indicating a high adaptability of the minor groove. Furthermore, the ligands also exert an influence on the B(I)/B(II) backbone conformational substate behavior. All together these results are important for the understanding of the binding process of sequence-specific ligands.

Complex of B-DNA with polyamides freezes DNA backbone flexibility.

The development of sequence-specific minor groove binding ligands is a modern and rapidly growing field of research because of their extraordinary importance as transcription-controlling drugs. We performed three molecular dynamics simulations in order to clarify the influence of minor groove binding of two ImHpPyPy-beta-Dp polyamides to the d(CCAGTACTGG)(2) decamer in the B-form. This decamer contains the recognition sequence for the trp repressor (5'-GTACT-3'), and it was investigated recently by X-ray crystallography. On one hand we are able to reproduce X-ray-determined DNA--drug contacts, and on the other hand we provide new contact information which is important for the development of potential ligands. The new insights show how the beta-tail of the polyamide ligands contributes to binding. Our simulations also indicate that complexation freezes the DNA backbone in a specific B(I) or B(II) substate conformation and thus optimizes nonbonded contacts. The existence of this distinct B(I)/B(II) substate pattern also allows the formation of water-mediated contacts. Thus, we suggest the B(I) <==> B(II) substate behavior to be an important part of the indirect readout of DNA.

The reaction rate constant of chlorine nitrate hydrolysis.

The first-order rate constant for the decomposition of chlorine nitrate (ClONO2) by water in a cyclic 1:3 complex at stratospheric temperatures is shown to be close to the values for the hydrolysis rate coefficient of chlorine nitrate on an ice surface determined in the laboratory. On the other hand the rate constants calculated for the cyclic 1:1 and 1:2 complexes are much lower than the experimental results. From the mechanistic point of view the reaction is found to be similar to a SN2 mechanism and coupled with water-mediated proton transfer in accordance with the intriguing findings of Bianco and Hynes [R. Bianco, J. T. Hynes. J. Phys. Chem. A 1998, 102, 309-314]. The function of additional water molecules is to act as a catalyst, that is, to accelerate the hydrolysis process. Quantum-mechanical tunneling is negligible above 125 K in the 1:3 complex and above 175 K in the 1:2 complex. At temperatures below these limits all involved protons tunnel through the barrier at energies at least 5 kcalmol(-1) below the barrier-top in a concerted, but asynchronous manner.

Prediction of the structure of human Janus kinase 2 (JAK2) comprising the two carboxy-terminal domains reveals a mechanism for autoregulation.

The structure of human Janus kinase 2 (JAK2) comprising the two C-terminal domains (JH1 and JH2) was predicted by application of homology modelling techniques. JH1 and JH2 represent the tyrosine kinase and tyrosine kinase-like domains, respectively, and are crucial for function and regulation of the protein. A comparison between the structures of the two domains is made and structural differences are highlighted. Prediction of the relative orientation of JH1 and JH2 was aided by a newly developed method for the detection of correlated amino acid mutations. Analysis of the interactions between the two domains led to a model for the regulatory effect of JH2 on JH1. The predictions are consistent with available experimental data on JAK2 or related proteins and provide an explanation for inhibition of JH1 tyrosine kinase activity by the adjacent JH2 domain.

Influence of netropsin's charges on the minor groove width of d(CGCGAATTCGCG)2.

The exact understanding of the interaction of minor groove binding drugs with DNA is of interest due to their importance as transcription controlling drugs. In this study we performed four molecular dynamics simulations, one of the uncomplexed d(CGCGAATTCGCG)(2) dodecamer and three simulations of the DNA complexed with the minor groove binder netropsin. The charged guanidinium and amidinium ends of the small ligand were in one simulation formally uncharged, in the second one normally charged, and in the third simulation we doubled the charges of the two ends. So we are able to filter out the influence the charges exert on the DNA structure. The positive charges reduce the width of the minor groove showing that charges are able to modify the groove width by charge neutralization of the negative phosphate groups. The quality of the used force field was successfully tested by comparing the results of the uncomplexed dodecamer with already reported NMR and x-ray studies. Thus our simulations should be able to describe the minor groove width of DNA in a correct manner underlying the validity of the results.

Exocyclic groups in the minor groove influence the backbone conformation of DNA.

Exocyclic groups in the minor groove of DNA modulate the affinity and positioning of nucleic acids to the histone protein. The addition of exocyclic groups decreases the formation of this protein-DNA complex, while their removal increases nucleosome formation. On the other hand, recent theoretical results show a strong correlation between the B(I)/B(II) phosphate backbone conformation and the hydration of the grooves of the DNA. We performed a simulation of the d(CGCGAATTCGCG)2 Drew Dickerson dodecamer and one simulation of the d(CGCIAATTCGCG)2 dodecamer in order to investigate the influence of the exocyclic amino group of guanine. The removal of the amino group introduces a higher intrinsic flexibility to DNA supporting the suggestions that make the enhanced flexibility responsible for the enlarged histone complexation affinity. This effect is attributed to changes in the destacking interactions of both strands of the DNA. The differences in the hydration of the minor groove could be the explanation of this flexibility. The changed hydration of the minor groove also leads to a different B(I)/B(II) substate pattern. Due to the fact that the histone preferentially builds contacts with the backbone of the DNA, we propose an influence of these B(I)/B(II) changes on the nucleosome formation process. Thus, we provide an additional explanation for the enhanced affinity to the histone due to removal of exocyclic groups. In terms of B(I)/B(II) we are also able to explain how minor groove binding ligands could affect the nucleosome assembly without disrupting the structure of DNA.

Toward elimination of discrepancies between theory and experiment: the rate constant of the atmospheric conversion of SO3 to H2SO4.

The hydration rate constant of sulfur trioxide to sulfuric acid is shown to depend sensitively on water vapor pressure. In the 1:1 SO3-H2O complex, the rate is predicted to be slower by about 25 orders of magnitude compared with laboratory results [Lovejoy, E. R., Hanson, D. R. & Huey, L. G. (1996) J. Phys. Chem. 100, 19911-19916; Jayne, J. T., Poschl, U., Chen, Y.-m., Dai, D., Molina, L. T., Worsnop, D. R., Kolb, C. E. & Molina, M. J. (1997) J. Phys. Chem. A 101, 10000-10011]. This discrepancy is removed mostly by allowing a second and third water molecule to participate. An asynchronous water-mediated double proton transfer concerted with the nucleophilic attack and a double proton transfer accompanied by a transient H3O+ rotation are predicted to be the fastest reaction mechanisms. Comparison of the predicted negative apparent "activation" energies with the experimental finding indicates that in our atmosphere, different reaction paths involving two and three water molecules are taken in the process of forming sulfate aerosols and consequently acid rain.

Elbow flexibility and ligand-induced domain rearrangements in antibody Fab NC6.8: large effects of a small hapten.

Four 700-ps molecular dynamics simulations were carried out to analyze the structural dynamics of the antigen-binding antibody fragment NC6.8, which is known to exhibit large structural changes upon complexation. The first simulation was started from the x-ray structure of the uncomplexed Fab and produced trajectory averages that closely match the crystallographic results. It allowed assessment of the flexibility of the Fab, revealing an elbow motion of the variable domains with respect to the constant domains. The second simulation was started from the uncomplexed x-ray structure after insertion of the ligand into the binding site. This perturbation resulted in a significantly altered trajectory, with quaternary structural changes corresponding in many aspects to the experimental differences between complexed and uncomplexed state. The observed trend toward a smaller elbow angle and a higher flexion of the H-chain could also be seen in the third simulation, which was started from the x-ray structure of the complex. The changes were revealed to be a clear consequence of the complexation with the ligand because in the fourth simulation (started from the experimental complex structure after removal of the hapten) the Fab remained close to its initial structure. Analyses of the quaternary structure and the binding site of Fab NC6.8 are presented for all four simulations, and possible interpretations are discussed.

Automated docking of ligands to antibodies: methods and applications.

Many approaches to studying protein-ligand interactions by computational docking are currently available. Given the structures of a protein and a ligand, the ultimate goal of all docking methods is to predict the structure of the resulting complex. This requires a suitable representation of molecular structures and properties, search algorithms to efficiently scan the configuration space for favorable interaction geometries, and accurate scoring functions to evaluate and rank the generated orientations. For many of the available methods, tests on experimentally known antibody-antigen or antibody-hapten complexes have appeared in the literature. In addition, some of them have been used in predictive studies on antibody-ligand interactions to provide structural insights where adequate experimental information is missing. The AutoDock program is presented as example of a method for flexibly docking ligands to antibodies. Applying parameters of the second-generation AMBER force field, three antibody-hapten complexes (AN02, DB3, NC6.8) are used as new test cases to analyze the ability of the method to reproduce experimental findings. The X-ray structures could be reconstituted and the corresponding solutions were ranked with best energy score in all cases. Docking to the free instead of the complexed NC6.8 structure indicated the limits of the rigid protein treatment, although fairly good guesses about the location of the binding site and the contact residues could still be obtained if conformational flexibility was allowed at least in the ligand.

Application of multivariate data analysis methods to comparative molecular field analysis (CoMFA) data: proton affinities and pKa prediction for nucleic acids components.

Multivariate data analysis methods (Principal Component Analysis (PCA) and Partial Least Squares (PLS)) are applied to the analysis of the CoMFA (Comparative Molecular Field Analysis) data for several nucleic acids components. The data set includes nitrogenated bases, nucleosides, linear nucleotides, 3', 5'-cyclic nucleotides and oligonucleotides. PCA is applied to study the structure of the CoMFA data and to detect possible outliers in the data set. PLS is applied to correlate the CoMFA data with either calculated AM1 proton affinities or with experimental pKa values. The possibility of making a prediction of pKa values directly from 3D structures of the monomers for polynucleotides is also shown. The influence of the superposition criteria and of conformational changes along the glycosidic bond on the pKa prediction are studied as well.

Helix morphology changes in B-DNA induced by spontaneous B(I)<==>B(II) substrate interconversion.

Investigations of spontaneous, i.e. not forced, B-DNA's B(I)<==>B(II) substate transitions are carried out on the d(CGCGAATTCGCG)2 EcoRI dodecamer sequence using Molecular Dynamics Simulations. Analysis of the resulting transition processes with respect to the backbone angles reveals concerted changes not only for backbone angles epsilon, zeta, and beta, but also for the 5'-delta and 5'-chi angles. For alpha and delta inside the interconverting base step, a change is seen in short lived B(II) conformers. With respect to base morphology distinct changes are observed for buckle, propeller twist, shift, roll and twist, as well as x-displacement and tip. The base mainly involved in the changes is identified as the base preceding the interconverting phosphate. Altogether single B(I)<==>B(II) interconversions result only in local distortions represented by the larger spread of most parameters. Comparison of the atomic positional fluctuations derived from the simulation with those obtained from the static X-ray structure results in striking similarities.

Unexpected BII conformer substate population in unoriented hydrated films of the d(CGCGAATTCGCG)2 dodecamer and of native B-DNA from salmon testes.

Conformational substates of B-DNA had been observed so far in synthetic oligonucleotides but not in naturally occurring highly polymeric B-DNA. Our low-temperature experiments show that native B-DNA from salmon testes and the d(CGCGAATTCGCG)2 dodecamer have the same BI and BII substates. Nonequilibrium distribution of conformer population was generated by quenching hydrated unoriented films to 200 K, and isothermal structural relaxation toward equilibrium by interconversion of substates was followed by Fourier transform infrared spectroscopy. BI interconverts into BII on isothermal relaxation at 200 K, whereas on slow cooling from ambient temperature, BII interconverts into BI. Our estimation of the dodecamer's BI-to-BII conformer substate population by curve resolution of the symmetrical stretching vibration of the ionic phosphate is 2.4 +/- 0.5 to 1 at 200 K, and it is 1.3 +/- 0.5 to 1 between 270 and 290 K. Pronounced spectral changes upon BI-to-BII interconversion are consistent with base destacking coupled with migration of water from ionic phosphate toward the phosphodiester and sugar moieties. Nonspecific interaction of proteins with the DNA backbone could become specific by induced-fit-type interactions with either BI or BII backbone conformations. This suggests that the BI-to-BII substate interconversion could be a major contributor to the protein recognition process.

Ligand binding by antibody IgE Lb4: assessment of binding site preferences using microcalorimetry, docking, and free energy simulations.

Antibody IgE Lb4 interacts favorably with a large number of different compounds. To improve the current understanding of the structural basis of this vast cross-reactivity, the binding of three dinitrophenyl (DNP) amino acids (DNP-alanine, DNP-glycine, and DNP-serine) is investigated in detail by means of docking and molecular dynamics free energy simulations. Experimental binding energies obtained by isothermal titration microcalorimetry are used to judge the results of the computational studies. For all three ligands, the docking procedure proposes two plausible subsites within the binding region formed by the antibody CDR loops. By subsequent molecular dynamics simulations and calculations of relative free energies of binding, one of these subsites, a tyrosine-surrounded pocket, is revealed as the preferred point of complexation. For this subsite, results consistent with experimental observations are obtained; DNP-glycine is found to bind better than DNP-serine, and this, in turn, is found to bind better than DNP-alanine. The suggested binding mode makes it possible to explain both the moderate binding affinity and the differences in binding energy among the three ligands.

Comparative molecular field analysis of artemisinin derivatives: ab initio versus semiempirical optimized structures.

Based on the belief that structural optimization methods, producing structures more closely to the experimental ones, should give better, i.e. more relevant, steric fields and hence more predictive CoMFA models, comparative molecular field analyses of artemisinin derivatives were performed based on semiempirical AM1 and HF/3-21G optimized geometries. Using these optimized geometries, the CoMFA results derived from the HF/3-21G method are found to be usually but not drastically better than those from AM1. Additional calculations were performed to investigate the electrostatic field difference using the Gasteiger and Marsili charges, the electrostatic potential fit charges at the AM1 level, and the natural population analysis charges at the HF/3-21G level of theory. For the HF/3-21G optimized structures no difference in predictability was observed, whereas for AM1 optimized structures such differences were found. Interestingly, if ionic compounds are omitted, differences between the various HF/3-21G optimized structure models using these electrostatic fields were found.

Ligand-induced domain movement in an antibody Fab: molecular dynamics studies confirm the unique domain movement observed experimentally for Fab NC6.8 upon complexation and reveal its segmental flexibility.

Two molecular dynamics simulations were carried out for the antibody Fab NC6.8, both with and without the guanidinium sweetener ligand NC174, in order to assess the segmental flexibility as well as the conformational changes upon ligand binding. Trajectory analyses of the simulation of the uncomplexed Fab suggest low-amplitude motions of the Ig domains with respect to each other, most clearly reflected by a periodic alteration of the elbow angle within a range of 11 degrees. Upon insertion of the hapten into the binding site, the quaternary structure of the Fab exhibits considerable rearrangements: the elbow angle changes by almost 30 degrees, the light chain is elongated and the heavy chain becomes more flexed. Comparison with experiment reveals some interesting agreements with X-ray crystallographic results published previously.

Quantitative analysis of the structural requirements for blockade of the N-methyl-D-aspartate receptor at the phencyclidine binding site.

Blockade of the N-methyl-D-aspartate receptor by uncompetitive antagonists has implications for symptomatic and neuroprotective therapy of various neuropsychiatric diseases. Since the three-dimensional (3D) structure of this ion channel is unknown, the structural requirements for uncompetitive inhibition were investigated by application of a three-step strategy: At first, Ki values were measured for a number of structurally diverse compounds at the phencyclidine (PCP) binding site in postmortem human frontal cortex. Second, a pharmacophore model was developed for this set of molecules employing a mathematical method called graph theory. The resulting pharmacophore provided a very good explanation for the ability of structurally diverse compounds to bind to the same binding site. Using the experimental data and the pharmacophore as a basis for the third step, a three-dimensional quantitative structure-activity relationship (3D-QSAR) applying comparative molecular field analysis (CoMFA) was performed. The QSAR proved to be highly consistent and showed good predictiveness for several additional molecules. The results give a deeper insight into the structural requirements for effective NMDA receptor antagonism and offer the opportunity for improved drug design. The study represents the first quantitative 3D-QSAR model for NMDA receptor blockade, and it comprises structurally very different molecules. That the alignment for a highly consistent CoMFA is based on an automated 3D pharmacophore analysis has important methodological implications.

Comparative molecular field analysis of haptens docked to the multispecific antibody IgE(Lb4)

Using comparative molecular field analysis (CoMFA), three-dimensional quantitative structure-activity relationships were developed for 27 haptens which bind to the monoclonal antibody IgE(Lb4). In order to obtain an alignment for these structurally very diverse antigens, the compounds were docked to a previously modeled receptor structure using the automated docking program AUTODOCK (Goodsell, D.S.; Olson, A.J. Proteins: Struct., Funct., Genet. 1990, 8, 195-202). Remarkably, this alignment method yielded highly consistent QSAR models, as indicated by the corresponding cross-validated r2 values (0.809 for a model with carbon as probe atom, 0.773 for a model with hydrogen as probe atom). Conventional alignment failed in providing a basis for self-consistent CoMFAs. Amino acids Tyr H 50, Tyr H 52, and Trp H 95 of the receptor appeared to be of crucial importance for binding of various antigens. These findings are consistent with earlier considerations of aromatic residues being responsible for the multispecificity of certain immunoglobulins.

Comparative docking studies on ligand binding to the multispecific antibodies IgE-La2 and IgE-Lb4.

A large comparative study is presented in which the binding of approximately 30 different ligands to two IgE antibodies (La2 and Lb4) is analyzed by means of an automated-docking procedure based on simulated annealing. The method is able to reproduce experimentally verified binding orientations, as shown by application to the Ig-AN02-hapten complex. The main address of the study is to investigate the concept of antibody multispecificity. Problems and usefulness of docking in this context are discussed. The results indicate reasons for multispecific binding properties and how they can be understood from the topology of the binding site. Though similar in general behaviour, the two antibodies show interesting differences in their binding characteristics. The binding sites of both antibodies are described and the main interacting residues revealed.

Prediction of IgE(Lb4)-ligand complex structures by automated docking.

A mouse monoclonal anti-2,4,6-trinitrophenyl IgE (clone Lb4) was screened with a random set of over 2000 compounds, and several ligands were found to bind with affinities comparable to that of the immunizing hapten (KD in the microM range). An automated docking algorithm was used for the prediction of complex structures formed by 2,4-dinitrophenyl (DNP) and non-DNP ligands in the fragment variable region of IgE(Lb4). All ligands were found to dock in an L-shaped cavity of 15 x 16 x 10 A, surrounded by complementary-determining regions L1, L3, H2 and H3. The ligands were found to occupy the same binding site in different orientations. For rigid ligands the most stable orientation could be predicted with high probability, based on the calculated energy of binding and the occurrence frequencies of identical complexes obtained by repeated simulations. The localization of a flexible ligand (cycrimine-R) was more ambiguous, but it still docked in the same site. The results support a model for heteroligating antibody (Ab) binding sites, where different ligands utilize the total set of available contacts in different combinations. It is suggested that although pseudoenergies calculated by the docking algorithm do not correlate with experimentally measured binding energies, the screening-and-docking procedure can be useful for the mapping of Ab and other receptor binding sites ligating small molecules.