Parallel double helices of DNA. Conformational analysis of regular helices with the second order symmetry axis

We have performed a conformational analysis of DNA double helices with parallel directed backbone strands connected with the second order symmetry axis being at the same time the helix axis. The calculations were made for homopolymers poly(dA).poly(dA), poly(dC).poly(dC), poly(dG) poly(dG), and poly(dT).poly(dT). All possible variants of hydrogen bonding of base pairs of the same name were studied for each polymer. The maps of backbone chain geometrical existence were constructed. Conformational and helical parameters corresponding to local minima of conformational energy of "parallel" DNA helices, calculated at atom-atom approximation, were determined. The dependence of conformational energy on the base pair and on the hydrogen bond type was analysed. Two major conformational advantageous for "parallel" DNA's do not depend much on the hydrogen-bonded base pair type were indicated. One of them coincided with the conformational region typical for "antiparallel" DNA, in particular for the B-form DNA. Conformational energy of "parallel" DNA depends on the base pair type and for the most part is similar to the conformational energy of "antiparallel" B-DNA.     

Structure-functional organization of the secretin molecule. II. Direct conformation problem

Spatial structure of peptide hormone secretin was investigated by the theoretical conformational method. A solution of the "direct conformational problem" for this hormone indicated that the possible structure of the secretin molecule under polar conditions may be described only by two families of low-energy conformations, possessing relatively conformational valid (Thr7-Leu22) and variable (His1-Phe6 and Leu23-Val27NH2) fragments. One of these families is comprised by five twists of the alpha-helix, while the second isoenergetic family possesses two short segments of the alpha-helix, divided by an irregular structure of the tetrapeptide.     

Parallel DNA double helices. II. Conformational analysis of regular helices having a second order axis of symmetry

Conformational analysis of double helices of DNA with parallel arranged sugar-phosphate chains connected by twofold symmetry has been performed. Homopolymers poly(dA).poly(dA), poly(dC).poly(dC), poly(dG).poly(dG) and poly(dT).poly(dT) were studied. For each of the homopolymers all variants of H-bond pairing were checked. The maps of closing of sugar-phosphate backbone were previously computed. By the optimization of potential energy the dihedral angles and helix parameters of relatively stable conformations of parallel stranded polynucleotides were calculated. The dependence of conformational energy on the nucleic base character and the base pair type were studied. Two main conformational regions for favourable "parallel" helix of polynucleotides were found. The former of these two regions coincide with the region of typical conformational parameters of B-DNA. On an average the conformational energy of "parallel" DNA is close to the energy of canonic "antiparallel" B-DNA.     

Conformational changes of glyceraldehyde-3-phosphate dehydrogenase induced by the binding of NAD. A unified model for positive and negative cooperativity.

The fluorescence of the natural coenzyme, NADH, is used to monitor the environment of the nicotinamide moiety at the active centre of rabbit muscle glyceraldehyde-3-phosphate dehydrogenase (EC 1.2.1.12). Changes of the fluorescence quantum yield and polarization of a small amount of NADH, totally bound by an excess of enzyme, show that at half-saturation of the oligomer with NAD a conformational change is induced which affects the active centre regions of the remaining subunits. This conformational transition is not effected by adenosine diphosphoribose, suggesting that the binding of the nicotinamide moiety of NAD to two subunits is essential for the change of tertiary structure of the remaining subunits that causes the observed changes of the fluorescence properties of the ADH "tracer probe". It is suggested that this conformational transition of the oligomer is responsible for the major decrease of affinity for NAD which occurs at half-saturation, and possibly for the activation by NAD+ of the reductive dephosphorylation reaction catalysed by the enzyme. It is also suggested, by analogy with haemoglobin, that the molecular basis of the negative cooperativity may be the creation of additional intersubunit bonds during the binding of the first two NAD molecules to the tetramer, and a change from a "relaxed" quaternary structure to a "tense" structure at half-saturation.     

Approach to studying the structure-activity relationship of native peptides. I. Structural organization of the opioid peptide Tyr-Glu-Gly-Phe-Met-Arg-Gly-Leu and its artificial analogs]

Theoretical conformational analysis was carried out for the octapeptide Tyr1-Gly2-Gly3-Phe4-Met5-Arg6-Gly7-Leu8. Possible structure of the opioid peptide under physiological conditions may be described by a set of low-energy conformations belonging to 14 different forms of the backbone. The solution of the "reverse conformational problem" for the opioid peptide enables one to predict the modified amino acid sequences (Ala2, D-Ala2, Ala3, D-Ala3, Ala7, D-Ala7, MeMet5, MeArg6-analogues) which may assume one of the low-energy states of the native hormone. The influence of the solute was not taken into account in our calculations.     

The multiple-minima problem in the conformational analysis of polypeptides. III. An Electrostatically Driven Monte Carlo Method: Tests on enkephalin

The three-dimensional conformation of Met-enkephalin, corresponding to the lowest minimum of the empirical potential energy function ECEPP/2 (empirical conformational energy program for peptides), has been determined using a new algorithm, viz. the Electrostatically Driven Monte Carlo Method. This methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. These features are included in the algorithm that produces a Monte Carlo search in the conformational hyperspace of the polypeptide, using electrostatic predictions and a random sampling technique to locate low-energy conformations. In addition, we have incorporated an alternative mechanism that allows the structure to escape from some conformational regions representing metastable local energy minima and even from regions of the conformational space with great stability. In 33 test calculations on Met-enkephalin, starting from arbitrary or completely random conformations, the structure corresponding to the global energy minimum was found inall the cases analyzed, with a relatively small search of the conformational space. Some of these starting conformations wereright orleft-handed -helices, characterized by good electrostatic interactions involving their backbone peptide dipoles; nevertheless, the procedure was able to convert such locally stable structures to the global-minimum conformation.On leave from the National University of San Luis, Faculty of Sciences and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Instituto de Matemática Aplicada, San Luis, Ejército de los Andes 950, 5700 San Luis, Argentina.     

Raman spectroscopy of the "potential sensor" of potential-dependent channels.

Raman spectroscopy (RS) study shows that the "potential sensor" responds to changes in intramembrane potential by conformational changes. The mechanism of regulation of the channel "gate" by the carotenoid potential sensor is discussed.     

UV resonance raman spectra of ligand binding intermediates of sol-gel encapsulated hemoglobin.

We report for the first time specific conformational changes for a homogeneous population of ligand-bound adult deoxy human hemoglobin A (HbA) generated by introducing CO into a sample of deoxy-HbA with the effector, inositol hexaphosphate, encapsulated in a porous sol-gel. The preparation of ligand-bound deoxy-HbA results from the speed of ligand diffusion relative to globin conformational dynamics within the sol-gel (1). The ultraviolet resonance Raman (UVRR) difference spectra obtained reveal that E helix motion is initiated upon ligand binding, as signaled by the appearance of an alpha14beta15 Trp W3 band difference at 1559 cm(-1). The subsequent appearance of Tyr (Y8a and Y9a) and W3 (1549 cm(-1)) UVRR difference bands suggest conformational shifts for the penultimate Tyralpha140 on the F helix, the "switch" region Tyralpha42, and the "hinge" region Trpbeta37. The UVRR results expose a sequence of conformational steps leading up to the ligation-induced T to R quaternary structure transition as opposed to a single, concerted switch. More generally, this report demonstrates that sol-gel encapsulation of proteins can be used to study a sequence of specific conformational events triggered by substrate binding because the traditional limitation of substrate diffusion times is overcome.     

A model for the behaviour of phosphorylase b. The generation of different binding sites via intermediate enzymatic states.

The model given in this paper can be applied to enzymatic systems which have more than two conformational states in equilibrium and which clearly exhibit heterogeneity in the binding of one ligand. The model we propose makes possible quantitative interpretation of our experimental results and of those of many other workers as well. In some cases calorimetric, dialysis and kinetic magnitudes, when plotted against ligand concentration, give multiregional or "stepwise" curves. We suggest that such a behaviour arises because total occupation of one class of binding sites completely moves the enzyme towards a different conformational state in which the affinity for the ligand is greatly increased by the formation of a new class of binding sites. Our calorimetric results for the interaction between some nucleotides and phosphorylase b closely conform to our model.     

Evaluation of conformational changes in the supercoiled DNA of eukaryotic cells by direct nucleoid fluorimetry. I. The comparative information value of various methods of studying structural disorders in the supercoiled DNA in thymocytes and macrophages

A direct fluorimetric method is elaborated for studying conformational distortions in the scDNA of "nucleoids" of macrophages and thymocytes. The method is more advantageous compared to methods of viscosimetry and sedimentation.     

Membrane structure of the human immunodeficiency virus gp41 fusion peptide by molecular dynamics simulation. II. The glycine mutants

In this work, molecular dynamics (MD) simulation of the interaction of three mutants, G3V, G5V and G10V, of the human immunodeficiency virus (HIV) gp41 16-residue fusion peptide (FP) with an explicit palmitoyloleoylphosphatidyl-ethanolamine (POPE) lipid bilayer was performed. The goals of this work are to study the correlation of the fusogenic activity of the FPs with the mode of their interaction with the bilayer and to examine the roles of the many glycine residues in the FP in the fusion process. The results of this work corroborate the main conclusion of our earlier MD work of the WT FP and several mutants with polar substitution. These two studies provide correlation between the mode of insertion and the fusogenic activity of these peptides and support the hypothesis that an oblique insertion of the fusion domain of the viral protein is required for fusogenic activity. Inactive mutants interact with the bilayer by a surface-binding mode. The results of this work, combined with the results of our earlier work, show that, while the secondary structures of the wild-type FP and its mutants do not affect the fusogenic activities, the conformational flexibility appears to be an important factor. The active WT FP and its partially active mutants, G3V and G5V, all have significant conformational transitions at one of the glycine sites. They occur at Gly(5) in FP-wt, at Gly(10) in FP-G5V and at Gly(13) in FP-G3V. Thus, a glycine site in each of these active (or partially active) FPs provides conformational flexibility. On the other hand, the inactive mutants FP-G10V, FP-L9R and FP-V2E do not have any conformational transitions except at either terminus and thus possess no conformational flexibility. Thus, the results of this work support the suggestion that the role of glycine residues in the fusion domain is to provide the necessary conformational flexibility for fusion activity.The glycines also form a "glycine strip" in the FP that locates on one (the less hydrophobic) face of the helix (the "sided helix"). However, whether this "glycine strip" is disrupted or not does not seem to correlate with the retention of fusogenic activities. Finally, although the FLGFL (8-12) motif is absolutely conserved in the HIV fusion domain, a well-structured motif stabilized by hydrogen bonding does not appear to be required for activity. In fact, hydrogen bonding in this motif was found to be missing in FP-G3V and FP-G5V. Both of these mutants are partially active.     

Symmetry-regulated dynamics of multi-enzyme complexes. A model of a pyruvate dehydrogenase complex from Escherichia coli

A dynamic model for quaternary structure of a multienzyme complex is considered. The model is based on the supposition of simultaneously existing similar subunits in a number of different conformational states in the "core" of the multienzyme complex. It is supposed that cyclic conformational transitions of the "core" subunits conserve the symmetry of the entire complex. Such transitions drive the core dynamics as well as the suprastructural multienzyme dynamics. The dynamic model is constructed for the pyruvate dehydrogenase complex from E. coli in a supposition of three different conformers existing in its "core" which correspond to the three steps of the cyclic catalytic process. The model is in accordance with the data from the literature.     

Modeling of the two wave response of ATP receptors to jumps of the agonist concentration in pheochromocytoma cells

A model for the kinetics of conformational transitions of ionotropic ATP receptors in pheochromocytoma cells was elaborated. The contribution of the states of ionotropic receptors (upon the blockage of the "open" channel state) to the kinetics of postsynaptic currents was estimated at mediator concentrations studied. The model enables one to determine the contribution of various conformational states of the receptor, in particular in the "closed" state, to the dynamics of ionic current that is registered upon stimulation of ATP receptors. It is shown that after the cessation of the agonist application, a secondary current wave can arise. The rate constants for conformational transitions of ATP receptors were determined.     

Analysis of the circular dichroism spectrum of proteins using the convex constraint algorithm: a practical guide.

Due to the time scale of circular dichroism (CD) measurements, it is theoretically possible to deconvolute such a spectrum if the pure CD spectra differ significantly from one another. In the last decade several methods have been published aiming at obtaining the conformational weights, or percentages (which are the coefficients for a linear combination) of the so-called typical secondary structural elements making up the three-dimensional structure of proteins. Two methods that can be used to determine the secondary structures of proteins are described here. The first method, called LINCOMB, is a simple algorithm based on a least-squares fit with a set of reference spectra representing the known secondary structures and yielding an estimation of weights attributed to alpha-helix, beta-pleated sheet (mainly antiparallel), beta-turns, unordered form, and aromatic/disulfide (or nonpeptide) contributions of the protein being analyzed. This method requires a "template" or reference curve set, which was obtained from the second method. The second method, "convex constraint analysis," is a general deconvolution method for a CD spectra set of any variety of conformational type. The algorithm, based on a set of three constraints, is able to deconvolute a set of CD curves to its common "pure"-component curves and conformational weights. To analyze a single CD spectrum with this method, the spectrum is appended to the data set used as a reference data set. As a way to determine the reliability of the algorithm and provide a guideline to its usage, some applications are presented.     

Characterization of stable conformational variants of erythrocyte carbonic anhydrases.

Limiting viscosity numbers of bovine and ovine erythrocytes carbonic anhydrase variants were calculated by the objective method of comparing viscosimetric data obtained from low-activity-human erythrocyte carbonic anhydrase and its natural variant. Shifts of mobilities and isoelectric points are shown for all species variants, but variations of limiting viscosity numbers were only detected for human and bovine variants. Results of the study are consistent with the observation that variants arise by deamidation of erythrocyte carbonic anhydrases, and that deamidation is responsible for changes in structure and hydration (i. e. "conformational" modifications). Thus, all the variants so far investigated are stable conformational variants or erythrocyte carbonic anhydrases.     

pH dependence of the circular dichroic bands of phenylmethanesulfonyl-mesentericopeptidase.

The effect of pH on the circular dichroism spectra of phenylmethanesulfonyl-mesentericopeptidase (peptidyl peptide hydrolase, EC 3.4.21) was studied. The ellipticity of the bands below 250 nm, which reflects the backbone conformation of the protein molecule, remains almost unchanged in the pH range 6.2--10.4. However, below pH 6.2 and above pH 10.4 a conformational transition occurs. The pH-dependent changes above 250 nm were also studied. The titration of the CD band at 296 nm reflects the ionization of the "exposed" tyrosines, which phenolic groups are fully accessible to the solvent. An apparent pK of 9.9 is calculated from the titration curve. It is concluded that ionization of the tyrosyl residues with normal pK's is complete before conformational changes in the protein molecule occur.     

Synchronous reversible alterations in enzymatic activity (conformational fluctuations) in actomyosin and creatine kinase preparations.

The phenomenon of synchronism of oscillations of actomyosin and creatine kinase activity in the whole volume of the enzyme preparations was analysed. The synchronous "conformational oscillations" were observed in concentrated gels of actomyosin and in diluted actomyosin and creatine kinase solutions (ATP-creatine N-phosphotransferase, EC 2.7.3.2). The macromolecules of proteins studied may be in two or four conformational states differing enzymatic activity. Large fluctuations become possible in a range of conditions wherein two or four different states, or conformers, are equiprobable. The synchronization of conformational changes of separate macromolecules is maintained with energy derived, for instance, from some oxidative process or dilution of the solution, the process being displayed as conformational oscillations.     

Role of substrate binding forces in exchange-only transport systems: II. Implications for the mechanism of the anion exchanger of red cells

Summary The transition-state theory of exchange-only membrane transport is applied to experimental results in the literature on the anion exchanger of red cells. Two central features of the system are in accord with the theory: (i) forming the transition state in translocation involves a carrier conformational change; (ii) substrate specificity is expressed in transport rates rather than affinities. The expression of specificity is consistent with other evidence for a conformational intermediate (not the transition state) formed in the translocation of all substrates. The theory, in conjunction with concepts derived from the chemistry of macrocyclic ion inclusion complexes, prescribes certain essential properties in the transport site. Separate substites are required for the preferred substrates. Cl^– and HCO ₃ ^– , to account for tight binding in the transition state (K _diss1m). Further, the following mechanism is suggested. A substrate anion initially forms a loose surface complex at one subsite, but in the transition state the subsites converge to form an inclusion complex in which the binding forces are greatly increased through a chelation effect. The conformational change at the substrate site, which is driven by the mounting forces of binding, sets in train a wider conformational change that converts the carrier from an immobile to a mobile form. Though simple, this composite-site mechanism explains many unsual features of the system. It accounts for substrate inhibition, partially noncompetitive inhibition of one substrate by another, and tunneling, which is net transport under conditions where exchange should prevail, according to other models. All three types of behavior result from the formation of a ternary complex in which substrate anions are bound at both subsites. The mechanism also accounts for the enormous range of substrate structures accepted by the system, for the complex inhibition by the organic sulfate NAP-taurine, and for the involvement of several cationic side chains and two different protein domains in the transport site.     

The conformational states of cytochrome oxidase in the mitochondrion.

THE Soret spectrum of "resting" cytochrome oxidase in cytochrome-c depleted mitochondria has been determined. The spectrum obtained is dependent on the rate at which the oxidase is turning over. In the least active preparations, the spectrum is almost pure "oxidized" oxidase. With increasing activity the spectrum is converted to a mixture of "oxidized" and "oxygenated" oxidases. It is concluded that the same conformational differences between the two non-reduced forms that are found in the purified enzyme also occur in these cytochrome-c depleted mitochondria.