New Methodology for Computer-Aided Modelling of Biomolecular Structure and Dynamics 2. Local Deformations and Cycles

Abstract

A new methodology for the conformational modelling of biomolecular systems (1) is extended to local deformations of chain molecules and to flexible molecular rings. It is shown that these two cases may be reduced to considering an equivalent molecular model with a regular tree-like topology. A simple procedure is developed to analyze any flexible rings (the five- and six-membered suguar rings of carbohydrates and nucleic acids, in particular) and local deformation regions by energy minimization. Dynamic equations are also derived for such molecular systems. As a result, a unified approach is proposed for the efficient energy minimization and simulation of dynamic behavior of multimolecular systems having any set of variable internal coordinates, local deformation regions and cycles. Advantages and domains of applicability of the approach are discussed.     

Molecular Dynamics and Free Energy Calculations of the B and Z Forms of C8-Arylguanine Modified Oligonucleotides

Abstract

Arylhydrazines found in the mushroom Agaricus bisporus have been shown to be carcinogenic. Upon metabolic activation, arylhydrazines are transformed into aryl radicals, forming 8-arylpurines, which may play a role in arylhydrazine carcinogenesis. These adducts are poorly read and inhibit chain extension but do alter the conformational preferences of oligonucleotides. We have shown that C8-phenylguanine modification of d(CGCGCG*CGCG) (G*= 8-phenylguanine) stabilizes it in the Z-DNA conformation (B/Z-DNA=1:1, 200 mM NaCl, pH 7.4). Here we have conducted molecular dynamics and free energy calculations to determine the sources(s) of these conformational affects and to predict the affect of the related C8- tolyl and C8-hydroxymethylphenyl guanine adducts on B/Z-DNA equilibrium. Force field parameters for the modified guanines were first developed using Guassian98 employing the B3LYP method and the standard 6–31G* basis set and fit to the Cornell 94 force field with RESP. Molecular dynamics simulations and free energy calculations, using the suite of programs contained in Amber 6 and 7 with the Cornell 94 force field, were used to determine the structural and thermodynamic properties of the DNA. The principal factors that drive conformation are stacking of the aryl group over the 5′-cytosine in the phenyl and tolyl modified oligonucleotides while hydrogen bonding opposes stacking in the hydroxymethylphenyl derivative. The phenyl and tolyl-modified DNA's favored the Z-DNA form as did the hydroxymethylphenyl derivative when hydrogen bonding was not present. The B-DNA conformation was preferred by the unmodified oligonucleotide and by the hydroxymethylphenyl-modified oligonucleotide when hydrogen bonding was considered. Z-DNA stability was not found to directly correlated with carcinogenicity and additional biological factors, such as recognition and repair, may also need to be considered in addition to Z-DNA formation.     

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.     

Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides

Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl α-d-glucopyranoside and methyl α-d-galactopyranoside, employing the glycam06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the ω (O-5–C-5–C-6–O-6) and χ (C-3–C-4–O-4–HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4–HO-4 and O-6–HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.     

Radioprotective effect of hydrogen in cultured cells and mice

Abstract

It has been demonstrated that hydrogen can selectively reduce hydroxyl and peroxynitrite in vitro. Since most of the ionizing radiation-induced cellular damage is caused by hydroxyl radicals, this study was designed to test the hypothesis that hydrogen may be an effective radioprotective agent. This paper demonstrates that treating cells with hydrogen before irradiation could significantly inhibit ionizing irradiation(IR)-induced Human Lymphocyte AHH-1 cells apoptosis and increase cells viability in vitro. This paper also shows that hydrogen can protect gastrointestinal endothelia from radiation-induced injury, decrease plasma malondialdehyde (MDA) intestinal 8-hydroxydeoxyguanosine (8-OHDG) levels and increase plasma endogenous antioxidants in vivo. It is suggested that hydrogen has a potential as an effective and safe radioprotective agent.     

Spatial configurations of polynucleotide chains. 3. Polydeoxyribonucleotides

The virtual bond scheme set forth in preceding papers for treating the average properties of polyriboadenylic acid (poly rA) is here applied to the calculation of the unperturbed mean-square end-to-end distance of polydeoxyriboadenylic acid (poly dA). The modifications in structure and in charge distribution resulting from the replacement of the hydroxyl group at C_2′ in the ribose residue by hydrogen in deoxyribose produce only minor modifications in the conformational energies associated with the poly dA chain as compared to those found for poly rA. The main difference is manifested in the energy associated with rotations about the C_3′–O_3′ bond of the deoxyribose residue in the C_2′-endo conformation; accessible rotations are confined to the range between 0° and 30° relative to the trans conformation, whereas in the ribose unit the accessible regions comprise two ranges centered at approximately 35° and 85°. The characteristic ratio 〈r²〉0/nl² calculated on the basis of the conformational energy estimates is ≈9 for the poly dA chain with all deoxyribose residues in the C_3′-endo conformation and ≈21 with all residues in the C_2′-endo form. Satisfactory agreement is achieved between the theoretical values and experimental results on apurinic acid by treating the poly dA chain as a random copolymer of C_3′-endo and C_2′-endo conformational isomers present in a ratio of ～1 to 9.     

Towards the simulation of biomolecules: optimisation of peptide-capped glycine using FFLUX

Abstract

The optimisation of a peptide-capped glycine using the novel force field FFLUX is presented. FFLUX is a force field based on the machine-learning method kriging and the topological energy partitioning method called Interacting Quantum Atoms. FFLUX has a completely different architecture to that of traditional force fields, avoiding (harmonic) potentials for bonded, valence and torsion angles. In this study, FFLUX performs an optimisation on a glycine molecule and successfully recovers the target density-functional-theory energy with an error of 0.89 ± 0.03 kJ mol^?1. It also recovers the structure of the global minimum with a root-mean-squared deviation of 0.05 Å (excluding hydrogen atoms). We also show that the geometry of the intra-molecular hydrogen bond in glycine is recovered accurately.     

Amyloglucosidase-catalyzed synthesis of eugenyl and curcuminyl glycosides

Glycosylation of the phenolic hydroxyl group of the phenyl propanoid systems, eugenol 1 and curcumin 2, using an amyloglucosidase from Rhizopus and a β-glucosidase from sweet almonds together with carbohydrates (d-glucose 3, d-mannose 4, maltose 5, sucrose 6 and d-mannitol 7) in di-isopropyl ether produced glycosides at 7–52% yields in 72 h. Spectral studies indicated that the reaction occurred between the phenolic OH groups and C-1 and/or 6-O-groups of the carbohydrates with curcumin exhibiting bis glycosylation.     

A salt-bridge switch in the molecular recognition between RS receptor and RGD ligand from the ABEEM σπ molecular dynamics simulations

Abstract

How the receptor and ligand recognise each other is a challenging subject in explaining the mechanism of recognition at the molecular level. As a starting point, here, a synthesised RS receptor and its RGD ligand were investigated as a proper model to simulate their recognition process in terms of ABEEMσπ/MM polarisable force field. It is found that a switch of forming up a salt bridge in the ligand triggers the recognition of the receptor and ligand. Through the salt-bridge switch that undergoes several cycles from on-state with parallel hydrogen bonds to off-state with bifurcated hydrogen bonds, the active site of ligand can flex easily to interact with the active site of the receptor. In addition, the water molecules form a decisive bridge connecting the active sites of the bound system. The salt-bridge switch and water-mediated movement are cooperative as the important factors for the receptor-ligand recognition. In addition, the properties, such as binding free energy, conformational flexibility and solvent accessible surface area have been calculated to provide adequate evidence for the whole recognition process. According to the simulation, a detailed mechanism was derived involving diffusion, a switch triggered cooperative water-mediated movement, and conformational folding, for the flexible recognition.     

Conformational Studies of Dextran [α-(1⇒6)-D-Glucan]

Abstract

A theoretical conformational study of dextran, a (l?6)-linked α-D-glucan polysaccharide, has been made to allow an explicit comparison with earlier results on pustulan, the corresponding (1 ?6)-linked β-D-glucan. The nonbonded, torsional and hydrogen bond contributions to potential energy were calculated as a function of rotational angles φ, ψ, and ω The (φ, ψ, ω)-space of the disaccharide and of helices contain many local energy minima with very small energy differences. A comparison of (1?6)-α-D-glucans with (1?6)-β-D-glucans indicates significant differences in conformational behavior. Specifically, our results shed light on the fact that dextran does not gel, whereas pustulan does. The difference in tendency to gel may be related to the fact that dextran has no particularly favored conformations with structural regularity whereas pustulan does.     

The effect of hydration upon the conformation and dynamics of neocarrabiose,a repeat unit of β-carrageenan

Molecular mechanics calculations have been performed for the disaccharide neocarrabiose, one of the repeat units of β-carrageenan, as a general model for the (1 → 3)-linkage in the carrageenans. An adiabatic conformational energy map for this molecule has been prepared by constrained energy minimization and compared to previously reported relaxed maps. Neither the experimentally determined crystal structure of neocarrabiose nor the fiber diffraction conformation of β-carrageenan is a low energy conformation on the relaxed Ramachandran map. Molecular dynamics simulations in vacuum produced trajectories consistent with this relaxed vacuum surface. However, a simulation with explicitly included solvent water molecules produced a trajectory that remained in the region of the two experimental structures. This dramatic solvation effect is apparently the result of the breaking of an interring hydrogen bond between the O2 hydroxyl groups of neocarrabiose as both groups hydrogen bond to solvent. © 1996 John Wiley & Sons, Inc.     

Polarization Force Fields for Peptides Implemented in ECEPP2 and MM2

Abstract

The empirical conformational energy program for peptides (ECEPP2) and molecular mechanics (MM2) have been used for the simulation of the For-Gly-NH₂ backbone. I propose two different methods for the calculation of the polarization energy term: the polarization procedure by non-interacting induced dipoles (NID) which assumes scalar isotropic point polarizabilities and the polarization scheme by interacting induced dipoles (ID) which calculates tensor effective anisotropic point polarizabilities (method of Applequist). I present a comparative study of ECEPP2 and MM2 + polarization. I discuss molecular mechanics results including the total energy differences, partitional analyses of the total steric energies and torsion dihedral angles. The γ global and the α, β and Δ local minima are stabilized by intramolecular hydrogen bonds. Although ECEPP2-based calculations rather under or over-estimate the relative energy of some local minima, the ID polarization energy term represents a significant correction to the total relative energy.     

Alcohols as Replacements of the Central Amide in β-Turns, Synthesis of Pro-Aib Hydroxyethylene Isostere and Analysis in Model β-Turn Peptides

Pro-Aib hydroxyethylene isosteres (S,R)- and (S,S)-7 were synthesized by cascade addition of 2-methyl-1-propenylmagnesium bromide to Boc-Pro-OMe in the presence of CuCN, followed by ketone reduction and olefin oxidation. By protecting the amine and hydroxyl groups in an oxazolidinone ring, hydroxyethylene isosteres 7 were successfully incorporated into Boc-Phe-Pro- ψ -[CH-(OH)-CH₂]-Aib-NHBn(α -Me) (S,R)-and (S,S)-11, which were characterized by ¹H NMR and IR spectroscopy. Examination of the NOESY spectra and the influence of solvent changes on the chemical shifts of the amide and carbamate proton signals for (S,R)-and (S,S)-11 indicated that both hydroxyethylene isosteres could adopt compact turn structures. The alcohol appears to act as a hydrogen donor in a seven-membered ring intramolecular hydrogen bond. In addition, analysis of the respective peptide (S,S)-16, in which the hydroxyl group was masked as a methyl ether, showed that the turn conformation was disrupted, and indicated the importance of the alcohol as a hydrogen-bond donor for turn stability. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

On the Bioactive Conformation of a Small Peptide and its Set of Thermodynamically Accessible Conformations

Abstract

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics framework using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.     

Location of proline derivatives in conformational space. I. Conformational calculations; optical activity and NMR experiments

In order to develop methods of analysis applicable to the determination of the conformation of biological polymers in solution, a series of proline derivatives was studied. The steric constraints of the pyrrolidine ring limit these compounds to a relatively small set of conformations. This set was further reduced by eliminating conformations with large computed conformational energy. Computations revealed that the conformational energy of the proline derivatives fits into one of three classes, depending on the bulk and the polarity of the C-terminal group. Three analogous classes of optical activity were observed. The optical activity data were analyzed in terms of conformations computed to be of low energy. In some cases qualitative theoretical considerations enabled molecular groups to be located. For example, solvent-dependent isomerization of the carboxyl hydrogen of N-acetyl-L -proline was detected. Nuclear magnetic resonance provided an experimental measure of the fraction of molecules which had cis unsymmetrically-substituted tertiary amide groups. This information aided and confirmed the other measures of molecular conformation.     

Electrostatic interaction-mediated conformational changes of adipocyte fatty acid binding protein probed by molecular dynamics simulation

Abstract

Adipocyte fatty acid binding protein (A-FABP) is a potential drug target for treatment of diabetes, obesity and atherosclerosis. Molecular dynamics (MD) simulations, principal component (PC) analysis and binding free energy calculations were combined to probe effect of electrostatic interactions of residues R78, R106 and R126 with inhibitors ZGB, ZGC and IBP on structural stability of three inhibitor/A-FABP complexes. The results indicate that mutation R126A produces significant influence on polar interactions of three inhibitors with A-FABP and these interactions are main force for driving the conformational change of A-FABP. Analyses on hydrogen bond interactions show that the decrease in hydrogen bonding interactions of residues R126 and Y128 with three inhibitors and the increase in that of K58 with inhibitors ZGC and IBP in the R126A mutated systems mostly regulate the conformational changes of A-FABP. This work shows that R126A can generate a significant perturbation on structural stability of A-FABP, which implies that R126 is of significance in inhibitor bindings. We expect that this study can provide a theoretical guidance for design of potent inhibitors targeting A-FABP.

Communicated by Ramaswamy H. Sarma     

Conformational changes due to vicinal glycosylation: The branched α-l-Rhap(1–2)[β-d-Galp(1–3)]-β-d-Glc1-OMe trisaccharide compared with its parent disaccharides*

Conformations of the α-l -Rhap(1-2)-β-d -Glc1-OMe and β-d -Galp(1-3)-β-d -Glc1-OMe disaccharides and the branched title trisaccharide were examined in DMSO-d₆ solution by ¹H-nmr. The distance mapping procedure was based on rotating frame nuclear Overhauser effect (NOE) constraints involving C- and O-linked protons, and hydrogen-bond constraints manifested by the splitting of the OH nmr signals for partially deuteriated samples. An “isotopomer-selected NOE” method for the unequivocal identification of mutually hydrogen-bonded hydroxyl groups was suggested. The length of hydrogen bonds thus detected is considered the only one motionally nonaveraged nmr-derived constraint. Molecular mechanics and molecular dynamics methods were used to model the conformational properties of the studied oligosaccharides. Complex conformational search, relying on a regular Φ,Ψ-grid based scanning of the conformational space of the selected glycosidic linkage, combined with simultaneous modeling of different allowed orientations of the pendant groups and the third, neighboring sugar residue, has been carried out. Energy minimizations were performed for each member of the Φ,Ψ grid generated set of conformations. Conformational clustering has been done to group the minimized conformations into families with similar values of glycosidic torsion angles. Several stable syn and anti conformations were found for the 1→2 and 1→3 bonds in the studied disaccharides. Vicinal glycosylation affected strongly the occupancy of conformational states in both branches of the title trisaccharide. The preferred conformational family of the trisaccharide (with average Φ,Ψ values of 38°, 17° for the 1→2 and 48°, 1° for the 1→3 bond, respectively) was shown by nmr to be stabilized by intramolecular hydrogen bonding between the nonbonded Rha and Gal residues. © 1998 John Wiley & Sons, Inc. Biopoly 46: 417–432, 1998     

A crystallographic and theoretical study of an (E)-2-Hydroxyiminoethanone derivative: prediction of cyclooxygenase inhibition selectivity of stilbenoids by MM-PBSA and the role of atomic charge

We recently reported that the hydroxyiminoethanone derivative, (E)-OXM, behaves as a highly selective COX-1 inhibitor (COX-1 SI = 833), and also an interesting scaffold with unique characteristics. In the current study, a comprehensive crystallographic and computational study was performed to elucidate its conformational stability and pharmacological activity. Its conformational energy was studied at the B3LYP/6-311G** level of theory and compared to the single-crystal X-ray diffraction data. In addition, computational studies of three structurally different stilbenoid derivatives used as selective COX-1 or COX-2 inhibitors were undertaken to predict their COX selectivity potentials. Flexible docking was performed for all compounds at the active site of both COX-1 and COX-2 enzymes by considering some of the key residues as flexible during the docking operation. In the next step, molecular dynamic simulation and binding free energy calculations were performed by MM-PBSA. Final results were found to be highly dependent on the atomic charges of the inhibitors and the choice of force field used to calculate the atomic charges. The binding conformation of the hydroxyiminoethanone derivative is highly correlated with the type of COX isoform inhibited. Our predictive approach can truly predict the cyclooxygenase inhibition selectivity of stilbenoid inhibitors.     

A CFF 91-based Continuum Solvation Model: Solvation Free Energies of Small Organic Molecules and Conformations of the Alanine Dipeptide in Solution

Abstract

For molecular mechanics simulations of solvated molecules, it is important to use a consistent approach for calculating both the force field energy and the solvation free energy. A continuum solvation model based upon the atomic charges provided with the CFF91 force field is derived. The electrostatic component of the solvation free energy is described by the Poisson-Bolzmann equation while the nonpolar comonent of the solvation energy is assumed to be proportional to the solvent accessible surface area of the solute. Solute atomic radii used to describe the interface between the solute and solvent are fitted to reproduce the energies of small organic molecules. Data for 140 compounds are presented and compared to experiment and to the results from the well-characterized quantum mechanical solvation model AM1-SM2. In particular, accurate results are obtained for amino acid neutral analogues (mean unsigned error of 0.3 kcal/mol). The conformational energetics of the solvated alanine dipeptide is discussed.