Optimized monopole expansions for the representation of the electrostatic properties of polypeptides and proteins

A procedure for refining a coiled coil protein structure using x-ray fiber diffraction and modeling

We describe a combined use of experimental and simulation techniques to configure side chains in a coiled coil structure. As already demonstrated in a previous work, x-ray diffraction patterns from hard alpha-keratin fibers in the 5.15 A meridian zone reflect the global configuration of the chi(1) dihedral angle of the coiled coil side chains. Molecular simulations, such as energy minimization and molecular dynamics, and rotameric representation in the PDB, are used here on a heterodimeric coiled coil to investigate the dihedral angle distribution along the sequence. Different procedures have been used to build the structure, the quality assessment was based on the agreement between the simulated diffraction patterns and the experimental ones in the fingerprint region of coiled coils (5.15 A). The best one for building a realistic coiled coil structure consists of placing the side chains using molecular dynamics (MD) simulations, followed by side chain positioning using SMD or SCWRL procedures. The side chains and the backbone are equilibrated during the MD until they reach an equilibrium state for the t/g(+) ratio. Positioning the side chains on the resulting backbone, using the above procedures, gives rise to a well-defined 5.15 A meridian reflection.

Unraveling proteins: a molecular mechanics study

An internal coordinate molecular mechanics study of unfolding peptide chains by external stretching has been carried out to predict the type of force spectra that may be expected from single-molecule manipulation experiments currently being prepared. Rather than modeling the stretching of a given protein, we have looked at the behavior of simple secondary structure elements (alpha-helix, beta-ribbon, and interacting alpha-helices) to estimate the magnitude of the forces involved in their unfolding or separation and the dependence of these forces on the way pulling is carried out as well as on the length of the structural elements. The results point to a hierarchy of forces covering a surprisingly large range and to important orientational effects in the response to external stress.

Prediction of the positioning of the seven transmembrane alpha-helices of bacteriorhodopsin. A molecular simulation study

We have applied a search strategy for determining the optimal packing of protein secondary structure elements to the rotational positioning of the seven transmembrane helices of bacteriorhodopsin. The search is based on the assumption that the relative orientations of the helices within the bundle are conditioned principally by inter-helix side-chain interactions and that the extra-helical parts of the protein have only a minor influence on the bundle conformation. Our approach performs conformational energy optimization using a predetermined set of side-chain rotamers and appropriate methods for sampling the conformational space of peptide fragments with fixed backbone geometries. The final solution obtained for bacteriorhodopsin places each of the seven helices to a precision of a few degrees in rotation around the helical axis and to a few tenths of an ångström in translation along the helical axis with respect to the best experimental structure obtained by electron diffraction, except for helix D, where our results support the suggestion that this helix should be displaced along its axis toward its N terminus. The perspectives of such an approach for the determination of the structures of other transmembrane helical bundles are discussed.

Packing and recognition of protein structural elements: a new approach applied to the 4-helix bundle of myohemerythrin

We present a novel search strategy for determining the optimal packing of protein secondary structure elements. The approach is based on conformational energy optimization using a predetermined set of side chain rotamers and appropriate methods for sampling the conformational space of peptide fragments having fixed backbone geometries. An application to the 4-helix bundle of myohemerythrin is presented. It is shown that the conformations of the amino acid side chains are largely determined at the level of helix pairs and that superposition of these results can be used to construct the full bundle. The final solution obtained, taking into account restrictions due to the lateral amphiphilicity of the helices, differs from the native structure by only a 20 degrees rotation of a single helix.

Static and dynamic conformational properties of AT sequences in B-DNA

A theoretical study of the optimal conformations of nucleic acid oligomers containing tracts of AT base pairs is presented. The oligomers are studied in isolation and complexed with netropsin, a minor groove binding ligand. The flexibility of the oligomers and of their complexes is calculated by adiabatic mapping with respect to the total winding angle. The results of this study show that in uncomplexed oligomers the dinucleotide junctions AA, AT and TA have very different structural parameters and different responses to winding stress. The TA junction is clearly the most flexible and is the principal site for accommodating the imposed overwinding. Complexation by netropsin leads to two important effects: firstly, the three junctions adopt more uniform structures, the largest changes again being observed for TA, secondly, the differences in flexibility as a function of sequence are strongly attenuated.

How can the aromatic side-chains modulate the conductance of the gramicidin channel? A new approach using non-coded amino acids

In order to elucidate the role of the aromatic side-chains in the mechanism of transduction of monovalent cations through the channel of linear gramicidin, two series of analogues containing non-coded aromatic amino acids were synthesized. In the first series, the four tryptophans were replaced by either four L-3-(8-quinolyl)alanyl or four L-3-(4-quinolyl)alanyl residues and single channel conductance measurements showed that these substitutions led to a strong lowering of the channel conductance, which is attributed to a modification of the orientation of the aromatic side-chains due to an increase of their hydrophobicity. In the second series, the analogues contained both tryptophyl and naphthylalanyl residues in various amounts and positions. The single channel conductance data indicated that the conductance was mainly governed by the number of polar residues (Trp) and not by their positions. The conformational consequences of these results are discussed together with their influence on the energy profile of the gramicidin channel.

A new approach to the rapid determination of protein side chain conformations

Two efficient algorithms have been developed which allow amino acid side chain conformations to be optimized rapidly for a given peptide backbone conformation. Both these approaches are based on the assumption that each side chain can be represented by a small number of rotameric states. These states have been obtained by a dynamic cluster analysis of a large data base of known crystallographic structures. Successful applications of these algorithms to the prediction of known protein conformations are presented.

Induction of DNA bending by bifunctional intercalating agents of the 7H-pyridocarbazole family

The structures and binding energetics of selected complexes formed between the deoxynucleotides d(CpGpGpCpG).d(CpGpCpCpG), d(CpGpApTpCpG)2, d(GpCpGpCpCpG).d(CpGpGpCpGpC), and d(CpGpCpCpCpG)2 with the DNA bifunctional intercalating agent ditercalinium and three of its rigid linking chain derivatives have been investigated theoretically by means of a molecular mechanics approach that takes into account nucleic acid flexibility, ligand flexibility and solvent dielectric effects (R. Lavery, in: Unusual DNA structures, eds S. Harvey and R. Wells (Pergamon, New York, 1988) p. 189; R. Lavery, in: DNA bending and curvature, eds W.K. Olson et al. (Adenine Press, New York, 1988) p. 191). The piperidinium chains of the bis-intercalating ligands are always located in the major groove of DNA. For the energy-minimized complexes the ligand proceeds to bind following preferentially the 5'-pyrimidine-purine-3' alternating sequence, thus dictating the number of internal exclusion sites. The complexes with three exclusion sites will present (i) a bending of the structure towards the major groove, and (ii) a non-ideal distribution of unwinding angles; complexes with less than three exclusion sites will remain essentially linear. The absence of a bend does not preclude other types of local deformations of the base-pairs such as inclination, buckle and tip. The proposed structures of the d(CpGpApTpCpG)2 complexes are in agreement with NMR structural results. The possible relevance of these findings to a previously proposed mode of interaction for ditercalinium-like DNA ligands is discussed.

The intrinsic molecular potential of glyceryl monooleate layers and its effect on the conformation and orientation of an inserted molecule: example of gramicidin A

It is shown by explicit calculation that the distribution of the atomic charges in the constituent molecules of a lipid monolayer or bilayer of glyceryl monooleate creates an intrinsic potential difference between the head region and the hydrocarbon region which tends to repel positive charges towards the exterior and attract negative charges to the interior. The analogies and differences between a bilayer and a monolayer are analyzed. The possible consequences of the intrinsic potential gradient in a lipid layer on the preferred orientation and conformation of a polar neutral molecule are illustrated on the case of a gramicidin A monomer.

Intercalative binding of ditercalinium to d(CpGpCpG)2: a theoretical study

The structure of the complex formed between ditercalinium, 2,2'-[4,4'-bipiperidine-1,1'-bis-(ethane-1,2-diyl)]bis(10-me thoxy-7H- pyrido[4,3-c]carbazolium) tetramethane sulfonate (NSC 366241), and the self-complementary tetranucleotide duplex d(CpGpCpG)2 has been investigated by means of a novel theoretical approach for modelling the conformational flexibility of nucleic acids. The methodology used is the JUMNA procedure, a molecular mechanics systematics capable of evaluating the internal energy and the interaction energy of a complex formed from a large number of fragments. In the best energy-minimized structures, the piperidinium chains of ditercalinium are located in the major groove of the right-handed oligonucleotide. Calculations show a distortion of the base-paired d(CpGpCpG)2 minihelix consisting of lateral dislocation of one base pair with respect to another along an axis parallel to the long axis; strong propeller twist and tilt of the end base pairs; a collective motion of all base pairs with respect to the helical axis towards the drug; and an overwinding at the exclusion site. The proposed structure of the complex is in good agreement with reported proton NMR data, supporting the feasibility of such model.

A possible model for the inner wall of the acetylcholine receptor channel

A structural model of the inner wall of the acetylcholine receptor (AChR) channel is developed using assumptions derived from the results of the recent labelling experiments of the MII helices by noncompetitive blockers. The assumption of steric blocking of the channel by chlorpromazine (CPZ) in the neighbourhood of the labelled serines imposes the MII helices to be in contact at this level and allows the calculation of their minimal interaxial distance. The assumption that CPZ diffuses to this position through the upper crowded part of the channel imposes that the helices are more distant in this region and permits the determination of a tilt of about 7 degrees with respect to the central axis. Electrostatic potentials are used to demonstrate the effect of the charged residues at the exit of the pore. A discussion is given on the possible aptitude of MI to satisfy the contacts necessary with the MII/s at the different heights of the model.

Conformation and pairing properties of the N-terminal fragments of trichorzianine and alamethicin: a theoretical study

Energy optimizations are carried out on the N-terminal fragment of trichorzianine in comparison to that of alamethicin. The results indicate that the helical character of the (Ac...Pro13) sequence of trichorzianine (TA IIIc) is essentially alpha with a bend in the helix axis in the end proline region, a structure comparable to the optimal alpha-helical structure of the corresponding segment (Ac...Pro14) of alamethicin AI. However, two weak n----n + 3 interactions coexist in trichorzianine with the alpha-helical n----n + 4 hydrogen bonds. The possible role of the glutamine side-chains in pairing such segments together is considered.

Theoretical study of the sequence selectivity of isolexins, isohelical DNA groove binding ligands. Proposal for the GC minor groove specific compounds

A theoretical study is presented of complex formation between DNA fragments of different base sequences and isolexins, "isohelical base reading polymers", formed of heteroaromatic pentagonal rings joined by appropriate linkers. Extensive computations are performed for the isolexin composed of the furan-pyrrole-furan sequence. They involve charged ligands with propioamidinium groups at both ends as well as neutral molecules with terminal methyl, carbonyl and amino groups. Two different groups (C = O and NH) are used as linkers between the base reading moieties. The role of these elements on the binding preference of the ligands has been examined. The results show that the mere possibility of formation of hydrogen bonds between a ligand and the nucleic acid bases is not sufficient to ensure its binding specificity which is determined largely by the interplay of electrostatic factors. Thus the dicationic isolexins uniformly prefer AT sequences. For the neutral isolexins the nature of the groups forming the linkers is a major factor in defining the specificity, although these groups do not participate directly in the interaction with DNA. The C = O linkers favour binding to AT sequence while the N-H linkers permit preferential binding to the GAG sequence. Finally, for the first time in theoretical computations, a ligand is proposed which should bind preferentially to the minor groove of GC sequences: this ligand is a neutral isolexin composed of three furan rings linked by two N-H groups. This ligand is considered as an improvable prototype. Altogether the results presented open the path for the designing of minor groove ligands specific for any desirable DNA base sequence.

Active site model of cytochrome P-450 LM2

Based on (i) a detailed analysis of the physicochemical properties of selected benzphetamine derived substrates and (ii) the identification of Tyr-380 as active site residue trans to thiolate theoretical studies (computer aided molecular design) revealed a model of the substrate binding site of cytochrome P-450 LM2. The results indicate that substrates with a butterfly-like bulky conformation exhibit the highest intrinsic activity. Those substrates which preferably exist in an extended conformation are sterically hindered to intensively interact with the binding site which is demonstrated by computer graphics.

A theoretical study of the sequence specificity in binding of lexitropsins to B-DNA

A theoretical study is presented on the binding to B-DNA of a series of lexitropsins, these ligands being netropsin derivatives in which one or both of the pyrrole rings have been replaced by imidazoles. The best complexes have been located by energy minimisation taking into account nucleic acid flexibility, ligand flexibility, explicit, mobile counterions and solvent dielectric effects. Calculations have been performed for two homopolymeric DNA receptor sequences, AT base sequence, which only decreases in the imidazole derivatives. These results emphasize the decisive role of the molecular electrostatic potential of the nucleic acid in determining the sequence selectivity of these ligands, as opposed to the postulated role of adenine C2 - pyrrole beta hydrogen contacts.

Sequence selectivity, a test of the nature of the covalent adduct formed between benzo[a]pyrene and DNA

A theoretical study is presented of the energetic and structural properties of covalent adducts of benzo[a]pyrene and a DNA fragment. Energy optimisation is performed with the use of minimiser with constraints and an advanced semiempirical energy formula. Three types of adducts are studied: an external complex with the benzopyrene located in the DNA minor groove and two types of intercalative complexes with the carcinogen situated on the 3' side and 5' side of the covalently bound guanine. For each of the adducts the effects of DNA base sequence are examined. It is shown that the results for the intercalative complex with the carcinogen situated on the 5' side of the modified guanine correlate with the experimentally determined sequence preference.

Theoretical studies on the interaction of proteins and nucleic acid. II. The binding of alpha-helix to B-DNA

Interactions between B-DNA and homopolymeric alpha-helices of glycine, alanine, serine, asparagine and aspartic acid have been studied theoretically. The complexation energy has been minimised taking into account the interactions between DNA and the polypeptides as well as the internal energy of the alpha-helix and the interaction energy of counterions with the complex. The results obtained indicate the important role of strong hydrogen bonds between the peptide side chains and nucleic acid phosphate groups, these bonds being much stronger than specific interactions with the base-pairs. The formation of these structural bonds depends on the size of the alpha-helix, which in turn determines whether bridging across the major groove is possible. The steric role of the methyl group of thymine in orienting the peptide helix and the role of DNA screening cations in complex stabilization are also significant.

A general approach to the optimization of the conformation of ring molecules with an application to valinomycin

A general and efficient methodology is presented which allows molecules containing one or many rings of any size to be manipulated within energy minimization procedures. Variables describing the conformation of the molecules concerned are limited to dihedral and ring valence angles and the ring closure conditions are treated as equality constraints. An application is made to the ion transporter valinomycin and its complexes with K+ and Na+ which illustrates the possibilities of the approach and leads to results which allow a better understanding of the conformational mechanics of this important ionophore.

A theoretical investigation of the sequence specificity in the binding of the antitumor drug anthramycin to DNA

A theoretical study is presented concerning DNA-anthramycin adducts. By explicit energy minimisations using a semi-empirical energy formula and an advanced algorithm the structural properties and the energetics of this system are analysed. The results obtained demonstrate that the formation of a covalently bound adduct in which anthramycin is attached to the N2 site of a guanine within a DNA fragment is accompanied by a considerable change in the nucleic acid conformation as confirmed by recent experimental evidence. With the use of the "SIR" methodology for treating DNA flexibility the general features of this change are characterised. The sequence specificity of anthramycin binding is investigated and the important role of sequence dependent nucleic acid flexibility brought to light. This theoretical treatment thus provides new elements for the interpretation of the origins of ligand binding specificities.

Binding of non-intercalating antibiotics to B-DNA: a theoretical study taking into account nucleic acid flexibility

A detailed theoretical study has been made for five antibiotics which all bind selectively to AT sequences in the minor groove of B-DNA: SN-18071, NSC-101327, distamycin-2, distamycin-3 and netropsin. The optimal complexes were found for systems in which the flexibility of DNA, as well as that of the antibiotics, was taken into account. Explicit, mobile counterions and a dielectric function modelling aqueous solution were also included. The binding geometries of the most strongly interacting antibiotics, distamycin-3 and netropsin, are compared in considerable detail and it is shown that notable differences exist between them. The results for netropsin are also discussed in the light of recent disagreements concerning its exact binding location within DNA.

The flexibility of the nucleic acids: (II). The calculation of internal energy and applications to mononucleotide repeat DNA

Results concerning the flexibility of mononucleotide repeat DNA are presented using a novel methodology, denoted "SIR", to describe continuous changes in the structure of the nucleic acid. This methodology, combined with internal energy calculations and analytical energy gradients allows us to determine optimal conformations of poly(dG).poly(dC) and poly (dA).poly(dT) in both the A and B forms, taking into account the influence of the solvent medium and explicit counterions. Subsequently, several different types of distorsion of these optimal structures are investigated. It is shown that excellent correlation with experimental results concerning coupled changes in structural variables is obtained and several new correlations are also detected.

Theoretical studies on the interaction of proteins and nucleic acids. I. The binding of beta-pleated sheets to A- and B-DNA

A theoretical investigation of the interaction between a beta-ribbon consisting of two glycine hexapeptides and DNA in its A and B conformations is presented. A refined semi-empirical energy formula and a sophisticated energy minimization technique are used to optimize the complex, taking into account the DNA-beta-ribbon interaction, the full flexibility of the oligopeptide chains and of the positions of the DNA screening counterions. A considerable flexibility of the beta-ribbon is demonstrated, which allows the polypeptide fragment to interact comfortably with both forms of DNA considered and with different base-pair sequences. The results are discussed in connection with the general problem of DNA-protein recognition.

The gramicidin A channel: energetics and structural characteristics of the progression of a sodium ion in the presence of water

The distribution of water molecules in the Gramicidin A (GA) channel is determined by theoretical computations, and the role of this water on the energetics of the system upon progression of a sodium cation through the channel is investigated. In the absence of the ion, water molecules form a chain along the channel, hydrogen bonded to one another and to the L carbonyl oxygens, while others stay at the entrances of the channel, hydrogen-bonded to the free carbonyl oxygens of the L-Tryptophan residues. According to the definition adopted for the "inside" and the "outside" of the channel, it is found to contain at most 7 or 9 water molecules. When a hydrated sodium cation approaches and enters the channel, the structural properties corresponding to the minimized total energy of the system GA-water-Na+ indicate a reorganization, but not a destruction, of the chain of water molecules. The "energy profile" for the system GA-Na+-(22 waters) is analyzed in terms of its components and in comparison to the corresponding intrinsic profile computed earlier in vacuo. It appears that the presence of water does not unduely modify the pathway or the qualitative features of the energetics of the cation passage, except at the entrance, where the partial and progressive dehydration of the cation plays an important role. The presence and characteristics of the minimum found earlier at 10.5 A from the center are conserved.

The effect of spermine binding on the reactivity of DNA towards carcinogenic alkylating agents

The effect of spermine binding on the electrostatic potential of DNA is evaluated. The calculations are performed for the essential reactive sites, atoms N7 and O6 of guanine, N3 and N7 of adenine, of the nucleic acid and for its surface envelope. An important weakening of the potential is found affecting all the important reactive sites in both grooves and spreading moreover along the polynucleotide chain far away from the site of binding of the ligand. These results are discussed in connection with the experimentally observed inhibitory effect of spermine binding on DNA methylation by carcinogenic agents.