Flexibility of the pyranose ring in α- and β-D-glucoses

Conformational energies of α- and β-D -glucopyranoses were computed by varying all the ring bond angles and torsional angles using semiempirical potential functions. Solvent accessibility calculations were also performed to obtain a measure of solvent interaction. The results indicate that the ⁴C₁ (D ) chair is the most favored conformation, both by potential energy and solvent accessibility criteria. The ⁴C₁ (D ) chair conformation is also found to be somewhat flexible, being able to accommodate variations up to 10° in the ring torsional angles without appreciable change in energy. Observed solid-state conformations of these sugars and their derivatives lie in the minimum-energy region, suggesting that the substituents and crystal field forces play a minor role in influencing the pyranose ring conformation. Theory also predicts the variations in the ring torsional angles, i.e., CCCC < CCCO < CCOC, in agreement with the experimental results. The boat and twist-boat conformations are found to be at least 5 kcal mol^?1 higher in energy compared to the ⁴C₁ (D ) chair, suggesting that these forms are unlikely to be present in a polysaccharide chain. The ¹C₄ (D ) chair has energy intermediate between that of the ⁴C₁ (D ) chair and that of the twist-boat conformation. The calculated energy barrier between ⁴C₁ (D ) and ¹C₄ (D ) conformations is high—about 11 kcal mol^?1.     

Molecular dynamics simulation and nmr study of a blood group H trisaccharide

Molecular dynamics simulations in vacuum and solution have been carried out on 2′-α-L -fucosyllactitol, a model for blood group H in conjunction with two-dimensional nmr measurements on the same compound. Three independent starting conformations for the dynamics were chosen from low energy conformations obtained by a ?/ψ grid search. Nine 5 ns vacuum simulations of the trisaccharide were performed, employing three different ways to treat electrostatic interactions for each starting conformation: distance-dependent dielectric with ε = r, constant dielectric with ε = 1, or constant dielectric with ε = 80. In vacuum, transitions of ? and ψ for the α-L -Fuc-(1 → 2)-β-D -Gal element occur in a cooperative manner. The virtual distance obtained for H1 in fucose to H2 in galactose from nuclear Overhauser effect spectroscopy experiments agree with one of the conformations of the trisaccharide in one of the three 100 ps aqueous simulations (?/ψ ca. ?100°/150°), indicating this may be a dominant solution conformation. The rms fluctuations of the ?- and ψ-dihedral angles were ～ 10° for a conformational state, both in the vacuum and the aqueous simulations. For the simulations in vacuum, the agreement with experimental NOE data is reasonable when a constant dielectric of 1 is used (major conformers having ?/ψ ca. ?100°/150° and ?140°/100°), whereas the agreement was poor with a constant dielectric of 80. Translational diffusion coefficients calculated from the simulation of the oligosaccharides were 0.12–0.18 × 10^?5 cm²/s and from nmr measurements 0.27 × 10^?5 cm²/s. © 1994 John Wiley & Sons, Inc.     

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

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

Conformational equilibrium of unsulphated iduronate in heparan sulphate tetrasaccharides

Proton-proton coupling constants for terminal -l-iduronate residues in tetrasaccharides obtained from heparan sulphates by complete nitrous acid deaminative cleavage were shown to vary with experimental conditions. It is proposed that the iduronate residue is in a conformational equilibrium between the¹C₄ chair and either the²S_o skewboat or possibly the²H₃ half-chair conformers. It was not possible to discriminate between the two non-chair forms empirically. The position of the equilibrium is sensitive to temperature, pH and sulphation of neighbouring residues. The likelihood of iduronate residues within glycosaminoglycans existing in the⁴C₁ conformer in addition to the¹C₄ and²S_o forms is discussed.     

Glycan flexibility: insights into nanosecond dynamics from a microsecond molecular dynamics simulation explaining an unusual nuclear Overhauser effect

An atomistic all-atom molecular dynamics simulation of the trisaccharide β-d-ManpNAc-(1→4)[α-d-Glcp-(1→3)]-α-l-Rhap-OMe with explicit solvent molecules has been carried out. The trisaccharide represents a model for the branching region of the O-chain polysaccharide of a strain from Aeromonas salmonicida. The extensive MD simulations having a 1-μs duration revealed a conformational dynamics process on the nanosecond time scale, that is, a ‘time window’ not extensively investigated for carbohydrates to date. The results obtained from the MD simulation underscore the predictive power of molecular simulations in studies of biomolecular systems and also explain an unusual nuclear Overhauser effect originating from conformational exchange.     

Preferences of AAA/AAG codon recognition by modified nucleosides, τm5s2U34 and t6A37 present in tRNALys

Deficiency of 5-taurinomethyl-2-thiouridine, τm⁵s²U at the 34th ‘wobble’ position in tRNA^Lys causes MERRF (Myoclonic Epilepsy with Ragged Red Fibers), a neuromuscular disease. This modified nucleoside of mt tRNA^Lys, recognizes AAA/AAG codons during protein biosynthesis process. Its preference to identify cognate codons has not been studied at the atomic level. Hence, multiple MD simulations of various molecular models of anticodon stem loop (ASL) of mt tRNA^Lys in presence and absence of τm⁵s²U₃₄ and N⁶-threonylcarbamoyl adenosine (t⁶A₃₇) along with AAA and AAG codons have been accomplished. Additional four MD simulations of multiple ASL mt tRNA^Lys models in the context of ribosomal A-site residues have also been performed to investigate the role of A-site in recognition of AAA/AAG codons. MD simulation results show that, ASL models in presence of τm⁵s²U₃₄ and t⁶A₃₇ with codons AAA/AAG are more stable than the ASL lacking these modified bases. MD trajectories suggest that τm⁵s²U recognizes the codons initially by ‘wobble’ hydrogen bonding interactions, and then tRNA^Lys might leave the explicit codon by a novel ‘single’ hydrogen bonding interaction in order to run the protein biosynthesis process smoothly. We propose this model as the ‘Foot-Step Model’ for codon recognition, in which the single hydrogen bond plays a crucial role. MD simulation results suggest that, tRNA^Lys with τm⁵s²U and t⁶A recognizes AAA codon more preferably than AAG. Thus, these results reveal the consequences of τm⁵s²U and t⁶A in recognition of AAA/AAG codons in mitochondrial disease, MERRF.     

Consequences of Sequential Ca2+ Occupancy for the Structure and Dynamics of CalbindinD9K: Computational Simulations and Comparison to Experimental Determinations in Solution

Abstract

Molecular dynamics (MD) simulations of the structures of calbindin_D9K (CAB) with different occupancies of the two Ca²⁺ binding sites were carried out to gain insight into structural and energetic consequences of sequential Ca²⁺ binding. The aim of the study is to identify effects of Ca-binding site occupancy that relate to the properties and functions of Ca-binding proteins containing EF-hand motifs. Two different models of solvation were employed, one defined by a linear, distance dependent dielectric permittivity (ε = r) and inclusion only of the 36 crystallographically observed water molecules, and the other with the protein immersed in a 9Å shell of explicit waters and ε = 1. Experimental results from x-ray crystallography, and insights from NMR and from measurements of hydrogen exchange rates in these systems served as tests and guides for assessing the quality, validity and mechanistic interpretation of the results from the computational study. The results of the MD simulations agree very well with earlier experimental observations that the structure of calbindin_D9k is rather insensitive to removal of Ca²⁺, and indicate that this insensitivity is not dependent on the order in which the ions are removed. The calculated values of the electrostatic potentials at the Ca²⁺ binding sites are very similar, in agreement with the small differences in the measured microscopic binding constants. Details of the dynamic mechanisms of molecular flexibility revealed by the MD simulations are also in good agreement with experimental findings, including the local properties identified from comparisons of hydrogen exchange rates in various parts of the structures of sequentially occupied forms of CAB. Estimation of the changes in configurational entropy from the rms fluctuations in the structures of CAB at various levels of Ca²⁺ occupancy in the EF-hands, supports earlier suggestions relating the dynamic properties of the protein to the observed cooperativity in the binding of Ca²⁺. The computational approaches and the results of the MD simulations are evaluated in relation to the study of effects of Ca²⁺ occupancy in calmodulin and troponin C where ion binding determines function and is known to trigger significant changes in structural and dynamic properties.     

The crystal structure of rice (Oryza sativa L.) Os4BGlu12, an oligosaccharide and tuberonic acid glucoside-hydrolyzing β-glucosidase with significant thioglucohydrolase activity

Rice Os4BGlu12, a glycoside hydrolase family 1 (GH1) β-glucosidase, hydrolyzes β-(1,4)-linked oligosaccharides of 3–6 glucosyl residues and the β-(1,3)-linked disaccharide laminaribiose, as well as certain glycosides. The crystal structures of apo Os4BGlu12, and its complexes with 2,4-dinitrophenyl-2-deoxyl-2-fluoroglucoside (DNP2FG) and 2-deoxy-2-fluoroglucose (G2F) were solved at 2.50, 2.45 and 2.40 Å resolution, respectively. The overall structure of rice Os4BGlu12 is typical of GH1 enzymes, but it contains an extra disulfide bridge in the loop B region. The glucose ring of the G2F in the covalent intermediate was found in a ⁴C₁ chair conformation, while that of the noncovalently bound DNP2FG had a ¹S₃ skew boat, consistent with hydrolysis via a ⁴H₃ half-chair transition state. The position of the catalytic nucleophile (Glu393) in the G2F structure was more similar to that of the Sinapsis alba myrosinase G2F complex than to that in covalent intermediates of other O-glucosidases, such as rice Os3BGlu6 and Os3BGlu7 β-glucosidases. This correlated with a significant thioglucosidase activity for Os4BGlu12, although with 200- to 1200-fold lower k_cat/K_m values for S-glucosides than the comparable O-glucosides, while hydrolysis of S-glucosides was undetectable for Os3BGlu6 and Os3BGlu7.     

Computer Aided Study of Ligand Binding With Catalytic Domain of Avian Sarcoma Virus Integrase and Its Ligand Binding Loops

Abstract

We have carried out 1 nanosecond (ns) Molecular Dynamics (MD) simulations of the drug Y3 (4-acetylamino-5-hydroxynaphthalene-2, 7-disulfonic acid) complexed with catalytic domain of Avian sarcoma virus Integrase (ASV-IN), both in vacuum and in the presence of explicit solvent. Starting models were obtained on the basis of PDB co-ordinates (1A5X) of ASV-IN-Y3 complex, by Lubkowski et al [1]. Mn²⁺ cation was present in the active site. To neutralize the positive charge in the presence of explicit solvent, eight Cl^? anions were added. Energy Minimization (EM) and MD simulations, for both the systems, were carried out using Sander's module of AMBER5.0 [2] with all atom force field. Analysis of ligand- protein interaction in both environments is discussed in the paper. We also carried out 1 ns MD simulation on two flexible loops—L1 (Gly54-Gln62) and L2 (Trp138-Met155) playing crucial role in interaction of IN with the drug [3], under differing environmental conditions (vacuum, aqueous and organic solvent methanol). Comparison of the conformational changes in the loops, monomer and dimer is presented in the paper. Our results showed that the conformation of the loop region was closest to crystallographic data in case of monomer and constrained loops in aqueous environment. However, the dimer in vacuum was more stable than monomer. The β sheet structure of the monomer in aqueous environment was unstable. Latter also took long time for equilibration. The box formed by loops L1 and L2 from two sub units IINA and INB) of the dimer satisfies prerequisites for ligand recognition site and seems to be the functional biological unit.     

Relating side-chain mobility in proteins to rotameric transitions: Insights from molecular dynamics simulations and NMR

The dynamic aspect of proteins is fundamental to understanding protein stability and function. One of the goals of NMR studies of side-chain dynamics in proteins is to relate spin relaxation rates to discrete conformational states and the timescales of interconversion between those states. Reported here is a physical analysis of side-chain dynamics that occur on a timescale commensurate with monitoring by ²H spin relaxation within methyl groups. Motivated by observations made from tens-of-nanoseconds long MD simulations on the small protein eglin c in explicit solvent, we propose a simple molecular mechanics-based model for the motions of side-chain methyl groups. By using a Boltzmann distribution within rotamers, and by considering the transitions between different rotamer states, the model semi-quantitatively correlates the population of rotamer states with ‘model-free’ order parameters typically fitted from NMR relaxation experiments. Two easy-to-use, analytical expressions are given for converting S²_axis’ values (order parameter for C–CH₃ bond) into side-chain rotamer populations. These predict that S²_axis’ values below 0.8 result from population of more than one rotameric state. The relations are shown to predict rotameric sampling with reasonable accuracy on the ps–ns timescale for eglin c and are validated for longer timescales on ubiquitin, for which side-chain residual dipolar coupling (RDC) data have been collected.     

Computer simulation of the cyclodextrin–phenylalanine complex

The results of the molecular dynamics simulations of the complexes of α-cyclodextrin-l-phenylalanine and β-cyclodextrine- L- phenylalanine in vacuo and in aqueous solution are presented. The trajectories of the insertion angle, rotation of the aromatic ring of the phenylalanine inside the macrocycle and the dihedral angle χ² (Cα–Cβ–Cγ–C_D2) describing the relative movement of the aromatic ring with respect to the polar region give detailed information of the dynamics of the complexes. It is found that the complex with α-cyclodextrin in water is not stable, in agreement with experimental data, while in all other situations studied the complex is stable within the computational limits. Comparing the different cases and the experimental evidence it comes out that a simulation of the complexes without an explicit treatment of the solvent gives unreliable results.     

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Boat conformations synthesis, NMR spectroscopy, and molecular dynamics of methyl 4,6-O-benzylidene-3-deoxy-3-phthalimido-alpha-D-altropyranoside derivatives

Addition of the elements of phthalimide to methyl 2,3-anhydro-4,6-O-benzylidene-alpha-D-mannopyranoside (1) under fusion conditions has yielded methyl 4,6-O-benzylidene-3-deoxy-3-phthalimido-alpha-D-altropyranoside (2). The conformation of the pyranose ring of 2 has been shown to be non-chair by 1H NMR spectroscopy, in contrast to the conformations of related derivatives having smaller substituents at C-3. Molecular dynamics simulations of 2 in explicit chloroform-d solvent have indicated four principal conformational possibilities. Of these, the 7C5/1S5 chair/skew boat form 2d has the lowest potential energy, and is largely consistent with the observed vicinal 1H-1H NMR coupling constants.     

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.     

Crystal structures of diketopiperazines containing α-aminoisobutyric acid: Cyclo(Aib-Aib) and cyclo(Aib-L-Ile)

The crystal and molecular structures of two α-aminoisobutyric acid (Aib)-containing diketopiperazines, cyclo(Aib-Aib) 1 and cyclo(Aib-L -Ile) 2 , are reported. Cyclo(Aib-Aib) crystallizes in the space group P1 with a = 5.649(3), b = 5.865(2), c = 8.363(1), α = 69.89(6), β = 113.04(8), γ = 116.0(3), and Z = 1, while 2 occurs in the space group P2₁2₁2₁ with a = 6.177(1), b = 10.791(1), c = 16.676(1), and Z = 4. The structures of 1 and 2 have been refined to final R factors of 0.085 and 0.086, respectively. In both structures the diketopiperazine ring shows small but significant deviation from planarity. A very flat chair conformation is adopted by 1, in which the C^α atoms are displaced by 0.07 Å on each side of the mean plane, passing through the other four atoms of the ring. Cyclo(Aib-Ile) favors a slight boat conformation, with Aib C^α and Ile C^α atoms displaced by 0.11 and 0.05 Å on the same side of the mean plane formed by the other ring atoms. Structural features in these two molecules are compared with other related diketopiperazines.     

Specific potassium ion interactions facilitate homocysteine binding to betaine‐homocysteine S‐methyltransferase

Betaine‐homocysteine S‐methyltransferase (BHMT) is a zinc‐dependent methyltransferase that uses betaine as the methyl donor for the remethylation of homocysteine to form methionine. This reaction supports S‐adenosylmethionine biosynthesis, which is required for hundreds of methylation reactions in humans. Herein we report that BHMT is activated by potassium ions with an apparent K_M for K⁺ of about 100 µM. The presence of potassium ions lowers the apparent K_M of the enzyme for homocysteine, but it does not affect the apparent K_M for betaine or the apparent k_cat for either substrate. We employed molecular dynamics (MD) simulations to theoretically predict and protein crystallography to experimentally localize the binding site(s) for potassium ion(s). Simulations predicted that K⁺ ion would interact with residues Asp26 and/or Glu159. Our crystal structure of BHMT bound to homocysteine confirms these sites of interaction and reveals further contacts between K⁺ ion and BHMT residues Gly27, Gln72, Gln247, and Gly298. The potassium binding residues in BHMT partially overlap with the previously identified DGG (Asp26‐Gly27‐Gly28) fingerprint in the Pfam 02574 group of methyltransferases. Subsequent biochemical characterization of several site‐specific BHMT mutants confirmed the results obtained by the MD simulations and crystallographic data. Together, the data herein indicate that the role of potassium ions in BHMT is structural and that potassium ion facilitates the specific binding of homocysteine to the active site of the enzyme. Proteins 2014; 82:2552–2564. © 2014 Wiley Periodicals, Inc.     

The structure of the antitumor complex cis-(diammino) (1,1-cyclobutanedicarboxylato)-Pt(II): X ray and nmr studies

X-ray crystallography shows that Pt(NH₃)₂(CBDCA) is a square-planar complex with the dicarboxylate chelate ring in the boat conformation and a planar cyclobutane ring. ¹H and ¹³C nmr studies suggest that rapid chelate ring flipping occurs in solution. The value of ¹⁹⁵Pt nmr combined with ¹⁵N labeling as an informative new method of studying carboxylate coordination is illustrated. nmr results are also reported for the analogous ethylmalonate complex.     

The STATIS method: Characterization of conformational states of flexible molecules from molecular dynamics simulations in solution

STATIS, a data analysis method used when data can be expressed as matrices, seems particularly well suited to characterize the internal molecular motions and conformational states extracted from MD trajectories. We first outline this method and the “adapted STATIS” method. Applications are presented for 18-crown-6 (simulated for 2 nsec in acetonitrile solution) and for the (L30)₂Cu⁺ catenate (stimulated for 150 psec in chloroform). STATIS should be valuable for the classification of molecular conformations and simplified visualization of MD trajectories.