Conformational characteristics of the RNA subunit ApA from energy minimization studies

In an attempt to better our understanding of the conformational stabilities in RNAs, an intensive theoraticl study has been carried out on one of its dimeric subunits, ApA, using an improved set of atom-atom interaction energy parameters and an improved version of energy-minimization technique. The C(3′)0endo and the C(2′)-endo sugar ApA units were sperately considered and 38 probable conformations have been analyzed in each case. The total potential energy, comprising nonbonded, electrostatic, and torsional contributions, was minimized by varying all seven relevant dihedral angles simumtaneously. The result reveal that 17 conformations in the case of C(3′)-endo sugar ApA and 7 confomations in the case of C(2′)-endo sugar ApA unit, the lowest energy conformation corresponds to a nonhelical structure and the A-RNA and the Watson-Crick-yype conformations lie at energy levels of about 0.5 and 1.0 Kcal/mo., respectively, above the lowest energy found. For ApA with the lops of different types in the backbone and they all differ in energies by about 3.5 Kcal/mol with refrence to the lowest energy founs. It is noted that the order ofmprefrence of the base stacking is observed in the A-RNA and the Watson-Crick type conformers. The ApA unit with C(2′)-endo sugar is forced to assume phosphodiester conformations with large deviations fom the expected staggered conformations compared to the ApA unit with C(3′)-endo sugar. The result obtained for ApA are discussed with refrence to those previously obtained for the dApdA unit. Te theoretical predictions are compared with the experimental data on the tRNA^Phe crystal, as well as those on fibrous RNAs and RNa subunitlike crystal structures. This study brings out many important aspects of the conformational stability of ApA which have been missed by studies made by others on this system.     

Investigation of fibrous structures. I. Computations for collagen

The structure of collagen is investigated by means of an energy minimization procedure. There is a rather broad trough of potential energy which involves a variety of structures (all collagen II type). It is seen that the unit height and the unit twist can strongly deviate from 2.86 Å, 36° values. Some of the good conformations are presented.     

Conformational characteristics of dApdA,dApdT, dTpdA,and dTpdT from energy minimization studies

The conformational characteristics of the deoxydinucleoside monophosphates with adenine and thymine bases in all possible sequences, namely, dApdA, dApdT, dTpdA, and dTpdT have been studied using an improved set of energy parameters to calculate the total potential energy and an improved set of energy parameters to calculate the total potential energy and an improved version of the minimization technique to minimize the total energy by allowing all seven dihedral angles of the molecular fragment to vary simultaneously. The results reveal that the most preferred conformation in all these units usually corresponds to one of the four helical conformations, namely, the A-DNA, B-DNA, C-DNA, and Watson-Crick DNA models. These helical conformations differ in energies by about 3 kcal/mol with respect to one another. The conformations which could promote a loop or bend in the backbone are, in general, less stable by about 3.5 kcal/mol with respect to the respective lowest-energy helical conformation. The results indicate that there is a definite influence of bases and their actual sequences on the preferred conformations of the deoxydinucleoside monophosphates. The lowest-energy structure, although corresponding to one of the four helical conformations, differ with the type of the deoxydinucleoside monophosphate. Good or reasonable base stacking is noted in dApdA and dTpdA with both C(3′)-endo and C(2′)-endo sugars and in dApdT and dTpdT with only C(3′)-endo sugar. The inversion of the base sequence in deoxydinucleoside monophosphates alters the order of preference of low-energy conformations as well as the base-stacking property of the unit. The paths linking the starting and final states in the (ω′, ω) plane show interesting features with regard to the energy spread, thus providing insight into the path of conformational movement ofthe molecule under slight perturbation. The stabilities of the A and B forms, including the internal energies of the C(3′)-endo ans C(2′)-endo sugar systems, indicate that for dTpdT the B → A transition is less probable. For dApdA, dApdT, and dTpdA this transition is probable in the same order of preference. We propose that the T-A sequence in the polynucleotide chain might serve as the site accessible for B ? A transitions. The theoretical predictions are in good agreement with the experimental observations.     

Stereochemical studies on nucleic acid analogues. I. Conformations of alpha-nucleosides and alpha-nucleotides: interconversion of sugar puckers via O4'-exo.

The preferred conformations of ribo and deoxyribo alpha-nucleosides and alpha-nucleotides, the stereoisomers of the naturally occurring beta-isomers, are worked out by minimizing the conformational energy as a function of all the major parameters including the sugar ring conformations along the pseudorotation path. The results of the studies bring out certain distinct conformational features that are at variance with their beta counterparts. The range of glycosyl conformations are found to be not only quite restricted here but favor predominantly the anti conformation. The syn glycosyl conformation for the entire region of P values are found to be energetically less favorable, with the barrier to anti in equilibrium with syn interconversion being higher especially in alpha-ribonucleosides. The energetically preferred range of pseudorotation phase angles (P) is also considerably restricted and P values corresponding to the C1'-exo range of sugars are highly unfavorable for alpha-nucleosides, in sharp contrast to the broad range of sugar ring conformations favored by beta-isomers. While both trans congruent to 180 degrees and skew congruent to 270 degrees conformations around the C3'-O3' (phi') bond are favored for alpha-3'-nucleotides with deoxyribose sugars, ribose sugars seem to favor only the skew values of phi'. Most interestingly and in sharp contrast to beta-stereoisomers, an energy barrier is encountered at P values corresponding to O4'-endo sugars. This suggests that the possible sugar pucker interconversion between C2'-endo/C3'-exo and C3'-endo/C2'-exo in alpha-anomers could take place only through the O4'-exo region. Likewise the possible path of anti in equilibrium with syn interconversion in alpha-nucleosides is not via high anti, in sharp contrast to beta-nucleosides. These observations should be borne in mind while assigning the sugar ring conformers in alpha-nucleosides and those containing them from nmr investigations. Comparison of the results with beta-anomers seem to suggest on the whole a lack of conformational variability or the restricted nature of alpha-stereoisomers. This could be one of the reasons for its nonselection in the naturally occurring nucleic acids.     

Stereochemical effects of non-ionic methylphosphonates on nucleotide conformations.

The preferred conformations of deoxyribo and ribonucleoside 3'-methylphosphonates are analysed by minimizing the conformational energy as a function of all the major parameters including the sugar ring for both the S- and R-isomers. The results show that neither the substitution nor the nature of the diastereomer affects significantly the preferred conformations compared to the naturally occurring nucleoside 3'-phosphates. The preferred range of C3'-O3' bond torsions or the phase angles of pseudorotation (P) of the sugar are unaffected. The chiral substitution on the phosphate always adopts a conformation distal to the secondary C3' carbon atom in the minimum energy conformational state. Further, it introduces certain restrictions on the preferred range of P-O3' torsions depending on the methylphosphonate configuration. Methylphosphonate, especially the S-isomer, renders the normal gauche- range of P-O3' bond torsions responsible for the stacked helical duplexes to be energetically unfavourable besides introducing a high energy barrier between trans and gauche conformations. Therefore it is suggested that duplexes with S-methylphosphonate may favour extended phosphodiester conformations. These factors explain the observed lower melting temperature as well as the downfield shifts in the 31P signals in duplexes containing the S-isomer.     

Structural studies on the bacterial cell wall peptidoglycan pseudomurein. I. Conformational energy calculations on the glycan strands in C1 conformation and comparison with murein

Conformational energy calculations have been used to explore the conformations which may be realized for the sugar moiety of murein and pseudomurein. For the building blocks of the pseudomurein sugar strands, i.e. for the monosaccharides beta-D-N-acetylglucosamine (NAG) and alpha-L-N-acetyltalosaminuronic acid (NAT), both in C1 ring conformation, as well as for their 1,3 and 1,4 linked disaccharides, the favoured conformations were obtained. The helical parameters of sugar strands of both linkage types, which describe the regular structure of the corresponding polysaccharides, poly-(1,3-NAT-NAG) and poly-(1,4-NAT-NAG), were calculated. Both types of polysaccharides poly-(NAG-NAT) considered in this study favoured extended conformations, which in the case of 1,3 linked polymers showed less gain of length per saccharide unit compared to 1,4 linked poly-(NAG-NAT) residues. For a 1,3 linked sugar moiety of pseudomurein every pair of neighbouring peptides attached to glycan chain pointed in favoured conformations approximately to opposite sides of the strands, whereas in a 1,4 linked poly-(NAG-NAT) the peptides protruded approximately to the same side of the glycan moiety. A comparison between pseudomurein and murein revealed that the sugar moieties of both peptidoglycans have similar features in respect to their overall structure, i.e. both favoured more or less extended structures. In contrast to these data the shapes of the resulting polysaccharide moieties were remarkably different. In poly-(1,3-NAG-NAT) the glycan chains possessed a zig-zag-like arrangement, whereas for glycan chains of the murein type relatively flat structures were preferred. These remaining differences in the conformational arrangement between both peptidoglycans depend strongly on the C1 chair conformation of NAT. It is, therefore, attractive to speculate about an hypothetical pseudomurein sugar chain configuration comprising beta-L-N-acetyltalosaminuronic acid in its 1C conformation.     

The crystal and molecular structure of 2'-deoxy-2'-fluoroinosine monohydrate.

The structure of the hydrate of 2'-deoxy-2'-fluoroinosine has been determined by single-crystal x-ray diffraction. The nucleoside crystallizes in space group P2(1)2(1)2(1) with unit cell dimensions, a = 33.291, b = 10. 871, c = 6.897A. There are two nucleosides and two water molecules in the asymmetric unit. The structure was solved by direct methods and refined to a residual R = 0.095. The two independent nucleosides in the asymmetric unit show different conformations about the glycosidic bond, while other structural details are similar. The base orientation to the sugar is syn in molecule A, whereas anti in molecule B. The exocyclic C(4')-C(5') bond conformation defined with respect to C(3')-C(4')-C(5')-O(5') is gauche+ in both molecules A and B. The sugar ring pucker defined by the pseudorotation phase angle P is a twisted conformation in both, C(3')-endo-C(4')-exo with P = 29 degrees in molecule A and C(4')-exo-C(3')-endo with P = 41 degrees in molecule B. It is shown by comparison with x-ray results of other 2'-fluoronucleosides and unmodified nucleosides including inosines that, in addition to a strong preference of the C(3')-endo type pucker, twisted conformations involving C(4')-exo puckering may be one of characteristic features of 2'-fluoronucleosides.     

Determination of the conformation of d(GGAAATTTCC)2 in solution by use of 1H NMR and restrained molecular dynamics

The conformation of the putative bent DNA d(GGAAATTTCC)2 in solution was studied by use of 1H NMR and restrained molecular dynamics. Most of the resonances were assigned sequentially. A total of 182 interproton distance restraints were determined from two-dimensional nuclear Overhauser effect spectra with short mixing times. Torsion angle restraints for each sugar moiety were determined by qualitative analysis of a two-dimensional correlated spectrum. Restrained molecular dynamics was carried out with the interproton distances and torsion angles incorporated into the total energy function of the system in the form of effective potential terms. As initial conformations for restrained molecular dynamics, classical A-DNA and B-DNA were adopted. The root mean square deviation (rmsd) between these two conformations is 5.5 A. The conformations obtained by use of restrained molecular dynamics are very similar to each other, the rmsd being 0.8 A. On the other hand, the conformations obtained by use of molecular dynamics without experimental restraints or restrained energy minimization depended heavily on the initial conformations, and convergence to a similar conformation was not attained. The conformation obtained by use of restrained molecular dynamics exhibits a few remarkable features. The second G residue takes on the BII conformation [Fratini, A. V., Kopka, M. L., Drew, H. R., & Dickerson, R. E. (1982) J. Biol. Chem. 257, 14686-14707] rather than the standard BI conformation. There is discontinuity of the sugar puckering between the eighth T and ninth C. The minor groove of the oligo(dA) tract is rather compressed. As a result, d(GGAAATTTCC)2 is bent.     

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by ±10° in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3⁷) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150–325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CX_nG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5′ end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10⁶ conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar–phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques. © 1993 John Wiley & Sons, Inc.     

Theoretical studies on the conformations of higher gangliosides

The possible conformations of higher gangliosides (GD3, GT1a. GT1b, GQ1b) have been determined by computing their potential energy using semi-empirical potential functions. The favoured conformation of the disialic acid fragment in these gangliosides is independent of its position (internal or terminal). The favoured conformations of these gangliosides have also been correlated to their biological activity. The results suggest that tetanus toxin and sendai virus may have a large binding site which can accommodate at least four sugar residues.     

Conformational studies of (2''-5'') polynucleotides: theoretical computations of energy, base morphology, helical structure, and duplex formation.

A detailed theoretical analysis has been carried out to probe the conformational characteristics of (2'-5') polynucleotide chains. Semi-empirical energy calculations are used to estimate the preferred torsional combinations of the monomeric repeating unit. The resulting morphology of adjacent bases and the tendency to form regular single-stranded structures are determined by standard computational procedures. The torsional preferences are in agreement with available nmr measurements on model compounds. The tendencies to adopt base stacked and intercalative geometries are markedly depressed compared to those in (3'-5') chains. Very limited families of regular monomerically repeating single-stranded (2'-5') helices are found. Base stacking, however, can be enhanced (but helix formation is at the same time depressed) in mixed puckered chains. Constrained (2'-5') duplex structures have been constructed from a search of all intervening glycosyl and sugar conformations that form geometrically feasible phosphodiester linkages. Both A- and B-type base stacking are found to generate non-standard backbone torsions and mixed glycosyl/sugar combinations. The 2'- and 5'-residues are locked in totally different arrangements and are thereby prevented from generating long helical structures.     

Conformational energy calculations of enzyme-substrate complexes of lysozyme. I. Energy minimization of monosaccharide and oligosaccharide inhibitors and substrates of lysozyme

Conformational energy calculations were performed on monosaccharide and oligosaccharide inhibitors and substrates of lysozyme to examine the preferred conformations of these molecules. A grid-search method was used to locate all of the low-energy conformational regions for N-acetyl-β-D -glycosamine (NAG), and energy minimization was then carried out in each of these regions. Three stable positions for the N-acetyl group have ben located, in two of which the plane of the amide unit is normal to the mean plane of the pyranosyl ring. Nine local energy minima were located for the —CH₂OH group. The positions of the two vicinal cis —OH groups are determined predominantly by interactions with either the —CH₂OH or the N-acetyl group. The most stable conformations of β-N-acetylmuramic acid (NAM) were determined from the study of the low-energy conformations of NAG. In the two stable orientations for the D -lactic acid side chain, the O—C—C′ plane (C′ being the carbon atom of the terminal carboxyl group) was found to be normal to the mean plane of the pyranosyl ring. The low-energy positions for the COOH group of NAM are determined mainly by interactions with neighboring groups. The conformational preferences of the α-anomers of NAG and NAM were also explored. The calculated conformation of the N-acetyl group for α-NAG was quite close to that determined by X-ray analysis. Two of the three lowest energy conformations of α-NAM are similar to the corresponding conformations of the β-anomer. A third low-energy structure, which has a hydrogen bond from the NH of the N-acetyl group to the C?O of the lactic acid group, corresponds very closely to the X-ray structure of this molecule. The preferred conformations of the disaccharides NAG–NAG, NAM–NAG and NAG–NAM were also investigated. Two preferred orientations of the reducing pyranosyl ring relative to the nonreducing ring were found for all of these disaccharides, both of which are close to the extended conformation. In one of these conformations, a hydrogen bond can form between the OH group attached to C₃ of the reducing sugar and the ring oxygen of the preceding residue. Each conformation can be stabilized further by a hydrogen bond between the CH₂OH (donor) of residue i + 1 and the C?O of residue i (acceptor). The interactions that determine conformations for all oligosaccharides containing both NAG and NAM are shown to be exclusively intraresidue and nearest neighbor interactions, so that it is possible to predict all stable conformations of oligosaccharides containing NAG and NAM in any sequence.     

Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A

The four Trp dipoles in the gramicidin A (gA) channel modulate channel conductance, and their side chain conformations should therefore be important, but the energies of different conformations are unknown. A conformational search for the right-handed helix based on molecular mechanics in vacuo yielded 46 conformations within 20 kcal/mol of the lowest energy conformation. The two lowest energy conformations correspond to the solid-state and solution-state NMR conformations, suggesting that interactions within the peptide determine the conformation. For representative conformations, the electrostatic potential of the Trp side chains on the channel axis was computed. A novel application of the image-series method of. Biophys. J. 9:1160-1170) was introduced to simulate the polarization of bulk water by the Trp side chains. For the experimentally observed structures, the CHARm toph19 potential energy (PE) of a cation in the channel center is -1.65 kcal/mol without images. With images, the PE is -1.9 kcal/mol, demonstrating that the images further enhance the direct dipole effect. Nonstandard conformations yielded less favorable PEs by 0.4-1.1 kcal/mol.     

The Study of Possible A and B Conformations of Alternating DNA Using a New Program for Conformational Analysis of Duplexes (CONAN)

Abstract

A new program, CONAN has been designed for CONformational ANalysis of oligonucleotide duplexes with natural and modified bases. It allows to model both regular DNA fragments with different types of symmetry and irregular ones including bends, junctions, mismatched pairs and base lesions. Computations and minimization of the energy are performed in a space of internal structural variables chosen to build start structure easier and conveniently analyze the results obtained. These internal structural variables determine mutual base-base and base-sugar arrangement and sugar puckering. The analytical closure procedure is applied both to sugar rings and to backbone fragments between adjacent sugars. For more effective energy minimization, analytical gradient is calculated. The CONAN was applied to the search for low-energy conformations of poly(dA-dT)·poly(dA-dT) and poly(dG-dC)·poly(dG-dC) duplexes. Extended regions of low-energy A and B conformations are revealed and characterized. These regions contain structures with different relative values of helical twist, τ, for pur-pyr and pyr-pur steps, namely, conformations with τ(pur-pyr)>τ(pyr-pur) and with τ(pur-pyr)<τ(pyr-pur). Two types of sugar puckering were found for B-form low-energy conformations, the first type with all C2′-endo sugar residues and the second one—;with C2′-endo purines and O1′-endo pyrimidines. The calculated conformations are compared with X-ray diffraction data for crystals and fibers and NMR data for solution.     

Conformational behaviour of the architectural units of peptides and proteins. Assessment of current understanding by ab initio quantum mechanical methods and refinement of the dipeptide energy surface.

The conformational energy surfaces of analogues of the dipeptide unit of polypeptides and proteins are calculated by ab initio methods using extended basis sets.The calculations are not particularly sensitive to the choice of (extended) basis set.The calculations are shown to support a particular empirical method parameterized with respect to crystal data. Non-hydrogen bonded conformations agree to within 3 kcal mol^?1, even for conformations in which quite considerable degrees of atomic overlap occur.Hydrogen bonded conformations, are, however, in less satisfactory agreement and it is the ab initio calculations which appear to be at fault.A simple correction is applied to the ab initio energy for hydrogen bonded conformations, and with the use of the empirical energy surface a full quantum mechanical conformational energy map is interpolated for the alanyl dipeptide.The effect of flexibility in the peptide backbone is taken into account, and supports recent empirical findings that distortions in valence angles must be considered in calculations of the conformational behaviour of peptides.     

Reversible scaling of dihedral angle barriers during molecular dynamics to improve structure prediction of cyclic peptides.

Peptide cyclization or chemical cross-linking has frequently been used to restrict the conformational freedom of a peptide, for example, to enhance its capacity for selective binding to a target receptor molecule. Structure prediction of cyclic peptides is important to evaluate possible conformations prior to synthesis. Because of the conformational constraints imposed by cyclization low energy conformations of cyclic peptides can be separated by large energy barriers. In order to improve the conformational search properties of molecular dynamics (MD) simulations a potential scaling method has been designed. The approach consists of several consecutive MD simulations with a specific lowering of dihedral energy barriers and reduced nonbonded interactions between atoms separated by three atoms followed by gradually scaling the potential until the original barriers are reached. Application to four cyclic penta- and hexa-peptide test cases and a protein loop of known structure indicates that the potential scaling method is more efficient and faster in locating low energy conformations than standard MD simulations. Combined with a generalized Born implicit solvation model the low energy cyclic peptide conformations and the loop structure are in good agreement with experiment. Applications in the presence of explicit water molecules during the simulations showed also improved convergence to structures close to experiment compared with regular MD.     

Molecular and crystal structure of Sp-thymidin-3'-yl 4-thiothymidin-5'-yl methylphosphonate.

The molecular structure of one diastereomer of the dinucleoside methylphosphonate Tp(Me)sT (1) has been determined by X-ray diffraction methods. The crystal asymmetric unit contains one molecule of 1 and one methanol in an orthorhombic unit cell of dimensions a = 13.241(4), b = 13.844(3), c = 14.944(7) A, space group P2(1)2(1)2(1). Both pyrimidine bases in 1 are oriented anti relative to the 2'-deoxyribose rings, and the sugar conformations are 2E and 2(3)T in the 4-thiothymidine and thymidine moieties, respectively. The deoxyribose-phosphonate backbone has an extended conformation with the bases completely unstacked and almost parallel. The absolute configuration at the phosphorus center in 1 is Sp.     

Stability of polypeptide conformational states. II. Folding of a polypeptide chain by the scanning simulation method, and calculation of the free energy of the statistical coil

The scanning simulation method suggested by Meirovitch is extended to a study of the stability of decaglycine at 100 and 300 K. The model is based on the potential energy function ECEPP (Empirical Conformational Energy Program for Peptides) with rigid geometry and without solvent. The free energy of the statistical coil, which is defined over the whole phase space excluding the region of the right-handed α-helix, is calculated. At 100 K, the molecule is found to be unstable in the statistical coil region, and the method generates (i.e., “folds”) conformations that are left-handed or right-handed α-helices with very high preference. Their free energy is found to be comparable with that obtained by another method developed in our previous paper (paper I) [H. Meirovitch, M. Vásquez, and H. A. Scheraga, (1987) Biopolymers 26 , 651–671]. At 300 K the statistical coil becomes the most stable state; sample conformations of the coil are generated efficiently with the scanning method and the free energy is calculated. It appears that both the scanning method and the method of paper I can be used to carry out a complete analysis of the stability of a polypeptide based on free energy considerations.     

Discrimination of the native from misfolded protein models with an energy function including implicit solvation.

An essential requirement for theoretical protein structure prediction is an energy function that can discriminate the native from non-native protein conformations. To date most of the energy functions used for this purpose have been extracted from a statistical analysis of the protein structure database, without explicit reference to the physical interactions responsible for protein stability. The use of the statistical functions has been supported by the widespread belief that they are superior for such discrimination to physics-based energy functions. An effective energy function which combined the CHARMM vacuum potential with a Gaussian model for the solvation free energy is tested for its ability to discriminate the native structure of a protein from misfolded conformations; the results are compared with those obtained with the vacuum CHARMM potential. The test is performed on several sets of misfolded structures prepared by others, including sets of about 650 good decoys for six proteins, as well as on misfolded structures of chymotrypsin inhibitor 2. The vacuum CHARMM potential is successful in most cases when energy minimized conformations are considered, but fails when applied to structures relaxed by molecular dynamics. With the effective energy function the native state is always more stable than grossly misfolded conformations both in energy minimized and molecular dynamics-relaxed structures. The present results suggest that molecular mechanics (physics-based) energy functions, complemented by a simple model for the solvation free energy, should be tested for use in the inverse folding problem, and supports their use in studies of the effective energy surface of proteins in solution. Moreover, the study suggests that the belief in the superiority of statistical functions for these purposes may be ill founded.     

Stereochemistry of nucleic acids and polynucleotides. VI. Minimum energy conformations of dimethyl phosphate

As part of a study on the conformation of polynucleotides and nucleic acids the preferred conformations of the model conpound dimethyl phosphate are worked out using potential energy functions. In calculating the total potential energy associated with the conformation, nonbonded, torsional, and electrostatic terms have been considered. The variation of the total conformational energy is represented as a function of two torsion angles ? and ψ which are the rotations about the two phosphoester bonds. The most stable conformations are found to be the gauche–gauche conformations about these bonds. The conformations observed for phosphodiesters in the solid state and in the proposed structures of polynucleotides and nucleic acids cluster around the minimum. Also, regions of minimum energy correspond well with the typical allowed regions of a representative dinucleotide.