Molecular Modeling Studies on Amphotericin B and its Complex with Phospholipid,

Molecular dynamics simulation studies on polyene antifungal antibiotic amphotericin B, its head-to-tail dimeric structure and lipid - amphotericin B complex demonstrate interesting features of the flexibilities within the molecule and define the optimal interactions for the formation of a stable dimeric structure and complex with phospholipid.

Analysis of sequence dependent variations in secondary and tertiary structure of tRNA molecules

The double helical regions of the five tRNA(Phe) and two tRNA(Asp) crystal structures have been analyzed using the local basepair step parameters. The sequence dependent effects in the mini double helices of tRNA are very similar to those observed in the crystal structures of oligonucleotides in the A-form, the purine-pyrimidine and purine-purine steps have small roll angles when compared to the fiber models of A-DNA as well as A-RNA, while the pyrimidine-purine doublet steps have large roll angles. The orientation of the basepairs in the D-stem is unusual and invariant i.e. they are different from the other three stems but are very similar in all the five tRNA(Phe) crystal structures, presumably due to tertiary interaction of the Watson-Crick basepairs with other bases, with all bases being highly conserved. The origin of the differences between the tertiary structures of tRNA(Phe) and tRNA(Asp) from yeast has also been investigated. It is found that even though the angle between the acceptor arm and the D-stem is very similar in the two structures, the angle subtended by the acceptor arm and the anticodon arm is smaller in the tRNA(Phe) structure (by more than 10 degrees). This is due to differences in the orientation of the two mini helices constituting the anticodon arm, which are inclined to each other by approximately 25 degrees in tRNA(Phe) and 16 degrees in tRNA(Asp). In addition, the acceptor arm, the D-stem and the anticodon stem are nearly coplanar in tRNA(Phe), while in tRNA(Asp) the anticodon stem projects out of the plane defined by the acceptor arm and the anticodon stem. These two features together lead to a larger separation between the acceptor and anticodon ends in tRNA(Asp) and indicate that the junction between the D-stem and the anticodon stem is quite variable, with features characteristic of a ball-and-socket type joint and determined for each tRNA molecule by the base sequence at the junction.

Computer modelling studies of ribonuclease A-pyrimidine nucleotide complexes

Different modes of binding of pyrimidine monophosphates 2'-UMP, 3'-UMP, 2'-CMP and 3'-CMP to ribonuclease (RNase) A are studied by energy minimization in torsion angle and subsequently in Cartesian coordinate space. The results are analysed in the light of primary binding sites. The hydrogen bonding pattern brings out roles for amino acids such as Asn44 and Ser123 apart from the well known active site residues viz., His12,Lys41,Thr45 and His119. Amino acid segments 43-45 and 119-121 seem to be guiding the ligand binding by forming a pocket. Many of the active site charged residues display considerable movement upon nucleotide binding.

Computer modeling studies on the binding of 2',5'-linked dinucleoside phosphates to ribonuclease T1-influence of subsite interactions on the substrate specificity

The modes of binding of Gp(2',5')A, Gp(2',5')C, Gp(2',5')G and Gp(2',5')U to RNase T1 have been determined by computer modelling studies. All these dinucleoside phosphates assume extended conformations in the active site leading to better interactions with the enzyme. The 5'-terminal guanine of all these ligands is placed in the primary base binding site of the enzyme in an orientation similar to that of 2'-GMP in the RNase T1-2'-GMP complex. The 2'-terminal purines are placed close to the hydrophobic pocket formed by the residues Gly71, Ser72, Pro73 and Gly74 which occur in a loop region. However, the orientation of the 2'-terminal pyrimidines is different from that of 2'-terminal purines. This perhaps explains the higher binding affinity of the 2',5'-linked guanine dinucleoside phosphates with 2'-terminal purines than those with 2'-terminal pyrimidines. A comparison of the binding of the guanine dinucleoside phosphates with 2',5'- and 3',5'-linkages suggests significant differences in the ribose pucker and hydrogen bonding interactions between the catalytic residues and the bound nucleoside phosphate implying that 2',5'-linked dinucleoside phosphates may not be the ideal ligands to probe the role of the catalytic amino acid residues. A change in the amino acid sequence in the surface loop region formed by the residues Gly71 to Gly74 drastically affects the conformation of the base binding subsite, and this may account for the inactivity of the enzyme with altered sequence i.e., with Pro, Gly and Ser at positions 71 to 73 respectively. These results thus suggest that in addition to recognition and catalytic sites, interactions at the loop regions which constitute the subsite for base binding are also crucial in determining the substrate specificity.

Polypurine-polypyrimidine sequences adopt unwound structure in pBR322 form V DNA as probed by single-hit analysis of HpaII sites.

Structure at the polypurine-polypyrimidine sequences flanking the HpaII sites (CCGG) in pBR322 form V DNA was probed employing single-hit analysis using HpaII restriction endonuclease. Reduced cleavage efficiency of HpaII sites flanked by polypurine-polypyrimidine sequences suggested that under high torsional stress these sequences adopt unwound structures rendering these sites insensitive to restriction enzyme cleavage. In addition to polypurine-polypyrimidine sequences. HpaII sites flanked by alternating purine-pyrimidine sequence, a potential motif of left handed Z-DNA, were also found to be resistant to HpaII cleavage. Results obtained from various studies implicating structure sensitivity of restriction endonucleases and methylases were compiled and a direct correlation was observed between the occurrence of altered sites in a domain and its G/C content in pBR322 form V DNA.

Stereochemistry of 2'-5' linked nucleic acids: crystal and molecular structure of ammonium adenylyl-2',5'-adenosine tetrahydrate: a core fragment of 2'-5' oligo A's produced by interferon induced adenylate synthetase

The preponderance of 3'-5' phosphodiester links in nucleic acids is well known. Albeit less prevalent, the 2'-5' links are specifically utilised in the formation of 'lariat' in group II introns and in the msDNA-RNA junction in myxobacterium. As a sequel to our earlier study on cytidylyl-2',5'-adenosine we have now obtained the crystal structure of adenylyl-2',5'-adenosine (A2'p5'A) at atomic resolution. This dinucleoside monophosphate crystallizes in the orthorhombic space group P2(1)2(1)2(1) with a = 7.956(3) A, b = 12.212(3) A and c = 36.654(3) A. CuK alpha intensity data were collected on a diffractometer. The structure was sloved by direct methods and refined by full matrix least squares methods to R = 10.8%. The 2' terminal adenine is in the commonly observed anti (chi 2 = 161 degrees) conformation and the 5' terminal base has a syn (chi 1 = 55 degrees) conformation more often seen in purine nucleotides. A noteworthy feature of A2'p5'A is the intranucleotide hydrogen bond between N3 and O5' atoms of the 5' adenine base. The two furanose rings in A2'p5'A show different conformations - C2' endo, C3' endo puckering for the 5' and 2' ends respectively. In this structure too there is a stacking of the purine base on the ribose O4' just as in other 2'-5' dinucleoside structures, a feature characteristically seen in the left handed Z DNA. In having syn, anti conformation about the glycosyl bonds, C2' endo, C3' endo mixed sugar puckering and N3-O5' intramolecular hydrogen bond A2'p5'A resembles its 3'-5' analogue and several other 2'-5' dinucleoside monophosphate structures solved so far. Striking similarities between the 2'-5' dinucleoside monophosphate structures suggest that the conformation of the 5'-end nucleoside dictates the conformation of the 2' end nucleoside. Also, the 2'-5' dimers do not favour formation of miniature classical double helical structures like the 3'-5' dimers. It is conceivable, 2-5(A) could be using the stereochemical features of A2'p5'A which accounts for its higher activity.

Calcium ionophore, A23187 and its amino acid complexes: spectroscopic and molecular modeling studies.

The circular dichroism, fluorescence, Nuclear Magnetic Resonance and BLM conductance studies indicate that A23187 forms a stable complex with amino acids at low ionophore concentrations (<10(-4)M). However, A23187 prefers to be in a dimeric structure with no significant binding to amino acids, at concentrations higher than 10(-4)M. It was also observed that at lower concentrations, at which the amino acids bind to the ionophore, the affinity for calcium ions was several orders of magnitude lower than that at higher ionophore concentrations. We have also conducted molecular modeling studies to examine the structure of the A23187 dimer and its amino acid complexes. The results of these modeling studies strongly support our experimental results and validate the formation of a hydrogen bonded and energetically stable A23187 dimer and its amino acid complexes.

Characterization of substrate UpA binding to RNase A--computer modelling and energetics approach.

In the past two decades RNase A has been the focus of diverse investigations in order to understand the nature of substrate binding and to know the mechanism of enzyme action. Although this system is reasonably well characterized from the view point of some of the binding sites, the details of interactions in the second base binding (B2) site is insufficient. Further, the nature of ligand-protein interaction is elucidated generally by studies on RNase A-substrate analog complexes (mainly with the help of X-ray crystallography). Hence, the details of interactions at atomic level arising due to substrates are inferred indirectly. In the present paper, the dinucleotide substrate UpA is fitted into the active site of RNase A Several possible substrate conformations are investigated and the binding modes have been selected based on Contact Criteria. Thus identified RNase A-UpA complexes are energy minimized in coordinate space and are analysed in terms of conformations, energetics and interactions. The best possible ligand conformations for binding to RNase A are identified by experimentally known interactions and by the energetics. Upon binding of UpA to RNase A the changes associated,with protein back bone, Side chains in general and at the binding sites in particular are described. Further, the detailed interactions between UpA and RNase A are characterized in terms of hydrogen bonds and energetics. An extensive study has helped in interpreting the diverse results obtained from a number of experiments and also in evaluating the extent of changes the protein and the substrate undergo in order to maximize their interactions.

Interaction of ecoP15I DNA methyltransferase with oligonucleotides containing the asymmetric sequence 5?-CAGCAG-3

EcoP15I DNA methyltransferase (Mtase) recognizes the asymmeteric sequence CAGCAG and catalyzes the transfer of a methyl group from S-adenosyl-L-methionine to the second adenine residue. We have investigated the DNA binding properties of EcoP15I DNA Mtase using gel mobility shift assays. EcoP15I DNA Mtase binds approximately threefold more tightly to DNA containing its recognition sequence, CAGCAG, than to non-specific sequences in the absence or presence of cofactors. Interestingly, in the presence of ATP the discrimination between specific and non-specific sequences increases significantly. These results suggest for the first time a role for ATP in DNA recognition by type III restriction-modification enzymes. In addition, we have shown that bromodeoxyuridine-containing oligonucleotides form complexes with EcoP15I DNA Mtase that are crosslinked upon irradiation. More importantly, we have shown that the crosslink site is at the site of DNA binding, since it can be suppressed by an excess of unmodified oligonucleotide. EcoP15I DNA Mtase exhibited Michaelis-Menten kinetics with both unmodified and bromodeoxyuridine-substituted DNA, with a higher specificity constant for the latter. Furthermore, gel mobility shift assays showed that proteolyzed EcoP15I DNA Mtase formed a specific complex with DNA, which had similar mobility as the native protein-DNA complex. Taken together these results form the basis fora detailed structure-function analysis of EcoP15I DNA Mtase.

Structure and stability of human telomeric sequence

Telomeric DNA of a variety of vertebrates including humans contains the tandem repeat d(TTAGGG)(n). We have investigated the structural properties of the human telomeric repeat oligonucleotide models d(T(2)AG(3))(4), d(G(3)T(2)A)(3)G(3), and d(G(3)T(2)AG(3)) using CD, gel electrophoresis, and chemical probing techniques. The sequences d(G(3)T(2)A)(3)G(3) and d(T(2)AG(3))(4) assume an antiparallel G quartet structure by intramolecular folding, while the sequence d(G(3)T(2)AG(3)) also adopts an antiparallel G quartet structure but by dimerization of hairpins. In all the above cases, adenines are in the loop. The TTA loops are oriented at the same end of the G tetrad stem in the case of hairpin dimer. Further, the oligonucleotide D(G(3)T(2)AG(3)) forms a higher order structure by the association of two hairpin dimers via stacking of G tetrad planes. Here we show that N-7 of adenine in the hairpin dimer is Hoogsteen hydrogen-bonded. The partial reactivity of loop adenines with DEPC in d(T(2)AG(3))(4) suggests that the intramolecular G quartet structure is highly polymorphic and structures with different loop orientations and topologies are formed in solution. Intra- and interloop hydrogen bonding schemes for the TTA loops are proposed to account for the observed diethyl pyrocarbonate reactivities of adenines. Sodium-induced G quartet structures differ from their potassium-induced counterparts not only in stability but also in loop conformation and interactions. Thus, the overall structure and stability of telomeric sequences are modulated by the cation present, loop sequence, and the number of G tracts, which might be important for the telomere function.

Lipid conformational nomenclature: a general method

We propose a conformational nomenclature for amphiphilic lipid molecules that is general and compatible with the stereospecific numbering scheme, in contrast to earlier methods in which discrepancies with the sn-scheme lead to contradictory assignments of the absolute configuration of the system. The present method can be rationally extended to different classes of lipids, both natural and synthetic. It is simple and provides a convenient framework for conformational studies on widely varying classes of lipids.

Modes of binding of 2'-AMP to RNase T1. A computer modeling study

The modes of binding of adenosine 2'-monophosphate (2'-AMP) to the enzyme ribonuclease (RNase) T1 were determined by computer modelling studies. The phosphate moiety of 2'-AMP binds at the primary phosphate binding site. However, adenine can occupy two distinct sites--(1) The primary base binding site where the guanine of 2'-GMP binds and (2) The subsite close to the N1 subsite for the base on the 3'-side of guanine in a guanyl dinucleotide. The minimum energy conformers corresponding to the two modes of binding of 2'-AMP to RNase T1 were found to be of nearly the same energy implying that in solution 2'-AMP binds to the enzyme in both modes. The conformation of the inhibitor and the predicted hydrogen bonding scheme for the RNase T1-2'-AMP complex in the second binding mode (S) agrees well with the reported x-ray crystallographic study. The existence of the first mode of binding explains the experimental observations that RNase T1 catalyses the hydrolysis of phosphodiester bonds adjacent to adenosine at high enzyme concentrations. A comparison of the interactions of 2'-AMP and 2'-GMP with RNase T1 reveals that Glu58 and Asn98 at the phosphate binding site and Glu46 at the base binding site preferentially stabilise the enzyme-2'-GMP complex.

Molecular dynamics studies on nucleoside 2',3'-cyclic phosphates

2',3'-cyclic nucleotides are intermediates and substrates of Ribonuclease (RNase)-catalysed reactions. The characterization of the equilibrium conformation as well as the flexibility inherent in these molecules helps in understanding the enzymatic action of RNases. The present study explores parameters like phase angle, glycosydic torsion angle and hydrogen bond to find possible interrelationship between them through Molecular Dynamics (MD) simulations on 3'-GMP, 3'-UMP, A>p, G>p, U>p, C>p, GpA>p and UpA>p. Interesting results of the effect of cyclisation and other constraints such as hydrogen bond between certain groups on the equilibrium ribose conformation have emerged from this study.

Functional analysis of conserved motifs in EcoP15I DNA methyltransferase

EcoP15I DNA methyltransferase recognizes the sequence 5'-CAGCAG-3' and transfers a methyl group to N-6 of the second adenine residue in the recognition sequence. All N-6 adenine methyltransferases contain two highly conserved sequences, FxGxG (motif I), postulated to form part of the S-adenosyl-L-methionine binding site and (D/N/S)PP(Y/F) (motif IV) involved in catalysis. We have altered the second glycine residue in motif I to arginine and serine, and substituted tyrosine in motif IV with tryptophan in EcoP15I DNA methyltransferase, using site-directed mutagenesis. The mutant enzymes were overexpressed, purified and characterized by biochemical methods. The mutations in motif I completely abolished AdoMet binding but left target DNA recognition unaltered. Although the mutation in motif IV resulted in loss of enzyme activity, we observed enhanced crosslinking of S-adenosyl-L-methionine and DNA. This implies that DNA and AdoMet binding sites are close to motif IV. Taken together, these results reinforce the importance of motif I in AdoMet binding and motif IV in catalysis. Additionally, limited proteolysis and UV crosslinking experiments with EcoP15I DNA methyltransferase imply that DNA binds in a cleft formed by two domains in the protein. Methylation protection analysis provides evidence for the fact that EcoP15I DNA MTase makes contacts in the major groove of its substrate DNA. Interestingly, hypermethylation of the guanine residue next to the target adenine residue indicates that the protein probably flips out the target adenine residue. (C) 1996 Academic Press Limited