A Computer-Aided Radiopharmaceutical Drug Design Study Using Ab Initio and Molecular Mechanics Methods

This is an investigation of technetium ligands and their complexes with [TcO]³⁺ using ab initio population analysis and molecular mechanics conformational searching methods. Calculated atomic electronic populations on the technetium atom in complexes with a number of ligands gauge the degree of covalent bonding between technetium and these ligands. Here a reduction in the positive charge on the [TcO]³⁺ moiety by complexation with a given ligand is correlated with covalent bonding. Our ab initio results suggest that ligands with more sulphur atoms have better covalent bonding to technetium than do other ligands. A conformational analysis of the uncomplexed ligands indicates that conformational reorganization before complexation correlates inversely with stable complex formation. This conformational analysis shows that ligands with ethylene carbonyl bridges have low energy conformations closer to the final complexation geometries than do ligands with ethylene, propylene or propylene carbonyl bridges. The presence of these low energy conformations facilitates a faster complexation of the ethylene carbonyl [TcO]³⁺ moiety. This result produces a kinetic explaination why ethylene carbonyl bridged ligands form stable complexes while many other ligands do not [1]. The conclusion is that kinetic and thermodynamic considerations play a role in stable complex formation between these ligands and technetium.     

Conformation of histidine model peptides. I. Conformational energy calculations for L-alanyl-L-histidine diketopiperazine

Semiempirical conformational energy calculations were carried out for the cyclic dipeptide L -alanyl-L -histidine diketopiperazine. The results indicate that electrostatic effects are probably significant in determining the conformation assumed by this molecule. When the imidazole group is in its uncharged state the most stable conformations of the molecule are those with the imidazole ring folded over the diketopiperazine ring (χ¹ = 60°). Upon protonation of the imidazole group the folded conformation may be destabilized relative to conformations characterized by χ¹ positions near 180°.     

Molecular Parameterisation and the Theoretical Calculation of Electrode Potentials

The difference in reduction potentials between ortho and para-benzoquinones has been calculated. The employs gas phase ab initio and semi-empirical computations in combination with free energy perturbation theory applied to gas and solution phase Monte Carlo simulations. The effects on calculated results of altering solute electrostatic parameterisation in solution phase simulations is examined. Atom centred charges derived from the molecular electrostatic potentials, MEPs, from optimised ab initio wavefunctions and charges generated by consideration of hydrogen bonded complexes are considered. Parameterisation of hydroxyl torsions in hydroquinone molecules is treated in a physically realistic manner. The coupled torsional system of the ortho-hydrobenzoquinone molecule is described by a potential energy surface calculated using gas phase AM1 semi-empirical computations rather than the simple torsional energy functions frequently employed in such calculations. Calculated differences in electrode potentials show that the electrostatic interactions of quinone and hydroquinone molecules in aqueous solution are not well described by atom centred charges derived from ab initio calculated MEPs. Moreover, results in good agreement with the experimental reduction potential difference can be obtained by employing high level ab initio calculations and solution phase electrostatic parameters developed by consideration of hydrogen bonded complexes.     

Calculations on the low energy conformers of N-acetyl-D-alanyl-D-alanine

In this article a conformational analysis of the D -alanyl-D -alanine dipeptide, both charged and neutral, has been carried out. The preferred conformations were determined by means of ab initio and semiempirical quantum, together with empirical force field calculations. The AMBER* force field and the 6-31 + G** and 6-31G** ab initio levels give rise to a coincident minimum energy structure, which, on the other hand, differs from that determined by AM1, 3-21 + G, and 3-21G. The solvent effect on the different charged and neutral conformations have been considered through the AMSOL semiempirical method. A quantification regarding the structural similarities between the different dipeptide conformations and the ampicillin has been performed. The results show that the best overlay is attained by the minimum structure energy obtained by using the 6-31 + G** methodology, which presents a planar amidic nitrogen. © 1998 John Wiley & Sons, Inc. Biopoly 45: 119–133, 1998     

The sensitivity of conformational free energies of the alanine dipeptide to atomic site charges

Different atomic point charge sets are obtained for the α_R and C_7.eq conformations of the alanine dipeptide by fitting the charges of each conformation to the respective ab initio electrostatic potential surfaces both individually and simultaneously, in both the united atom and the all-atom representations. Using these charge sets, the sensitivity of the relative conformational aqueous free energies to the atomic site charges is investigated. For this particular system, we find that the solute-water contributions to the conformational free energy differences have a rather weak dependence on site charges; the calculated intramolecular contributions, however, show a rather strong dependence on the atomic site charges. It is suggested that the calculated results for the alanine dipeptide using a single, simultaneously fit set of charges for both conformations are in better agreement with experiments than the calculations carried out with charges determined individually for each conformation. © 1997 John Wiley & Sons, Inc.     

The conformational and vibrational behavior of the inhibitory neuropeptide derived from beta-endorphin

In this study, conformational behavior, structural, and vibrational characterization of the carboxy terminal dipeptide of β-endorphin (glycy-l-glutamine, glycyl-glutamine, beta-endorphin_30-31), which is an inhibitory neuropeptide synthesized from beta-endorphin_1-31 in brain stem regions, has been investigated. The theoretically possible stable conformers were searched by means of molecular mechanics method to determine their energetically preferred conformations. The 360 different conformations were calculated with the φ, Ψ, χ dihedral angles using the Ramachandran maps. The most stable conformation of the title molecule is characterized by the extended backbone shape (e) in the BR conformational range with ?.78 kcal/mol energy. The cis- and trans-dimeric forms of the dipeptide were also formed and energetically preferred conformations of dimers were investigated. The experimental methods (FT-IR, micro-Raman spectroscopies) coupled with quantum chemical calculations based on density functional theory (DFT) have been used to identify the geometrical, energetic, and vibrational characteristics of the dipeptide. The assignment of the vibrational spectra was performed based on the potential energy distribution of the vibrational modes. To investigate the electronic properties, such as nonlinear optical properties, the electric dipole moment, the mean polarizability, the mean first hyperpolarizability, and HOMO–LUMO energy gaps were computed using the DFT with the B3LYP/6-31++G(d,p) basis set combination. The second-order interaction energies were derived from natural bonding orbital analysis. The focus of this study is to determine possible stable conformation on inhibitory neuropeptide and to investigate molecular geometry, molecular vibrations of monomeric and dimeric forms, and hydrogen bonding interactions of glycy-l-glutamine dipeptide.     

Conformational studies on substance P. C-terminal pentapeptide pGlu-Phe-Phe-Gly-Leu-Met NH2

The conformational behavior of the active C-terminal pentapeptide of substance P(SP), pGlu-Phe-Phe-Gly-Leu-Met NH₂ [pGlu-SP(7–11)] was investigated using empirical energy calculations. A sequential approach was used to display the specific contribution of each residue to induce stable conformations of the whole pentapeptide. The most stable conformations include the α_R helix and some partially helical structures; some conformations with glycyl residue in a C₇^eq and C₇^ax configurations (γ and $\text{γ}$ turns) are also favoured. Helical conformations provide a good accessibility of side-chains which play an important role in interacting with the receptor. Fully extended structures and β turns are not specially stable. Such helical stable structures would favour a “lock and key” model of binding.     

Conformational behavior of α,α-dialkylated peptides

The preferred conformations of N-acetyl-N′-methyl amides of some dialkylglycines have been determined by empirical conformational-energy calculations; minimum-energy conformations were located by minimizing the energy with respect to all the dihedral angles of the molecules. The conformational space of these compounds is sterically restricted, and low-energy conformations are found only in the regions of fully extended and helical structures. Increasing the bulkiness of the substituents on the C^α, the fully extended conformation becomes gradually more stable than the helical structure preferred in the cases of dimethylglycine. This trend is, however, strongly dependent on the bond angles between the substituents on the C^α atom: In particular, helical structures are favored by standard values (111°) of the N-C^α-C′ angle, while fully extended conformations are favored by smaller values of the same angle, as experimentally observed, for instance, in the case of α,α-di-n-propylglycine.     

Neural networks as a tool for compact representation of ab initio molecular potential energy surfaces

Ab initio quantum chemical calculations of molecular properties such as, e.g., torsional potential energies, require massive computational effort even for moderately sized molecules, if basis sets with a reasonable quality are employed. Using ab initio data on conformational properties of the cofactor (6R,1′R,2′S)-5,6,7,8-tetrahydrobiopterin, we demonstrate that error backpropagation networks can be established that efficiently approximate complicated functional relationships such as torsional potential energy surfaces of a flexible molecule. Our pilot simulations suggest that properly trained neural networks might provide an extremely compact storage medium for quantum chemically obtained information. Moreover, they are outstandingly comfortable tools when it comes to making use of the stored information. One possible application is demonstrated, namely, computation of relaxed torsional energy surfaces.     

Ab initio refinement of an orbital-centred force field for biomolecules: Orbital localisation and parameterisation of the COP(O2)OC fragment of nucleotides

Classical calculations of conformational potential surfaces, based on simple analytical functions of the interactions between atomic centres, continue to be of considerable importance. However, it has become apparent that not all interactions of importance can be included as interactions between positions which are effectively those of the nuclei. Thus, there has been recent interest in including lone pairs and special functions for interactions involving excited state orbitals. A particularly interesting test case is the COP(O₂)OC fragment of the nucleic acid backbone, which would seem to be the most flexible “hinge point” of polynucleotides when the classical type of calculation is carried out. In contrast, ab initio quantum mechanical calculations show the conformational space of this fragment to be much more restricted. The disagreement is such that it calls into doubt the validity of the bench-top modelling of nucleotide behaviour. In the following study, a variety of ab initio calculations are carried out to localise, in an objective manner, non-core orbitals. Coulombic interactions are introduced between these localised orbitals, with charge parameters optimised to reproduce the total ab initio potential surface. The results imply an interesting disgreement with other authors, concerning the importance of lone pair interactions in the nucleotide backbone, and the origins of this disagreement are analysed in some detail.     

Peptide Free Energy Landscapes Calibrated by Molecular Orbital Calculations

Free energy landscapes of peptide conformations werecalibrated by ab initiomolecular orbital calculations, after enhancedconformational sampling using the multicanonical molecular dynamicssimulations. Three different potentials of mean force for an isolateddipeptide were individually obtained using the conventional force fields,AMBER parm94, AMBER parm96, and CHARMm22. Each potential ofmean force was calibrated based on the umbrella sampling algorithm fromthe adiabatic energy map that was calculated separately by the abinitiomolecular orbital method. All the calibrated potentials of mean forcecoincided well. The calibration was applied to a peptide in explicit water,and the calibrated free energy landscapes did not depend on the force fieldused in conformational sampling, as far as the conformational space waswell sampled.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Application of ab initio molecular-orbital calculations to the structural moieties of carbohydrates. Part VI

Ab initio RHF/4–31G molecular-orbital calculations have been conducted on methoxymethyl formate and methoxymethyl acetate as models for examining the anomeric effect and stereochemistry of 1-O-acetylglycopyranoses. The results indicate that, as with the methyl glycopyranosides, the α-⁴C₁(D) configurations are more stable than the β-⁴C₁(D), except that the energy difference is more dependent on the disposition about the glycosidic bond. The lowest-energy conformations occur with glycosidic torsion-angles of ?  180°, where the anomeric energy is about 4 kcal/mol. There is a secondary energy-minimum at ?  90°, for which the anomeric energy is less, about 2 kcal/mol. This orientation corresponds to the conformation most commonly observed in the crystal structures of peracetylated glycopyranoses. Small differences in the CO single-bond lengths, which are observed experimentally in both the α and β anomers, are reproduced by the theoretical calculations.     

Peptides and Peptoids - A Systematic Structure Comparison

A systematic analysis of the conformational space of the basic structure unit of peptoids in comparison to the corresponding peptide unit was performed based on ab initio MO theory and complemented by molecular mechanics (MM) and molecular dynamics (MD) calculations both in the gas phase and in aqueous solution.The calculations show three minimum conformations denoted as C_7ß, a_D and a that do not correspond to conformers on the gas phase peptide potential energy hypersurface. The influence of aqueous solvation was estimated by means of continuum models. The MD simulations indicate the a_D form as the preferred conformation in solution both in cis and trans peptide bond orientations.     

Conformations in solution of the stereoisomeric,peracetylated aldohexose dimethyl acetals and diethyl dithioacetals

The conformations in solution of acyclic carbohydrate derivatives having four contiguous asymmetric centers in all eight diastercoisomeric forms have been studied by ¹H-n.m.r. spectroscopy. The 250-MHz, ¹H-n.m.r. spectra for solutions in chloroform-d of eight penta-O-acetylaldohexose dimethyl acetals, and the corresponding diethyl dithioacetals, furnished a complete set of chemical shifts and proton-proton spin-couplings that are interpreted in terms of conformational compositions at room temperature. The galacto and manno derivatives adopt planar, extended conformations, whereas the other six stereoisomers all adopt one or more non-extended (“sickle”) conformations. The results are interpreted on the basis of the avoidance of parallel 1,3-interactions of substituents. The conformational assignments are correlated with observations made previously for aldopentose analogs. An assessment is made of the extent to which valid conformational predictions may be advanced for four-center, and longer, asymmetrically substituted chains, based on observations made for shorter-chain analogs.     

Towards the Understanding of the Folding of Methylene Units in the Glutamine Residue

The conformational preferences of the methylenic sequence in the side chain of the glutamine residue were investigated by ab initio and semi-empirical quantum mechanical calculations and examination of both the Brookhaven Protein Databank and Cambridge Structural Data Base. The results were analysed on the basis of our previous findings about the folding of methylene groups in aliphatic segments. Both energy calculations and the crystallographic structure of small peptides indicate that methylene units of the glutamine residue tend to fold in a gauche conformation. In contrast, such groups usually adopt an all-trans conformation in proteins due basically to the entropic and solvent contributions. These results have been demonstrated by computing the entropic correction to the free energy and evaluating the solvent effects through SCRF calculations     

Study of the Tautomerism by Experimentally and Theoretically Estimated 13C and 15N Chemical Shifts

2-Substituted 5-Me-7OH-1,2,4-Triazolo[1,5-a]pyrimidines can exist as different tautomeric forms. The tautomeric equilibrium of these compounds were determined by various NMR spectroscopic methods and quantum-chemical calculations. The tautomeric equilibria of 2-substituted 1,2,4-Triazolopyrimidines depend on the used solvent. In order to take solvent effects into account, ab initio calculations using the SCI-PCM method were carried out. The solvation energy depends on the polarity of the substituents in position 2. The comparison of chemical shifts determined by experimental methods and theoretical ab initio methods were used to definitely find out the the tautomeric equilibria of these compounds.     

A Molecular Mechanics and Database Analysis of the Structural Preorganization and Activation of the Chromophore-Containing Hexapeptide Fragment in Green Fluorescent Protein

Abstract

We propose that heterologous posttranslational chromophore formation in green fluorescent protein (GFP) occurs because the chromophore-forming amino acid residues 65SYG67 are preorganized and activated for imidazolinone ring formation. Based on extensive molecular mechanical conformational searching of the precursor hexapeptide fragment (64FSYGVQ69), we suggest that the presence of low energy conformations characterized by short contacts (～3Å) between the carbonyl carbon of Ser65 and the amide nitrogen of Gly67 accounts for the initial step in posttranslational chromophore formation. Database searches showed that the tight turn required to establish the key short contact is a unique structural motif that is rarely found, except in other FSYG tetrapeptide sequences. Additionally, ab initio calculations demonstrated that an arginine side chain can hydrogen bond to the carbonyl oxygen of Ser65, activating this group for nucleophilic attack by the nearby lone pair of the Gly67 amide nitrogen. We propose that GFP chromophore-formation is initiated by a unique combination of conformational and electronic enhancements, identified by computational methods.     

Monte Carlo calculations on the conformations of models for the glycopeptide linkage of glycoproteins

In order to search for probable conformations of the peptide, the amino acid side chain, and the carbohydrate linkage in glycoproteins, conformational energy surfaces of glycopeptide model compounds were studied by Monte Carlo methods using the Metropolis algorithm. The potential energies were composed of empirical energy functions which include nonbonded interactions, electrostatics, hydrogen bonding, and torsional energies specified by parameters which have been used for peptides. Calculations were performed on 1-N-acetyl-2-acetamido-beta-D-glucopyranosyl amine and the glycosylated dipeptide N-acetyl-delta-N-(2-acetamido-beta-D-glucopyranosyl)-L-asparaginyl-N'-methyl amide as models for N-glycosylated peptides and on methyl-2-acetamido-alpha-D-galactopyranoside as well as the glycosylated dipeptides N-acetyl-gamma-O-(2-acetamido-alpha-D-galactopyranosyl)-L-threonyl-N'-methyl amide and its seryl analog as models for O-glycosylated glycoproteins. The probable conformations of these compounds were analyzed by single-angle probability tables and by two-dimensional conformation density maps projected from the Markov chains which contained up to six independently varied conformational dihedral angles. The presence of high barriers to rotation required the use of search strategies which resulted in a rather low acceptance rate for new conformations in the Metropolis algorithm in order to avoid trapping of the Markov chain in local energy minima. This problem contributed to the failure of these calculations to attain complete convergence to the thermodynamic limit for the glycosylated dipeptide models in which six dihedral angles were independently varied. Analysis of the results shows that the conformational space available to the highly crowded axial glycosides of the alpha-O-GalNAc type is much more restricted than that for the N-asparaginyl glycopeptides. The most probable conformation for the O-glycosylated peptides is is a beta-turn while N-glycosylated peptides may be either in a beta-turn or an extended conformation.