The molecular electrostatic potential and steric accessibility of A-DNA.

The molecular electrostatic potential and steric accessibility of A-DNA are computed for base sequences (dG.dC)n and (dA.dT)n. An interpretation of the results in terms of the structure of A-DNA is provided and differences with respect to other forms of DNA, namely B-DNA and Z-DNA, are discussed.     

The electrostatic molecular potential of tRNAPhe. IV. The potentials and steric accessibilities of sites associated with the bases.

The sites of the 76 nucleic acid bases of tRNAPhe potentially reactive towards electrophiles are studied by calculations on the associated molecular electrostatic potentials and the static steric accessibilities. Each of these sites is treated in its environment within the macromolecule. The influence of various schemes of screening by countercations of the backbone phosphates on the electrostatic potentials is investigated. The possible significance of the potentials and accessibilities in connection with observed chemical reactivities is discussed.     

The molecular electrostatic potential and steric accessibility of poly (dI.dC). Comparison with poly (dG.dC)

The influence of the amino group of guanine on the molecular electrostatic potential and the accessibility to reactive sites of B-DNA is investigated by comparing the two model double helices poly (dI.dC) and poly (dG.dC). The calculations clarify the “disruptive” role of the guanine amino group on nucleic acid-polypeptide interactions.     

The molecular electrostatic potential and steric accessibility of poly (dA-dT). poly (dA-dT) in various conformations: B-DNA, D-DNA and 'alternating-B' DNA.

The influence of conformational changes on the molecular electrostatic potential and the steric accessibility of the double stranded polynucleotide poly (dA-dT). poly (dA-dT) are investigated by calculating these properties for three different conformations : B-DNA, D-DNA and alternating-B DNA.     

Glycocalyx electrostatic potential profile analysis: ion, pH, steric, and charge effects

The Poisson-Boltzmann equation is modified to consider charge ionogenicity, steric exclusion, and charge distribution in order to describe the perimembranous electrostatic potential profile in a manner consistent with the known morphology and biochemical composition of the cell's glycocalyx. Exact numerical and approximate analytical solutions are given for various charge distributions and for an extended form of the Donnan potential model. The interrelated effects of ionic conditions, bulk pH, ion binding, local dielectric, steric volume exclusion, and charge distribution on the local potential, pH, and charge density within the glycocalyx are examined. Local charge-induced, potential-mediated pH reductions cause glycocalyx charge neutralization. Under certain conditions, local potentials may be insensitive to ionic strength or may decrease in spite of increasing charge density. The volume exclusion of the glycocalyx reduces the local ion concentration, thereby increasing the local potential. With neutral lipid membranes, the Donnan and surface potential agree if the glycocalyx charge distribution is both uniform and several times thicker than the Debye length (approximately 20 A in thickness under physiological conditions). Model limitations in terms of application to microdomains or protein endo- and ectodomains are discussed.     

Mapping the electrostatic potential within the ribosomal exit tunnel

Electrostatic potentials influence interactions among proteins and nucleic acids, the orientation of dipoles and quadrupoles, and the distribution of mobile charges. Consequently, electrostatic potentials can modulate macromolecular folding and conformational stability, as well as rates of catalysis and substrate binding. The ribosomal exit tunnel, along with its resident nascent peptide, is no less susceptible to these consequences. Yet, the electrostatics inside the tunnel have never been measured. Here we map both the electrostatic potential and accessibilities along the length of the tunnel and determine the electrostatic consequences of introducing a charged amino acid into the nascent peptide. To do this we developed novel probes and strategies. Our findings provide new insights regarding the dielectric of the tunnel and the dynamics of its local electric fields.     

The properties of steric gate mutants reveal different constraints within the active sites of Y-family and A-family DNA polymerases

Y-family (lesion-bypass) DNA polymerases show the same overall structural features seen in other members of the polymerase superfamily, yet their active sites are more open, with fewer contacts to the DNA and nucleotide substrates. This raises the question of whether analogous active-site side chains play equivalent roles in the bypass polymerases and their classical DNA polymerase counterparts. In Klenow fragment, an A-family DNA polymerase, the steric gate side chain (Glu710) not only prevents ribonucleotide incorporation but also plays an important role in discrimination against purine-pyrimidine mispairs. In this work we show that the steric gate (Phe12) of the Y-family polymerase Dbh plays a very minor role in fidelity, despite its analogous role in sugar selection. Using ribonucleotide discrimination to report on the positioning of a mispaired dNTP, we found that the pyrimidine of a Pu-dPyTP nascent mispair occupies a similar position to that of a correctly paired dNTP in the Dbh active site, whereas in Klenow fragment the mispaired dNTP sits higher in the active site pocket. If purine-pyrimidine mispairs adopt the expected wobble geometry, the difference between the two polymerases can be attributed to the binding of the templating base, with the looser binding site of Dbh permitting a variety of template conformations with only minimal adjustment at the incoming dNTP. In Klenow fragment the templating base is more rigidly held, so that changes in base pair geometry would affect the dNTP position, allowing the Glu710 side chain to serve as a sensor of nascent mispairs.     

The electrostatic potential of B-DNA

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson-Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence-specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate-phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.     

Predicting potential miRNA target sites within gene promoters

Synthetic small duplex RNAs that are complementary to gene promoters can activate or inhibit target gene expression. The potency and robustness of gene modulation by these RNAs suggests that natural mechanisms may exist to facilitate recognition of sequences within gene promoters by endogenous small RNAs. Here, we describe computational methods for identifying potential miRNA target sites within gene promoters. These methods will facilitate investigations of whether miRNAs interact with sequences outside of 3′-untranslated regions and suggest new targets for the design of synthetic modulators of gene expression.     

An unusual functional group interaction and its potential to reproduce steric and electrostatic features of the transition states of peptidolysis

The donor-acceptor interaction between a tertiary amine and an aldehyde, first observed among a select class of alkaloids, was deliberately established in a peptidomimetic (1a-c) to mimic features of the two principal transition states of peptide hydrolysis. Compounds 1a-c show preferential adoption in methanol and water of a 'folded' conformation displaying the interaction. Proportions of the folded form in MeOH range from 45% to 70% and can reach 84% in buffer. Significantly, three tendencies for the folded/unfolded equilibrium are observed: increasing solubility and polarity of the medium and decreasing temperature results in a higher extent of folding. In the absence of any parameter set available for this weak bond, no modeling studies were conducted to aid in the design of 1a-c. The successful straightforward synthesis of 1 and its folding and inhibition results with HIV-1 peptidase using FRET technology encourage studies to further pre-organize candidate molecules and to screen the structure space by modeling and parallel combinatorial chemistry.     

Blue fluorescent antibodies as reporters of steric accessibility in virus conjugates

Nonenveloped viruses provide the chemist with large, preassembled polyvalent protein scaffolds for modification. These structures are typically porous to small molecules but not to large ones. The solution-phase structures and reactivities of such assemblies may be substantially different than indicated by X-ray crystal structures. Here, the attachment of organic compounds to either the inside or outside surface of the cowpea mosaic virus (CPMV) coat protein was verified with an indicating antibody-antigen interaction. Antibody binding was subsequently blocked by the installation of poly(ethylene glycol) chains. These results typify the type of site-specific control that is available with CPMV and related virus building blocks.     

Unusual class of Alu sequences containing a potential Z-DNA segment.

A potential Z-DNA sequence, (dA-dC)9, has been found to replace the customary A-rich region in the middle of an Alu family member in the African green monkey genome. This Alu, bounded by imperfect direct repeats, also contains an unusual 3' end and may be a member of a large subfamily of such sequences.     

The electrostatic molecular potential of yeast tRNAPhe. (I). The potential due to the phosphate backbone.

We present a preliminary theoretical study of the electrostatic potential surrounding yeast tRNAPhe by computing the component of this potential due to all the 76 phosphate groups of the molecule. The general features of the results obtained are discussed and deductions concerning the binding of Mg2+ ions to the molecule in relation to the potential are presented.     

The electrostatic potential of Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.     

Chemical detection of Z-DNA within the maize Adh1 promoter

Z-DNA is a left-handed helix which can form within tracts of alternating purines and pyrimidines. Tracts of potential Z-DNA identified by sequence inspection are often noted within regulatory portions of genes, but evidence that these tracts of sequence actually exist as Z-DNA is very limited, and not available for any plant gene. In this study, the chemical probes osmium tetroxide, diethylpyrocarbonate and hydroxylamine were used to show that a tract of alternating purines and pyrimidines in the Adh1 promoter (from -311 to -325) actually assumes a Z-DNA conformation under superhelical stress in vitro.     

Supercoil-induced Z-DNA formation within 5'-end of chicken myb proto-oncogene

We have analyzed the recently sequenced and characterized 2.9 kb fragment derived from the 5'-end of chicken myb proto-oncogene with respect to structural perturbations induced by DNA supercoiling. Within the first intron a 50 bp sequence stretch was localized, starting approximately 450 nucleotides downstream from putative ATG initiation codon, which forms a non-B-DNA structure. Fine mapping with structural probes revealed the three adjacent regions with imperfect purine-pyrimidine alternation creating together relatively long Z-forming tract, parts of which may undergo a B-Z DNA transition at different superhelical densities.