A possible binding path of ergosterol within elicitins revealed by molecular dynamics.

Elicitins, produced by most of the phytopathogenic fungi of the genus Phytophthora, provoke in the tobacco plant both remote leaf necrosis and the induction of a resistance against subsequent attack by various micro-organisms. The crystal structure of b-cryptogein (CRY), secreted by Phytophthora cryptogea, was previously reported as well as the first structure of a SCP/sterol complex, the ergosterol-complexed, mutated CRY (K13H). In K13H, the ergosterol molecule is encapsulated in a large internal hydrophobic cavity which is not present in CRY. This binding induces a minor conformational change in the protein structure. Molecular dynamics studies were undertaken to precise the structural behaviour of CRY and K13H with respect to the complexation of the ergosterol. Although it is not possible to simulate the entrance of the ergosterol in the protein, we assume that capture and release of the ligand possibly both occur following the same path. Our results show that, in the complex K13H, the ergosterol molecule is pushed towards the residue 13 which play a key role in the necrotic activity of the protein. It is likely that the polarity of residue 13, favouring the binding of the hydroxyl of the ligand, would be involved in the recognition of the sterol and in an optimisation of its orientation. Thus, in a first step, the molecule of ergosterol would be rotated around itself to a position which makes possible, in a second step, its translation to the internal cavity, as a key in a keyhole.

The uteroglobin fold.

Uteroglobin (UTG) forms a fascinating homodimeric structure that binds small- to medium-sized ligands through an internal hydrophobic cavity, located at the interface between the two monomers. Previous studies have shown that UTG fold is not limited to the UTG/CC10 family, whose sequence/structure relationships are highlighted here, but can be extended to the cap domain of Xanthobacter autotrophicus haloalkane dehalogenase. We show here that UTG fold is adopted by several other cap domains within the alpha/beta hydrolase family, making it a well-suited "geode" structure allowing it to sequester various hydrophobic molecules. Additionally, some data about a new crystal form of oxidized rabbit UTG are presented, completing previous structural studies, as well as results from molecular dynamics, suggesting an alternative way for the ligand to reach the internal cavity.

A theoretical and experimental study of two thiazole orange derivatives with single- and double-stranded oligonucleotides, polydeoxyribonucleotides and DNA.

The effect of interaction with DNA and oligonucleotides on the photophysical properties of two thiazole orange (TO) derivatives, with different side chains (-(CH2)3-N+(CH3)3 and -(CH2)6-I)) linked to the nitrogen of the quinoline ring of the thiazole orange, is presented here. The first one called TO-PRO1 is a commercially available dye, whereas the second one called TO-MET has been specially synthesized for further covalent binding to oligonucleotides with the aim of being used for specific in situ detection of biomolecular interactions. Both photophysical measurements and molecular calculations have been done to assess their possible mode of interaction with DNA. When dissolved in buffered aqueous solutions both derivatives exhibit very low fluorescence quantum yields of 8 x 10(-5) and 2 x 10(-4), respectively. However, upon binding to double-stranded DNA, large spectroscopic changes result and the quantum yield of fluorescence is enhanced by four orders of magnitude, reaching values up to phi F = 0.2 and 0.3, respectively, as a result of an intercalation mechanism between DNA base pairs. A modulation of the quantum yield is observed as a function of the base sequence. The two derivatives also bind with single-stranded oligonucleotides, but the fluorescence quantum yield is not so great as that when bound to double-stranded samples. Typical fluorescence quantum yields of 7 x 10(-3) to 3 x 10(-2) are observed when the dyes interact with short oligonucleotides, whereas the fluorescence quantum yield remains below 10(-2) when interacting with single-stranded oligonucleotides. This slight but significant quantum-yield increase is interpreted as a folding of the single strand around the dye, which reduces the internal rotation of the two heterocycles around the central methine bridge that links the two moieties of the dye. From these properties, it is proposed to link monomer covalently to oligonucleotides for the subsequent detection of target sequences within cells.

Conformational study based on simulated annealing and 1H NMR of peptides comprising the consensus sequence Arg5-Asp-Val-Arg-Gly9: effects of the substitution Ser 12 by Ala 12 or of 3 residues deletion at the N-terminus.

A conformational search by simulated annealing has been performed on two peptides derivated from the tetradecapeptide used to isolate the Xenopus laevis skin maturation RXVRG-endoprotease. The Ala 12 derivative, obtained by substitution in the hydrophobic C terminal fragment and the undecapeptide 4-14, obtained by deletion of an acidic rich tripeptide, were studied. No unique structure has been found for the tetradecapeptide Ala 12. This structural disorganization could explain the loss of activity of the endoprotease towards the substituted peptide. For the undecapeptide, two different models in accordance with the NMR data were found. The conformational differences between these two models are located in the consensus sequence and in each case an hairpin-like conformation is observed. These results could be related to the enhanced cleavage activity of the maturation enzyme. The obtained structures are also compared with those of the original tetradecapeptide.

Structure determination of a tetradecapeptide mimicking the RXVRG consensus sequence recognized by a Xenopus laevis skin endoprotease: an approach based on simulated annealing and 1H NMR.

The tetradecapeptide of sequence H-Asp-Val-Asp-Glu-Arg5-Asp-Val-Arg-Gly9-Phe-Ala-Ser-Phe-Leu- NH2 is recognized by a putative maturation endoprotease of the Xenopus laevis skin, which cleaves between Arg8 and Gly9. A conformational search has been performed on this peptide by simulated annealing calculations. Two different models in agreement with the NMR data were found. The conformational difference between the two types of model is located in the consensus sequence, i.e., from Arg5 to Gly9.

A conformational study of Lys-Arg-Asp-Ser and analogs, a series of potent antithrombotic peptides. An approach based on simulated annealing and 1H NMR.

Simulated annealing techniques were used to explore the conformational space of the potent antithrombotic peptide L.Lys-L.Arg-L.Asp-L.Ser (KRDS) and of two analogs: D.Lys-L.Arg-L.Asp-L.Ser (KDRDS), which is inactive, and L.Lys-L.Arg-L.Glu-L.Glu (KREE), which exhibits a strong biological activity. For each peptide, a set of initial conformations was generated and submitted to simulated annealing, including a heating to 1000 K followed by a cooling to 300 K. 200 resulting conformations of each compound were analyzed and classified according to the network of electrostatic interactions involving charged side chains and charged C- and N-terminal groups. A reduced number of conformational classes was obtained and conformations corresponding to predominant classes were found to be in qualitative agreement with structural parameters deduced from 1H NMR spectra. A comparison between the classes of the active and non active peptide was achieved. Some conformations were found to be specific of active peptides.

Theoretical study of the complexation of amphotericin B with sterols.

The aim of this present work was the study of the intermolecular complexes between amphotericin B (AmB) and either cholesterol or ergosterol. In such complexes the intermolecular interaction energy mainly proceeds from both Van der Waals and H-bonding (via water molecules) forces. Our calculations have shown that the Van der Waals forces slightly favor the AmB-ergosterol complex. Several relative positions of the sterol with regard to AmB lead to energy minima: sterol may be either in contact with the AmB polar head or repelled towards the end of the macrolide ring. It appeared that the role played by some water molecules was to maintain the sterol close to the AmB polar head.

Composite cylinder models of DNA: application to the electrostatics of the B-Z transition.

We develop and test a Poisson-Boltzmann model of the electrostatics of the B-Z transition of DNA. Starting from the detailed geometries of the two forms, we compute at each radius the fractions of DNA matter, of volume forbidden (for nonpoint-like ions), and of volume accessible to the center of ions. These radial distributions are incorporated in a composite cylinder model; availability to ions (porosity) and the dielectric constant at each radial distance are then obtained. The phosphate charge is distributed with cylindrical symmetry on two layers at the appropriate radial distances. The porous sheath, between the axis and the charge distribution, provides much more room for ions in B-DNA than in Z-DNA. By using previously developed methods, the Poisson-Boltzmann problem of such cylinders is easily solved. The computational load is small, so that results can be obtained for a large set of salt concentrations and for a number of ionic radii. The variation of the electrostatic free energy difference with salt concentration compares favorably with the experimental value (it is half as large). There is also qualitative agreement with experiments on supercoiled DNA, including a maximum of the free energy difference at submolar salt concentrations. The results for this cylinder with porous sheath are in line with those of the earlier simple planar model and of a plain cylinder with sheath, which is also presented here. They are thus insensitive to details of the model. They support the proposition that the main electrostatic feature of the B-Z transition is the better immersion of the B-DNA phosphates into the solution. They also give confidence in the validity of the Poisson-Boltzmann approach, despite the large salt concentrations involved. Prior studies using an approach based on the potential of mean force are discussed.

A unique structural feature of a phospholipase A2 is probed by molecular dynamics.

The unusual catalytic network, revealed by the crystal structure of one of the two phospholipases A2 (PLA2) from the venom of the crotalid A.p.piscivorus has been probed using molecular dynamics. The catalytic network has been remodeled to a conformation similar to that found in all other PLA2, and the modeled structure has been submitted to energy minimization and molecular dynamics simulation, to explore the conformational space of the network. The calculations have yielded a large reorganization of the catalytic network, which gets a conformation close to that of the crystal structure. These results suggest that the unusual catalytic network observed in the studied PLA2 is a structural feature of the protein and not a crystal artifact.

A simple explanation of the electrostatics of the B-to-Z transition of DNA.

Whereas the phosphates of B-DNA jut out into the solution, those of Z-DNA, being closer to DNA matter, are less subject to electrostatic screening by counterions. We present simple planar models of B- and Z-DNA that reflect these geometric features. The ionic strength dependence of the difference in the Poisson-Boltzmann electrostatic free energy of the models agrees with that measured by Pohl [Pohl, F. M. (1983) Cold Spring Harbor Symp. Quant. Biol. 47, 113-118]. This indicates that the electrostatics of the B-to-Z transition are primarily controlled by a qualitative geometrical difference and not by details of the DNA geometry or by complex electrostatic properties of the ionic solution.

Arginine-395 is required for efficient in vivo and in vitro aminoacylation of tRNAs by Escherichia coli methionyl-tRNA synthetase.

We have previously shown that the anticodon of methionine tRNAs contains the major recognition site required for aminoacylation of tRNAs by Escherichia coli methionyl-tRNA synthetase (MetRS) and have located part of the anticodon binding domain on the enzyme at a site close to Trp461 [Schulman, L. H., & Pelka, H. (1988) Science 242, 765-768; Ghosh, G., Pelka, H., & Schulman, L.H. (1990) Biochemistry 29, 2220-2225]. In order to gain information about other possible sites of contact between MetRS and its tRNA substrates, we have examined the effects of mutations at a series of positively charged residues on the surface of the C-terminal domain of the enzyme. Conversion of Arg356, Arg366, Arg380, or Arg453 to Gln had little or no effect on enzyme activity. Similarly, conversion of Lys402 or Lys439 to Asn failed to significantly alter aminoacylation activity. Conversion of Arg380 to Ala or Arg442 to Gln produced a 5-fold reduction in kcat/Km for aminoacylation of tRNAfMet, with no effect on methionine activation, indicating a possible minor role for these residues in interaction of the enzyme with the tRNA substrate. In contrast, mutation of a phylogenetically conserved residue, Arg395, to Gln increased the Km for aminoacylation of tRNAfMet about 30-fold and reduced kcat/Km by 25,000-fold. The mutant enzyme was also shown to be highly defective by its inability to complement a strain of E. coli having an altered chromosomal MetRS gene.(ABSTRACT TRUNCATED AT 250 WORDS)

Molecular dynamics simulations of phospholipases A2.

An extensive molecular dynamics study of phospholipases A2 from pancreatic bovine and Crotalus atrox venom has shown that the well-conserved homologous core of the phospholipases A2, including the so called catalytic network, is very stable during the course of the calculations. The fluctuations which occur are located in segments which have significantly different three-dimensional conformations in the two phospholipases A2 studied, suggesting that a particularly stable core conformation gives rise to a large homologous family of similar three-dimensional structure. The calcium ion, which exhibits a crucial structural role in the monomeric phospholipases A2, appears not to be required to stabilize the C.atrox dimer. Moreover, the behaviour of the dimeric structure during the dynamics raises the question of a possible dissociation of the two subunits into functional monomers.

Dimeric character of a basic phospholipase A2 from cobra venom: experimental and modelling study.

When it is gel filtered on Sephadex in the absence of calcium ions, basic phospholipase A2 from Naja nigricollis venom elutes as a dimer. In order to study the possibility of this dimerization from a structural point of view, three-dimensional models of both monomeric and dimeric N. nigricollis phospholipases A2 have been graphically built on the basis of homologies with the phospholipases A2 from pancreatic bovine and Crotalus atrox venom. The building of a dimeric model is made possible by the deletion of a particular loop of the bovine structure. The predicted models of N. nigricollis phospholipase A2 have been checked using molecular mechanics and molecular dynamics techniques according to a suitable protocol which has been developed starting from refined X-ray structures of phospholipases A2 as the test case. The observed stability of the dimeric model, in the absence of calcium, agrees with the hypothesis of the dimerization of the basic phospholipase A2. Particularly, Arg31, which replaces the hydrophobic residue present in pancreatic bovine and C.atrox venom phospholipases A2, contributes to this stability.

Molecular mechanics and dynamics of DNA-furocoumarin complexes: effect of the aromatization of the pyrone ring on the intercalation geometry.

Results of molecular mechanics and dynamics calculations on intercalation complexes of DNA with various furocoumarins (psoralen, angelicin, 7-methylpyrido[3,4-c]psoralen and 7-methylpyrido[4,3-c]psoralen) and their corresponding aromatized derivatives are presented. These calculations were undertaken with the aim to elucidate the roles of the pyrone and pyridine moieties in the interactions which tend to orient the furocoumarins and pyridopsoralens between DNA base pairs. It appears that the intercalation geometries are very similar for the furocoumarins and related aromatized compounds. Therefore the oxygen and nitrogen atoms of the pyrone and pyridine moieties are not important in the orientation of the drug within the oligonucleotide.

Geometry of intercalation of psoralens in DNA approached by molecular mechanics.

The results of molecular mechanical calculations on intercalation complexes of 3-carbethoxypsoralen, 5-methoxypsoralen, 8-methoxypsoralen, 7-methylpyrido[3,4-c]psoralen (MepyPs) and 7-methylpyrido[4,3-c]psoralen (2N-MePyPs) with the double stranded duodecanucleotide d(CGCGATATCGCG)2 are presented. In the energy-minimized structures, the psoralens are intercalated with their plane orthogonal to the helix axis. Stacking interactions between the furan ring of the psoralen and the adjacent bases are maximized in most derivatives studied, whereas the effect of the various substituents of the psoralen ring is to specifically push part of the molecule towards either the minor or the major groove, preventing a symmetrical intercalation (with respect to the two strands of the DNA). The relative position of the psoralen ring and of the adjacent thymine foreshadows the formation of furan-side monoadducts in 3-CPs, MePyPs and 2N-MePyPs, whereas the formation of a pyrone-side monoadduct appears as geometrically more favourable in 5-MOP and both furan- and pyrone-side monoadducts can be geometrically envisaged in 8-MOP. A good correlation therefore exists between the more or less favourable equilibrium geometries and the experimentally observed photoreactions. The present study is the first attempt to characterize the geometrical parameters as part of a complex set of geometrical, dynamical and excited state parameters governing the overall DNA-psoralen photoreaction.