Sequence-dependent molecular conformation of polynucleotides: right and left-handed helices

Based upon a stereochemical guideline, two topologically distinct types of helicalduplexes have been deduced for a polynucleotide duplex with alternating purine pyrimidine sequence (PAPP): (a) right-handed uniform (RU) helix and (b) left-handed zig-zag (LZ) helix. Both structures have trinucleoside diphosphate as the basic unit wherein the purine pyrimidine fragment has a different conformation from the pyrimidine-purine fragment. Thus, RU and LZ helices represent two different classes of sequence-dependent molecular conformations for PAPP. The conformationalf eatures of an RU helix of PAPP in B-form and three LZ-helices for B-, D- and Z-forms are discussed.

Sterochemical studies on cyclic peptides-VIII. Conformational analysis of hydrogen bonded cyclohexaglycyl molecule with a centre of inversion symmetry

A study on the conformational aspects of cyclo-hexaglycyl having inversion symmetry has been made. The cyclic backbone has been assumed to have two internal 4→1 types of NH... O hydrogen bonds. This molecule has been found to take up two types of conformations designated asA* andB* having nearly the same energy values. The theoretical conformations have been compared with the conformations of cyclohexaglycyl hemihydrate observed in the crystal structure. Two molecules with an approximate inversion symmetry are close to the conformation of the typeB* and two other molecules with exact inversino symmetry correspond nearly to the typesB* andA*. comparison with the theoretically possible conformations of cyclohexaglycyl molecule with 2-fold symmetry has been made. The preference of inversion symmetry and preferred ranges ofψ for glycyl molecules is discussed.

Intermolecular potentials in the dimer, the excimers, and the dimer ions of ethylene

Calculations are reported on the interaction energies in the dimer, the excimers, and the dimer ions of ethylene. The various a- and u-electron terms for different conformations of the dimeric species are determined by using the exchange perturbation method. The results predict that the singlet excimer and the dimer cation are stable primarily because of the large magnitude of the exciton-resonance and charge-resonance terms, respectively, while the neutral dimer, the triplet excimer, and the dimer anion are weakly stable. The variations of the various energy terms with conformations suggest that these dimeric species cannot have identifical structure.

Application of ab initio molecular orbital calculations to the structural moieties of carbohydrates. 3

Ab initio RHF/4-31G level molecular orbital calculations have been carried out on dimethoxymethane as a model compound for the acetal moiety in methyl pyranosides. The calculations are consistent with the predictions of the anomeric effect and the exo-anomeric effect. They reproduce very successfully the differences in molecular geometry observed by x-ray and neutron diffraction of single crystals of the methyl cy-D- and methyl 0-D-pyranosides. Calculations carried out at the 6-3 1G* level for methanediol confirm the earlier calculations at the 4-31G level, with smaller energy differences between the four staggered conformations.

A single point mutation disrupts the capsid assembly in Sesbania Mosaic Virus resulting in a stable isolated dimer

Protein-protein interactions play a Crucial role in Virus assembly and stability. With the view of disrupting capsid assembly and capturing smaller oligomers, interfacial residue mutations were carried Out in the coat protein gene of Sesbania Mosaic Virus, a T=3 ss (+) RNA plant virus. A single point mutation of a Trp 170 present at the five-fold interface of the virus to a charged residue (Glu or Lys) arrested assembly of virus like particles and resulted in stable Soluble dimers of the capsid Protein. The X-ray crystal structure of one of the isolated dimer mutants - rCP Delta N65W170K was determined to a resolution of 2.65 angstrom. Detailed analysis of the dimeric mutant protein structure revealed that a number of Structural changes take place, especially in the loop and interfacial regions during the course of assembly. The isolated chiller was ``more relaxed'' than the dimer found in the T=3 or T=1 capsids. The isolated dimer does not bind Ca2+ ion and consequently four C-terminal residues are disordered. The FG loop, which interacts with RNA in the Virus, has different conformations in the isolated dimer and the intact Virus Suggesting its flexible nature and the conformational changes that accompany assembly. The isolated choler mutant was much less stable when compared to the assembled capsids, suggesting the importance of inter-subunit interactions and Ca2+ mediated interactions in the stability of the capsids. With this study, SeMV becomes the first icosahedral virus for which X-ray crystal Structures of T=3, T=1 capsids as well as a smaller oligomer of the capsid protein have been determined.

Conformationally Restricted Nucleocyclitols: a Study into their Conformational Preferences and Supramolecular Architecture in the Solid State

With the intent of probing the feasibility of employing annulation as a tactic to engender axial rich conformations in nucleoside analogues, two adenine-derived, ``conformationally restricted'' nucleocylitols, 9 and 10, have been conceptualized as representatives of a hitherto unexplored class of nucleic acid base-cyclitol hybrids. A general synthetic strategy, with an inherent scope for diversification, allowed rapid functionalization of indane and tetralin to furnish 9 and 10 respectively in fair yield. Single-crystal X-ray diffraction analysis revealed that the two nucleocyclitols under study, though homologous, present completely dissimilar modes of molecular packing, marked, in particular, by the nature of involvement of the adenynyl NH2 group in the supramolecular assembly. In addition, the crystal structures of 9 and 10 also exhibit two different conformations of the functionalized cyclohexane ring. Thus, while the six-membered carbocycle in cyclopenta-annulated 9 exists in the expected chair (C) conformation that in cyclohexaannulated 10, which crystallizes as a dihydrate, shows an unusual twist-boat (TB) conformation. From a close analysis of the (HNMR)-H-1 spectroscopic data recorded for 9 and 10 in CD3OD, it was possible to put forth a putative explanation for the uncanny conformational preferences of crystalline 9 and 10.

A Case Study of the Conformation of Poly(alpha-aminoisobutyric acid): alpha- or 310-Helix

The relative stabilities of a- and Blo-helical structures for polymers of a-aminoisobutyric acid (Aib) have been worked out, using the classical potential energy functions. To make a comparative study, we have used Buckingham "6-exp" and Kitaigorodsky's potential functions. Conformational analysis of the dipeptide segment with Aib residue indicates the necessity for nonplanar distortion of the peptide unit, which is a common feature in the observed crystal structures with Aib residues. In the range of Aw -10 to +loo studied, a-helical conformations are preferred in the region -3" < Aw < +loo, and Blo-helical conformations are preferred in the region -3" > Aw > -10'. Minimum energy conformations for right-handed structures are found in the +ue region of Aw and correspondingly for left-handed structures in the -ue region of Aw. For Aw - 6", a-helical structures have four- or near fourfold symmetry with h - 1.5 A. Such a helix with n = 4 and h = 1.5 A is termed an a'-helix. This structure is found to be consistent with the electron diffraction data of Malcolm3 and energetically more favorable than the standard 310-helix.

Conformational analysis of cyclo (l-cystine

Conformational analysis of cyclo(L-cystine) shows that the diketopiperazine ring has to exist only in the boat form. With this geometry, the molecule can adopt two distinct forms differing mainly in the chirality of the disulphide bridge. In both the P- and M-models, corresponding to dihedral angles of nearly + 90° and —90° respectively about the S-S bond, the molecule displays an approximate two-fold symmetry. According to our semi-empirical energy calculations, the minimum energy of the M-model is —9.2 kcal/mol, only 0.3 kcal/mol lower than that of the P-model. Because the difference between the two minima is so small, neither form is clearly superior to the other. However, the number of low energy conformations of the M-model in the allowed conformational space is significantly larger than that of the P-model by a ratio of 3 to 1, and therefore the former is likely to be thermodynamically favoured.

Conformational Choices for the Stereochemically Constrained gamma-Amino Acid Residue Gabapentin: Theoretical Studies and Correlation With Experimental Results

Gabapentin (1-aminomethylcyclohexaneacetic acid, Gpn) is an achiral, conformationally constrained gamma amino acid residue. A survey of available crystal structures of Gpn peptides reveals that the torsion angles about the C-gamma-C-beta (theta(1)) and C-beta-C-alpha(theta(2)) bonds are overwhelmingly limited to gauche, gauche (g(+)g(+)/g(-)g(-)) conformations. The Gpn residue forms C-7 and C-9 hydrogen bonds in which the donor and acceptor atoms come from the flanking peptide units. In combination with alpha amino acid residues alpha gamma and gamma alpha segments can adopt C-12 hydrogen bonded structures. The conformational choices available to the Gpn residue have been probed using energy calculations, adopting a grid search strategy. Ramachandran phi-psi maps have been constructed for fixed values of theta(1) and theta(2), corresponding to the gauche and trans conformations. The sterically allowed and energetically favorable regions of conformational space have been defined and experimental observations compared. C-7 and C-9 hydrogen bonded conformational families have been identified using a grid search approach in which theta(1) and theta(2) values are varied over a range of +/- 10 degrees about ideal values at 1 degrees intervals. The theoretical analysis together with experimental observations for 59 Gpn residues from 35 crystal structures permits definition of the limited range of conformational possibilities at this gamma amino acid residue. .

Conformational studies on cyclic dipeptides

An analysis of 11 crystal structures of cyclic dipeptides so far reported in the literature is made, with main reference to the internal parameters of these molecules. Preferred conformations of the side chains of cyclic dipeptides with different α-amino acid residues have been studied by classical energy calculations. The possible conformations of the DKP ring are also studied. The significance of the non-bonded interaction in deciding the pathway for conformational change has also been investigated. The agreement between theoretical results and experimental observations is quite good, both with respect to the conformation of these molecules as well as the enthalpy difference as estimated from n.m.r. studies between different conformers.

Computational analyses of the catalytic and heparin-binding sites and their interactions with glycosaminoglycans in glycoside hydrolase family 79 endo-d-glucuronidase (heparanase)

Mammalian heparanase is an endo-β-glucuronidase associated with cell invasion in cancer metastasis, angiogenesis and inflammation. Heparanase cleaves heparan sulfate proteoglycans in the extracellular matrix and basement membrane, releasing heparin/heparan sulfate oligosaccharides of appreciable size. This in turn causes the release of growth factors, which accelerate tumor growth and metastasis. Heparanase has two glycosaminoglycan-binding domains; however, no three-dimensional structure information is available for human heparanase that can provide insights into how the two domains interact to degrade heparin fragments. We have constructed a new homology model of heparanase that takes into account the most recent structural and bioinformatics data available. Heparin analogs and glycosaminoglycan mimetics were computationally docked into the active site with energetically stable ring conformations and their interaction energies were compared. The resulting docked structures were used to propose a model for substrates and conformer selectivity based on the dimensions of the active site. The docking of substrates and inhibitors indicates the existence of a large binding site extending at least two saccharide units beyond the cleavage site (toward the nonreducing end) and at least three saccharides toward the reducing end (toward heparin-binding site 2). The docking of substrates suggests that heparanase recognizes the N-sulfated and O-sulfated glucosamines at subsite +1 and glucuronic acid at the cleavage site, whereas in the absence of 6-O-sulfation in glucosamine, glucuronic acid is docked at subsite +2. These findings will help us to focus on the rational design of heparanase-inhibiting molecules for anticancer drug development by targeting the two heparin/heparan sulfate recognition domains.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

X-ray studies on crystalline complexes involving amino acids and peptides. XXIII. Variability in ionization state, conformation and molecular aggregation in the complexes of succinic acid with DL- and L-lysine

Crystalline complexes of succinic acid with DL- and L-lysine have been prepared and analysed by X-ray diffraction. DL-Lysine complex: C6HIsN202 + 1 2- 1 ~C4H404 .~C4H604, Mr -- 264"2, PI, a = 5"506 (4), =8.070(2), c=14.089(2) A,, a=92.02(1), /3= 100"69 (3), y = 95"85 (3) ~>, Z = 2, Dx = 1"44 g cm -3, R = 0.059 for 2546 observed reflections. Form I of the e-lysine complex: C6HIsN20-, ~ .C4H504, Mr = 264.2, P1, a = 5" 125 (2), b = 8"087 (1), c = 8"689 (1) A,, a = 112.06 (1), /3 = 99.08 (2), y = 93"77(2) °, Z--l, D,,,=1"34(3), Dx=l"34gcm 3 R = 0.033 for 1475 observed reflections. Form II of + I 2- the e-lysine complex: C6H15N202 .,iC4H404 .- 1 I ") 4C4H604.4(C4HsO4""H'"CaH404)" , Mr = 264"2, P1, a = 10.143 (4), b = 10.256 (2), c = 12"916 (3) A,, a = 105.00 (2),/3 = 99-09 (3), y = 92"78 (3)::, Z = 4, Dm= 1"37(4), D,.= 1.38gcm 3, R=0.067 for 2809 observed reflections. The succinic acid molecules in the structures exhibit a variety of ionization states. Two of the lysine conformations found in the complexes have been observed for the first time in crystals containing lysine. Form II of the L-lysine complex is highly pseudosymmetric. In all the complexes, unlike molecules aggregate into separate alternating layers. The basic element of aggregation in the lysine layer in the complexes is an S2-type head-to-tail sequence. This element combines in different ways in the three structures. The basic element of aggre gation in the succinic acid layer in the complexes is a hydrogen-bonded ribbon. The ribbons are interconnected indirectly through amino groups in the lysine layer.