Conformational characteristics of mixed sugar puckered deoxytetranucleoside triphosphate d-GpCpGpC from energy minimization studies

As a continuation of our theoretical studies on nucleic acid subunit systems, in this article we consider the case of the tetranucleoside d-GpCpGpC, the minimally ideal representative unit for analyzing the relative stabilities of different forms of homo- and mixed helical conformation of polynucleotides. The four sugar rings are kept so as to generate B-genus, B+A genus and Z-genus conformations. Twenty five helical conformational states which resulted from judicious mixing of A-, B-, C-, W-, and Z-, states locally are subjected to energy minimization permitting the 19 dihedral angles to vary simultaneously. Conformational states corresponding to regular helical forms and mixed helical forms, when analyzed provide valuable information as to the local conformational flexibility and transitions available to polynucleotides.     

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.     

Conformational Characteristics of Mixed Sugar Puckered Deoxydinucleoside Triphosphate Units d-pCpGp and d-pGpCp from Energy Minimization Studies

Abstract

The deoxydinucleoside triphosphate units d-pCpGp and d-pGpCp were subjected to a rigorous theoretical investigation with a view to describing their distinctive conformational characteristics. For each unit 216 probable three-dimensional forms defined by the backbone-base dihedral angles and sugar pucker modes were considered for conformational energy minimization process and scrutinized with reference to properties, such as base-stacking, hydrogen-bonding, internal flexibility and base sequence-phosphate influence. The P-O bond torsions and the phosphate groups were treated with special attention. The results reveal a number of preferred conformational states other than the known helical forms, such as, A-, B-, C-, Z-, and Watson-Crick conformation. Many interesting one-step (change in only one of the dihedral angles or sugar puckers) conformational transitions which involve just about a kcal/mol of energy came to light. The two base sequences CG and GC were noted to differ strikingly in many of their conformational characteristics.     

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.     

Conformational characteristics of dimeric subunits of RNA from energy minimization studies. Mixed sugar-puckered ApG, ApU, CpG, and CpU.

Following the procedure described in the preceding article, the low energy conformations located for the four dimeric subunits of RNA, ApG, ApU, CpG, and CpU are presented. The A-RNA type and Watson-Crick type helical conformations and a number of different kinds of loop promoting ones were identified as low energy in all the units. The 3E-3E and 3E-2E pucker sequences are found to be more or less equally preferred; the 2E-2E sequence is occasionally preferred, while the 2E-3E is highly prohibited in all the units. A conformation similar to the one observed in the drug-dinucleoside monophosphate complex crystals becomes a low energy case only for the CpG unit. The low energy conformations obtained for the four model units were used to assess the stability of the conformational states of the dinucleotide segments in the four crystal models of the tRNAPhe molecule. Information on the occurrence of the less preferred sugar-pucker sequences in the various loop regions in the tRNAPhe molecule has been obtained. A detailed comparison of the conformational characteristics of DNA and RNA subunits at the dimeric level is presented on the basis of the results.     

Conformational studies of cyclic peptide structures in solution from 1H-Nmr data by distance geometry calculation and restrained energy minimization

The three-dimensional structure of a cyclic bouvardin analogue, cyclo (-Pro-MeTyr-Ala-MeTyr-MeTyr-D-Ala-) has been determined by distance geometry calculation and restrained energy minimization from nmr data. The preparation of the input for the distance geometry calculations, the modification of the amino acid library, and the analysis of the structures were done with the aid of a recently developed software package, GEOM. A great variety of different initial structures were explored to check the uniqueness of the determined solution structure. Calculations with 500 different initial structures and two different strategies led to a uniquely determined backbone conformation with a root mean square deviations value of 0.4 A. The backbone structure consists of two beta-turns, a beta-II turn at Pro1-MeTyr2, and a beta-VI turn at MeTyr4-MeTyr5. The efficiency of the two calculation strategies were compared in order to propose an optimal means for performing distance geometry calculations with cyclic structures.     

Conformational analysis of cyclolinopeptide A, a cyclic nonapeptide: nuclear Overhauser effect and energy minimization studies

The conformation of cyclolinopeptide A [cyclo(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val)], a naturally occurring cyclic nonapeptide has been investigated in dimethylsulfoxide solution by 270 MHz 1H-nmr. A complete assignment of all C alpha H and NH resonances has been accomplished using two-dimensional correlated spectroscopy and nuclear Overhauser effects (NOEs). Analysis of interresidue NOEs and JHNC alpha H values permit construction of a molecular model for the cyclic peptide backbone. The crude model derived from nmr has been used as a starting point for energy minimization, which yields a refined structure largely compatible with nmr observations. The major features of the conformation of cyclolinopeptide A are a Type VI beta-turn centered at Pro(1)-Pro(2), with a cis peptide bond between these residues and a gamma-turn (C7 structure) centered at Ile(6). Two intramolecular hydrogen bonds Val(9) CO--Phe(3)NH (4----1) and Leu(5) CO--Ile(7)NH (3----1) are observed in the low-energy conformation. The limited solvent accessibility observed for the Val(9) and Leu(5) NH groups in the nmr studies are rationalized in terms of steric shielding.     

Conformational analysis of phosphatides: Mapping and minimization of the intramolecular energy

Empirical intramolecular energy calculations were carried out on molecular fragments related to phosphatides in order to find the preferred conformations. The energy was mapped as a function of several pairs of torsional angles in progressively larger molecular fragments, with energy minimization being carried out at each map point with respect to other significant variables. The energy mapping results were used as starting points for energy minimization on diheptanoyl L-α-phosphatidic acid-C, which consisted of the named molecule plus a carbon atom attached to one of the phosphate oxygens. It was found that there are 6 pairs of values for 2 of the torsional angles at the 3-way branch point in the glyceryl group which give sterically acceptable conformations; only 4 of these are compatible with lipid bilayer structure in that they can give a parallel arrangement of the acyl chains. The several acceptable conformations of the phosphate and acyl ester groups within each of these conformational classes are enumerated. The results obtained may be used as a guide for further experimental and theoretical work on phosphatide structures.     

Computation of mixed phosphatidylcholine-cholesterol bilayer structures by energy minimization.

The energetically preferred structures of dimyristoylphosphatidylcholine (DMPC)-cholesterol bilayers were determined at a 1:1 mole ratio. Crystallographic symmetry operations were used to generate planar bilayers of cholesterol and DMPC. Energy minimization was carried out with respect to bond rotations, rigid body motions, and the two-dimensional lattice constants. The lowest energy structures had a hydrogen bond between the cholesterol hydroxyl and the carbonyl oxygen of the sn-2 acyl chain, but the largest contribution to the intermolecular energy was from the nonbonded interactions between the flat alpha surface of cholesterol and the acyl chains of DMPC. Two modes of packing in the bilayer were found; in structure A (the global minimum), unlike molecules are nearest neighbors, whereas in structure B (second lowest energy) like-like intermolecular interactions predominate. Crystallographic close packing of the molecules in the bilayer was achieved, as judged from the molecular areas and the bilayer thickness. These energy-minimized structures are consistent with the available experimental data on mixed bilayers of lecithin and cholesterol, and may be used as starting points for molecular dynamics or other calculations on bilayers.     

Conformational characteristics of the trinucleoside diphosphate d(ApApA) from energy-minimization studies

Energy-minimization studies were carried out on the trinucleoside diphosphate d(ApApA). The potential energy contributions from nonbonded, electrostatic, hydrogen-bonding, and torsional interactions were minimized by treating the 13 relevant dihedral angles as simultaneous variables. For the C(3′)-endo trimer, 14 low-energy conformations are within 10 kcal/mol above the lowest energy found, compared to only 3 in the case of the C(2′)-endo trimer. This result shows the flexible character of the C(3′)-endo unit. The hairpin-type, loop-promoting conformer with (ω′,ω) = (101°, 59°) was found to be the most favored structure at the 3′-terminus of d(ApApA). The predicted U- and L-type bend conformers were found to lie within 5 kcal/mol, compared to the lowest energy B-DNA structure. The A-DNA and Watson-Crick DNA types of helical conformers also lie within very small energy barriers. The phosphate group at the 5′-end of the nucleotide residue has a definite influence on the base of the corresponding nucleotide, keeping it in the normal anti-region, and hence on the base-stacking property. The results are compared with the presently available experimental data, mainly with the tRNA^Phe crystal.     

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K⁺ ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K⁺ channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K⁺ to Na⁺ substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na⁺ locally contracts the SF around Na⁺, releasing ligands that do not participate in Na⁺ coordination, and occluding the permeation pathway. K⁺ selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na⁺ ions due to basic structural differences between liquid water and the “quasi-liquid” SF matrix.     

Conformational characteristics of luliberin. Circular dichroism and fluorescence studies.

By circular dichroism and fluorescence spectroscopy, the conformation of luliberin (luteinizing hormone-releasing hormone) has been investigated under various conditions of pH and solvents. Several structural parameters have been defined which seem predominant for the maintenance of the hormone in some privileged conformation(s). Formation of an intramolecular hydrogen bond between CO (His) and NH (Ser) seems likely when dissolving the hormone in organic solvent such as dioxane. Energy transfer has been demonstrated between Tyr and Trp residues. Calculation of the energy-transfer efficiency at different pH's allowed us to estimate in the range of 10 A the distance which separates these residues. Evidence is also provided for a charge-transfer interaction between protonated histidine and tryptophan. These data suggest that, when luliberin has organized structure (under appropriate surrounding conditions), its conformational pattern would resemble that of beta-turn structure in which a beta bend would exist at the level of the aromatic residues.     

Influence of environment on the antifolate drug trimethoprim: energy minimization studies

Environmental effects on trimethoprim (TMP), an inhibitor of bacterial dihydrofolate reductase (DHFR), were investigated with energy minimizations in vacuo, in the crystal, and in aqueous solution. The conformations, harmonic dynamics, and energetics of the antibacterial drug calculated in these environments were compared with each other and with those of two enzyme-bound drugs. Valence and torsion angles and their energies and overall intra- and intermolecular energies compensated one another in the minimized TMP structures. The conformations of the isolated and aqueous molecules were similar to that of TMP bound to chicken liver DHFR, while the structures from the TMP crystal and from the Escherichia coli DHFR complex were unique. Since neither the small-molecule crystal nor a local minimum of the isolated molecule gave the conformation of TMP bound to the bacterial enzyme, a combination of several experimental and theoretical techniques may be necessary to probe accessible conformations of a molecule.     

Conformational requirements for the polymerization of hemoglobin S: studies of mixed liganded hybrids

Gelation experiments with artificially formed half-liganded hybrid tetramers of hemoglobin S demonstrate that when either the α chains or the β^s chains are fixed in the cyanmet (CNmet) liganded state, gelation occurs upon deoxygenation of the ferrous chains. The minimum concentration of hemoglobin required for gelation is equivalent for both hybrids (α₂^cnmetβ₂^s and α₂β₂^scnmet), is considerably higher than the concentration required to gel deoxy-Hb S (α₂β₂^s), and can be restored to the lower minimum gelling point of α₂β₂^s by reduction of the CNmet chains with dithionite. These results suggest that the most important conformational determinant of the deoxy state for polymerization of Hb S is the quaternary deoxy structure rather than the tertiary structural effect of the ligand state of the α or the β^s chains, and are furthermore consistent with the notion that asymmetric deoxy-CNmet hybrid tetramers assume a conformation which resembles, but is not identical to that of deoxyhemoglobin.The results of gelation experiments with mixtures of hemoglobins S and A in which selected chains of one or both hemoglobins are in the CNmet form support the concept that certain non-S hemoglobins may participate in the sickling process by forming hybrid tetramers with Hb S (such as α₂β^aβ^s). The conformational requirement for participation of these hybrids in polymers also appears to be a quaternary deoxy-like structure.     

Studies on hydrogen bonds. Part V--Hydrogen bonding in energy minimization studies of peptides

An energy term, representing the N-H...O type of hydrogen bond, which is a function of the hydrogen bond length (R) and angle (theta) has been introduced in an energy minimization program, taking into consideration its interpolation with the non-bonded energy for borderline values of R and theta. The details of the mathematical formulation of the derivatives of the hydrogen bond function as applicable to the energy minimization have been given. The minimization technique has been applied to hydrogen bonded two and three linked peptide units (gamma-turns and beta-turns), and having Gly, Ala and Pro side chains. Some of the conformational highlights of the resulting minimum energy conformations are a) the occurrence of the expected 4----1 hydrogen bond in all of the burn-turn tripeptide sequences and b) the presence of an additional 3----1 hydrogen bond in some of the type I and II tripeptides with the hydrogen bonding scheme in such type I beta-turns occurring in a bifurcated form. These and other conformational features have been discussed in the light of experimental evidence and theoretical predictions of other workers.     

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.     

Conformational aspects of polypeptides. XVIII. Conformational studies of oligopeptides derived from L-alanine

The conformations of oligopeptides derived from L -alanine and co-oligomers of L -alanine with γ-methyl-L -glutamate were studied in several solvents via optical rotation and far-ultraviolet spectroscopy. Calculated values for optical rotation based on model compounds were compared with experimental values for the oligomers. In trifluoroacetic and dichloroacetic acids, the oligomers and co-oligomers exhibit rotations in close agreement with predicted values based on model compounds. Thus, in these solvents only nonhelical conformations exist. In trifluoroethanol, the experimental points of molar rotation for the pentamer and larger oligomers no longer follow the predicted values. In addition, the benzyloxycarbonyl and acetyl cononamers show b₀ values of about ?150, which demonstrates the presence of stable helical forms for these peptides. We also examined the molar extinction coefficients of oligopeptides in the 190 mμ region and determined the values for nonhelical peptide groups. The molar extinction coefficients per amide bond for the benzyloxycarbonyl and acetyl cononamers show extensive hypo-chromism, once again indicating the presence of stable helices for these compounds in trifluoroethanol.