Aqueous Molecular Dynamics Research Articles

Context Effects on Human Milk Oligosaccharide Linkage Conformation and Dynamics Revealed by MA'AT Analysis.

An emerging NMR method, MA'AT analysis, has been applied to investigate context effects on the conformational properties of several human milk oligosaccharides (HMOs). The MA'AT model of the β-(1→4) linkage in the disaccharide, methyl β-lactoside (MeL), was compared to those obtained for the same linkage in the HMO trisaccharides, methyl 2'-fucosyllactoside (Me2'FL) and methyl 3-fucosyllactoside (Me3FL), and in the tetrasaccharide, methyl 2',3-difucosyllactoside (Me2',3DFL). MA'AT analysis revealed significant context effects on the mean values and circular standard deviations (CSDs) of the psi (ψ) torsion angles in these linkages. α-Fucosylation at both O2'Gal and O3Glc of MeL to give Me2',3DFL significantly constrained librational motion about ψ (70% reduction in the CSD) and shifted its mean value by ∼18°. α-Fucosylation at the O3Glc of MeL to give Me3FL constrained ψ more than α-fucosylation at the O2Gal to give Me2'FL. These effects can be explained by the expected solution conformation of Me3FL, which closely resembles the Lewisx trisaccharide. Comparisons of MA'AT models of ψ to those obtained by 1 μs aqueous molecular dynamics simulation (GLYCAM06) revealed identical trends, that is, MA'AT analysis was able to recapitulate molecular behavior in solution that was heretofore only available from MD simulation. The results highlight the capabilities of MA'AT analysis to determine probability distributions of molecular torsion angles in solution as well as degrees of librational averaging of these angles.

N-Acetyl Side-Chain Conformation in Saccharides: Solution Models Obtained from MA'AT Analysis.

MA'AT analysis has been applied to model the conformational properties of N-acetyl side-chains in biologically important GlcNAc and ManNAc monosaccharides and in a βGlcNAc-(1→4)-βGlcNAc disaccharide. Density functional theory calculations were conducted to obtain parameterized equations that relate the magnitudes and signs of 10 spin-coupling constants to conformations of the C2-N2 bonds of GlcNAc and ManNAc. Six of these equations were used with experimental J-couplings, measured in H2O/2H2O and DMSO-d6 solvents in selectively 13C-labeled compounds, to model the C1-C2-N2-C1' torsion angle (θ1) in GlcNAc and ManNAc residues. MA'AT analysis gave mean values of θ1 of 106° for αGlcNAc and ∼116° for βGlcNAc residues, with circular standard deviations (CSDs) of 21-22° in aqueous solution, in excellent agreement with those obtained by aqueous molecular dynamics (MD) simulation. Parameter space plots revealed unique MA'AT fits of the data, and root mean squared deviations (<0.2 Hz) were twofold smaller than those back-calculated from MD models, indicating that the MA'AT models better fit the experimental J-couplings. Context effects on both θ1 values were found to be small in a βGlcNAc-(1→4)-βGlcNAc disaccharide. MA'AT analysis gave a mean value of θ1 of 249° for αManNAc in H2O/2H2O, with a CSD of ∼19°, with both values in good agreement with MD. MA'AT models of N-acetyl side-chains are similar to those obtained previously for O-acetyl side-chains (J. Phys. Chem. B 2017, 121, 66-77). Both O- and N-acetylation conformationally constrain the C-O or C-N bonds relative to the same bonds in unsubstituted compounds. The present work confirms the ability of MA'AT analysis to reveal relatively small changes in mean molecular torsion angles in solution and provides additional evidence of the method as an experimental tool complementary to MD simulation.

MA'AT Analysis of Aldofuranosyl Rings: Unbiased Modeling of Conformational Equilibria and Dynamics in Solution.

MA'AT analysis has been applied to methyl β-d-ribofuranoside (3) and methyl 2-deoxy-β-d-erythro-pentofuranoside (4) to demonstrate the ability of this new experimental method to determine multi-state conformational equilibria in solution. Density functional theory (DFT) was used to obtain parameterized equations for >20 NMR spin-coupling constants sensitive to furanose ring conformation in 3 and 4, and these equations were used in conjunction with experimental spin-couplings to produce unbiased MA'AT models of ring pseudorotation. These models describe two-state north-south conformational exchange consistent with results obtained from traditional treatments of more limited sets of NMR spin-couplings (e.g., PSEUROT). While PSEUROT, MA'AT, and aqueous molecular dynamics models yielded similar two-state models, MA'AT analysis gives more reliable results since significantly more experimental observables are employed compared to PSEUROT, and no assumptions are needed to render the fitting tractable. MA'AT models indicate a roughly equal distribution of north and south ring conformers of 4 in aqueous (2H2O) solution compared to ∼80% north forms for 3. Librational motion about the mean pseudorotation phase angles P of the preferred north and south conformers of 3 in solution is more constrained than that for 4. The greater rigidity of the β-ribo ring may be caused by synergistic stereoelectronic effects and/or noncovalent (e.g., hydrogen-bonding) interactions in solution that preferentially stabilize north forms of 3. MA'AT analysis of oligonucleotides and other furanose ring-containing biomolecules promises to improve current experimental models of sugar ring behavior in solution and help reveal context effects on ring conformation in more complex biologically important systems.

Aqueous Molecular Dynamics for Understanding Glycosaminoglycan Recognition by Proteins.

Glycosaminoglycans (GAGs) are biopolymers that exist in most organisms. GAGs are known to bind to hundreds of proteins and partake in multiple biological processes such as growth, morphogenesis, inflammation, infection, and others. Their intrinsic structural heterogeneity and conformational variability introduce major challenges in experimental studies. On the other hand, recent advances in force field development and computational technology have yielded phenomenal opportunity to study thousands of GAG sequences simultaneously to understand recognition of target protein(s). Here, we describe experimental setup for conventional molecular dynamics simulations of GAGs to position an experimental biologist favorably in performance, analysis and interpretation of stability, specificity, and conformational properties of GAGs, while also elucidating their interactions with amino acid residues of a protein at an atomistic level in presence of water.

How the water-soluble hemicarcerand incarcerates guests at room temperature decoded with modular simulations

Molecular dynamics simulations of hemicarcerands and related variants allow the study of constrictive binding and offer insight into the rules of molecular complexation, but are limited because three-dimensional models of hemicarcerands are tedious to build and their atomic charges are complicated to derive. There have been no molecular dynamics simulations of the reported water-soluble hemicarcerand (Octacid4) that explain how Octacid4 encapsulates guests at 298 K and keeps them encapsulated at 298 K in NMR experiments. Herein we report a modular approach to hemicarcerand simulations that simplifies the model building and charge derivation in a manner reminiscent of the approach to protein simulations with truncated amino acids as building blocks. We also report that in aqueous molecular dynamics simulations at 298 K apo Octacid4 adopts two clusters of conformations one of which has an equatorial portal open but the guest-bound Octacid4 adopts one cluster of conformations with all portals closed. These results explain how Octacid4 incarcerates guests at room temperature and suggest that the guest-induced host conformational change that impedes decomplexation is a previously unrecognized conformational characteristic that promotes strong molecular complexation.

Use of Circular Statistics To Model αMan-(1→2)-αMan and αMan-(1→3)-α/βMan O-Glycosidic Linkage Conformation in 13C-Labeled Disaccharides and High-Mannose Oligosaccharides.

A new experimental method, MA' AT analysis, has been applied to investigate the conformational properties of O-glycosidic linkages in several biologically important mannose-containing di- and oligosaccharides. Methyl α-d-mannopyranosyl-(1→2)-α-d-mannopyranoside (2), methyl α-d-mannopyranosyl-(1→3)-α-d-mannopyranoside (3), and methyl α-d-mannopyranosyl-(1→3)-β-d-mannopyranoside (4) were prepared with selective 13C-enrichment to enable the measurement of NMR scalar couplings across their internal O-glycosidic linkages. Density functional theory (DFT) was used to parameterize equations for JCH and JCC values in 2-4 that are sensitive to phi (ϕ) and psi (ψ). The experimental J-couplings and parameterized equations were treated using a circular statistics algorithm encoded in the MA' AT program. Conformations about ϕ and ψ treated using single-state von Mises models gave excellent fits to the ensembles of redundant J-couplings. Mean values and circular standard deviations (CSDs) for each linkage torsion angle ϕ (CSD) and ψ (CSD) in 2, -29° (25°) and 20° (22°); in 3, -36° (36°) and 8° (27°); in 4, -37° (34°) and 10° (26°); ϕ = H1'-C1'-O1'-CX and ψ = C1'-O1'-CX-HX (CX = aglycone carbon) were compared to histograms obtained from 1 μs aqueous molecular dynamics (MD) simulations and X-ray database statistical analysis. MA' AT-derived models of ψ were in very good agreement with the MD and X-ray data, but not those of ϕ, suggesting a need for force field revision. The effect of structural context on linkage conformation was also investigated in four selectively 13C-labeled homomannose tri- and tetrasaccharides using the MA' AT method. In the cases examined, context effects were found to be small.

Conformational Populations of β-(1→4) O-Glycosidic Linkages Using Redundant NMR J-Couplings and Circular Statistics.

Twelve disaccharides containing β-(1→4) linkages and displaying systematic structural variations in the vicinity of these linkages were selectively labeled with 13C to facilitate measurements of multiple NMR spin-spin (scalar; J) coupling constants (JCH and JCC values) across their O-glycosidic linkages. Ensembles of spin-couplings (2JCOC, 3JCOCH, 3JCOCC) sensitive to the two linkage torsion angles, phi (ϕ) and psi (ψ), were analyzed by using parametrized equations obtained from density functional theory (DFT) calculations, Fredholm theory, and circular statistics to calculate experiment-based rotamer populations for ϕ and ψ in each disaccharide. With the statistical program MA'AT, torsion angles ϕ and ψ were modeled as a single von Mises distribution, which yielded two parameters, the mean position and the circular standard deviation (CSD) for each angle. The NMR-derived rotamer populations were compared to those obtained from 1 μs aqueous molecular dynamics (MD) simulations and crystallographic database statistical analyses. Conformer populations obtained exclusively from the MA'AT treatment of redundant J-couplings were in very good agreement with those obtained from the MD simulations, providing evidence that conformational populations can be determined by NMR for mobile molecular elements such as O-glycosidic linkages with minimal input from theory. The approach also provides an experimental means to validate the conformational preferences predicted from MD simulations. The conformational behaviors of ϕ in the 12 disaccharides were very similar, but those of ψ varied significantly, allowing a classification of the 12 disaccharides based on preferred linkage conformation in solution.

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors.

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (ELJ) and electrostatic (EC) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pKi = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R2 = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.

Tacticity effects on conformational structure and hydration of poly-(methacrylic acid) in aqueous solutions-a molecular dynamics simulation study

The influence of tacticity on chain dimensions, backbone and side-group conformational states and relaxation dynamics, intermolecular hydrogen bonding and its relaxation dynamics was investigated for 30 repeat unit poly(methacrylic acid) (PMA) as a function of the degree-of-neutralisation, f (i.e. charge density) [0 < f < 1 range] in explicit water with Na+ neutralising counter-ions, using molecular dynamics simulations. Chain expansion with increase in charge density is observed. Simulation results, where applicable, compare well with experimental results in the literature. Isotactic (i-PMA) chain exhibits the largest expansion (89% change in Rg), and syndiotactic (s-PMA-RR-0.7) chain shows 31%. For fully neutralised chain (high charge density), the probability for trans conformation at backbone bonds follows the trend i-PMA > a-PMA > s-PMA. At high charge density, a higher probability of trans state is obtained at meso dyads as compared to racemic dyads and the side group in i-PMA relative to other tacticity is conformationally more extended. RR triads impart better hydrogen bonding with water. In the COOH side group, greater rotational mobility is observed for i-PMA as compared to s-PMA. Water coordination to PMA is invariant with tacticity at low f. For f = 1, s-PMA exhibits the best level of H-bonding, due to its 3D chemical structural configuration. Binding distance between water and PMA atoms and their coordination number remain unchanged with respect to ionisation of COOH groups. Isotactic (i-PMA) exhibits the fastest relaxation of polymer-water H-bonds. While counter-ion coordination on to PMA is unaffected by tacticity, the intensity of ion condensation as well as hydrogen bond relaxation time at f = 1 varies in the order i-PMA > a-PMA > s-PMA.

Comparative simulation of pneumococcal serogroup 19 polysaccharide repeating units with two carbohydrate force fields

Streptococcus pneumoniae causes meningitis, pneumonia and severe invasive disease (IPD) in young children. Although widespread infant immunisation with the PCV7 seven-valent pneumococcal conjugate vaccine has led to a dramatic decrease in IPD, infections due to non-vaccine serotypes, particularly serotype 19A, have increased. As the 19F polysaccharide differs from 19A at a single linkage position, it was assumed that PCV7 (containing 19F) would cross-protect against 19A disease. However, vaccination with PCV7 results in only 26% effectiveness against IPD caused by 19A. We explored the conformations and dynamics of the polysaccharide repeating units from serotypes 19F and 19A, comparing free energy surfaces for glycosidic linkages with 100ns aqueous molecular dynamics simulations of the di- and trisaccharide components. All calculations were performed with both the CHARMM and the GLYCAM carbohydrate force fields to establish whether the choice of model affects the predicted molecular behaviour. Although we identified key differences between the force fields, overall they were in agreement in predicting a 19F repeating unit with a wider range of conformation families than the more restricted 19A trisaccharide. This suggests a probable conformational difference between the 19F and 19A polysaccharides, which may explain the low cross-protection of 19F vaccines against 19A disease.

Shaping up for structural glycomics: a predictive protocol for oligosaccharide conformational analysis applied to N-linked glycans

The human glycome comprises a vast untapped repository of 3D-structural information that holds the key to glycan recognition and a new era of rationally designed mimetic chemical probes, drugs, and biomaterials. Toward routine prediction of oligosaccharide conformational populations and exchange rates at thermodynamic equilibrium, we apply hardware-accelerated aqueous molecular dynamics to model μs motions in N-glycans that underpin inflammation and immunity. In 10μs simulations, conformational equilibria of mannosyl cores, sialyl Lewis (sLe) antennae, and constituent sub-sequences agreed with prior refinements (X-ray and NMR). Glycosidic linkage and pyranose ring flexing were affected by branching, linkage position, and secondary structure, implicating sequence dependent motions in glycomic functional diversity. Linkage and ring conformational transitions that have eluded precise quantification by experiment and conventional (ns) simulations were predicted to occur on μs timescales. All rings populated non-chair shapes and the stacked galactose and fucose pyranoses of sLea and sLex were rigidified, suggesting an exploitable 3D-signature of cell adhesion protein binding. Analyses of sLex dynamics over 25μs revealed that only 10μs were sufficient to explore all aqueous conformers. This simulation protocol, which yields conformational ensembles that are independent of initial 3D-structure, is proposed as a route to understanding oligosaccharide recognition and structure–activity relationships, toward development of carbohydrate-based novel chemical entities.

Conformational Analysis of Trimannoside and Bisected Trimannoside Using Aqueous Molecular Dynamics Simulations

The conformational properties of oligosaccharides are important to understand carbohydrate-protein interactions. A trimannoside, methyl 3,6-di-O-(<TEX>$\alpha$</TES>-D-Man)-<TEX>$\alpha$</TES>-D-Man (TRIMAN) is a basic unit of N-linked oligosaccharides. This TRIMAN moiety was further modified by GlcNAc (BISECT), which is important to biological activity of N-glycan. To characterize the trimannoside and its bisecting one we performed a molecular dynamics simulation in water. The resulting models show the conformational transition with two major and minor conformations. The major conformational transition results from the <TEX>$\omega$</TEX> angle transition; another minor transition is due to the <TEX>$\psi$</TEX> angle transition of <TEX>$\alpha$</TES> (1 <TEX>$\rightarrow$</TEX> 6) linkage. The introduction of bisecting GlcNAc on TRIMAN made the different population of the major and minor conformations of the TRIMAN moiety. Omega (<TEX>$\omega$</TEX>) angle distribution is largely changed and the population of gt conformation is increased in BISECT oligosaccharide. The inter-residue hydrogen bonds and water bridges via bisecting GlcNAc residue make alterations on the local and overall conformation of TRIMAN moiety. These changes of conformational distribution for TRIMAN moiety can affect the overall conformation of N-glycan and the biological activity of glycoprotein.

Electron Donation in the Water–Water Hydrogen Bond

Hydrogen bonding is central to aqueous systems, from small water clusters and microsolvated ions to bulk water, and solvated biomolecules. Despite numerous experimental and theoretical studies, the physical nature of hydrogen bonding is still debated. One issue is the extent of intermolecular charge transfer (CT) in hydrogen bonding. Natural bond orbital (NBO) analysis and natural energy decomposition analysis suggest that CT is predominant because, if CT is neglected, NBO analysis shows no binding at the water-dimer equilibrium geometry. However, other earlier decomposition methods estimated that CT contributes only around 20% of the overall binding energy. This question has practical significance for aqueous molecular dynamics simulations, where models based on purely electrostatic potentials, such as Coulomb plus Lennard-Jones with perhaps polarizability, seem to be very successful in reproducing many structural and thermodynamic properties of water. Is such good agreement soundly based or fortuitous? Recent X-ray absorption and X-ray Raman scattering experiments have challenged the accepted locally tetrahedral structure of liquid water. The failure of classical molecular dynamics simulations to reproduce the “chain and ring” local structure inferred from these experiments has generated questions about the reliability of existing water potentials. 12] This fact, combined with the CT character of the hydrogen bonding suggested by NBO analysis, has led to proposals to incorporate CT effects into empirical water potentials. However, the “chain and ring” interpretation of the X-ray experiments is highly controversial and has been challenged on many fronts. Herein, we rationalize the success of empirical electrostatic potentials by clarifying intermolecular CT effects in the simplest water cluster, the water dimer. We have used our recently developed energy-decomposition analysis (EDA) and charge-transfer analysis (CTA) based on absolutely localized molecular orbitals (ALMOs), which are ideal for separating CT from frozen density and polarization interactions. In ALMO EDA, the frozen density (FRZ) term is calculated as the interaction energy of the unrelaxed electron densities on the molecules. The polarization (POL) term is due to the deformation (or polarization) of the electron clouds of the molecules in the field of each other. Quantum mechanically, it is described as the energy lowering due to the intramolecular relaxation of each molecule s ALMOs in the field of the other molecule. CT is calculated as the energy lowering due to the intermolecular relaxation of the molecular orbitals corrected for the basis set superposition error (BSSE). Like related earlier methods, the ALMO atomic-tomolecular orbital transformation is constrained to be blockdiagonal in terms of the molecular fragments (prohibiting CT). Unlike those earlier methods, ALMO EDA and CTA treat the optimization of the ALMO s in a variationally optimal way. CT effects (energy lowering due to electron transfer from occupied orbitals on one molecule to virtual orbitals of another molecule, and then any further repolarization or higher-order relaxation) are corrections to the optimal polarized reference system, and cannot be over or underestimated. The ALMO charge transfer scale, DQ, provides a measure of the distortion of the electronic clouds upon formation of an intermolecular bond and is such that all CT terms, that is, forward-donation, back-donation, and higher order relaxation, have well defined energetic effects. The water dimer geometry with Cs symmetry was optimized at the MP2/aug-cc-pVQZ level. All calculations were [a] Dr. R. Z. Khaliullin Department of Chemistry, University of California Berkeley Berkeley, CA 94720 (USA) Fax: (+1)510-643-1255 E-mail : rustam@khaliullin.com [b] Prof. A. T. Bell Department of Chemical Engineering, University of California Berkeley Berkeley, CA 94720 (USA) [c] Prof. M. Head-Gordon Department of Chemistry, University of California Berkeley Berkeley, CA 94720 (USA) Fax: (+1)510-643-1255 E-mail : mhg@bastille.cchem.berkeley.edu Supporting information for this article is available on the WWW under http://dx.doi.org/10.1002/chem.200802107.

Comparison of Aqueous Molecular Dynamics with NMR Relaxation and Residual Dipolar Couplings Favors Internal Motion in a Mannose Oligosaccharide

An investigation has been performed to assess how aqueous dynamical simulations of flexible molecules can be compared against NMR data. The methodology compares state-of-the-art NMR data (residual dipolar coupling, NOESY, and (13)C relaxation) to molecular dynamics simulations in water over several nanoseconds. In contrast to many previous applications of residual dipolar coupling in structure investigations of biomolecules, the approach described here uses molecular dynamics simulations to provide a dynamic representation of the molecule. A mannose pentasaccharide, alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->2)-D-Manp, was chosen as the model compound for this study. The presence of alpha-linked mannan is common to many glycopeptides, and therefore an understanding of the structure and the dynamics of this molecule is of both chemical and biological importance. This paper sets out to address the following questions. (1) Are the single structures which have been used to interpret residual dipolar couplings a useful representation of this molecule? (2) If dynamic flexibility is included in a representation of the molecule, can relaxation and residual dipolar coupling data then be simultaneously satisfied? (3) Do aqueous molecular dynamics simulations provide a reasonable representation of the dynamics present in the molecule and its interaction with water? In summary, two aqueous molecular dynamics simulations, each of 20 ns, were computed. They were started from two distant conformations and both converged to one flexible ensemble. The measured residual dipolar couplings were in agreement with predictions made by averaging the whole ensemble and from a specific single structure selected from the ensemble. However, the inclusion of internal motion was necessary to rationalize the relaxation data. Therefore, it is proposed that although residual dipolar couplings can be interpreted as a single-structure, this may not be a correct interpretation of molecular conformation in light of other experimental data. Second, the methodology described here shows that the ensembles from aqueous molecular dynamics can be effectively tested against experimental data sets. In the simulation, significant conformational motion was observed at each of the linkages, and no evidence for intramolecular hydrogen bonds at either alpha(1-->2) or alpha(1-->3) linkages was found. This is in contrast to simulations of other linkages, such as beta(1-->4), which are often predicted to maintain intramolecular hydrogen bonds and are coincidentally predicted to have less conformational freedom in solution.

Quantitative conformational analysis of the core region of N-glycans using residual dipolar couplings, aqueous molecular dynamics, and steric alignment.

A method is described for quantitatively investigating the dynamic conformation of small oligosaccharides containing an alpha(1 --> 6) linkage. It was applied to the oligosaccharide Man-alpha(1 --> 3) [Man-alpha(1 --> 6)] Man-alpha-O-Me, which is a core region frequently observed in N-linked glycans. The approach tests an aqueous molecular dynamics simulation, capable of predicting microscopic dynamics, against experimental residual dipolar couplings, by assuming that alignment is caused purely by steric hindrance. The experimental constraints were heteronuclear and homonuclear residual dipolar couplings, and in particular those within the alpha(1 --> 6) linkage itself. Powerful spin-state-selective pulse sequences and editing schemes were used to obtain the most relevant couplings for testing the model. Molecular dynamics simulations in water over a period of 50 ns were not able to predict the correct rotamer population at the alpha(1 --> 6) linkage to agree with the experimental data. However, this sampling problem could be corrected using a simple maximum likelihood optimisation, indicating that the simulation was modelling local dynamics correctly. The maximum likelihood prediction of the residual dipolar couplings was found to be an almost equal population of the gg and gt rotamer conformations at the alpha(1--> 6) linkage, and the tg conformation was predicted to be unstable and unpopulated in aqueous solution. In this case all twelve measured residual dipolar couplings could be satisfied. This conformer population could also be used to make predictions of scalar couplings with the use of a previously derived empirical equation, and is qualitatively in agreement with previous predictions based on NMR, X-ray crystallography and optical data.

Glycosaminoglycan conformation: do aqueous molecular dynamics simulations agree with x-ray fiber diffraction?

Glycosaminoglycan-protein interactions are biologically important and require an appreciation of glycan molecular shape in solution, which is presently unavailable. In previous studies we found strong similarity between aqueous molecular dynamics (MD) simulations and published x-ray diffraction refinements of hyaluronan. We have applied a similar approach here to chondroitin and dermatan, attempting to clarify some of the issues raised by the x-ray diffraction literature relating to chondroitin and dermatan sulfate. We predict that chondroitin has the same beta(1-->4) linkage conformation as hyaluronan, and that their average beta(1-->3) conformations differ. This is explained by changes in hydrogen-bonding across this linkage, resulting from its axial hydroxyl, causing a different sampling of left-handed helices in chondroitin (2.5- to 3.5-fold) as compared with hyaluronan (3.0- to 4.0-fold). Few right-handed helices, which lack intramolecular hydrogen-bonds, were sampled during our MD simulations. Thus, we propose that the 8-fold helix observed in chondroitin-6-sulfate, represented in the literature as an 8(3) helix (right-handed), though it has never been refined, is more likely to be 8(5) (left-handed) helix. Molecular dynamics simulations implied that (4)C(1) and (2)S(O), but not (1)C(4), forms of iduronate could be used in refinements of dermatan x-ray fiber diffraction patterns. Current models of 8-fold dermatan sulfate chains containing (4)C(1) iduronate refine to right-handed helices, which possess no intramolecular hydrogen-bonds. However, MD simulations predict that models containing (2)S(O) iduronate could provide better (8(5) helix) starting structures for refinement. Thus, the 8-fold dermatan sulfate refinement (8(3) helix) could be in error.

Restrained molecular dynamics of solvated duplex DNA using the particle mesh Ewald method.

Restrained and unrestrained aqueous solution molecular dynamics simulations applying the particle mesh Ewald (PME) method to DNA duplex structures previously determined via in vacuo restrained molecular dynamics with NMR-derived restraints are reported. Without experimental restraints, the DNA decamer, d(CATTTGCATC).d(GATGCAAATG) and trisdecamer, d(AGCTTGCCTTGAG).d(CTCAAGGCAAGCT), structures are stable on the nanosecond time scale and adopt conformations in the B-DNA family. These free DNA simulations exhibit behavior characteristic of PME simulations previously performed on DNA sequences, including a low helical twist, frequent sugar pucker transitions, BI-BII(epsilon-zeta) transitions and coupled crakshaft (alpha-gamma) motion. Refinement protocols similar to the original in vacuo restrained molecular dynamics (RMD) refinements but in aqueous solution using the Cornell et al. force field [Cornell et al. (1995) J. Am. Chem. Soc., 117, 5179-5197] and a particle mesh Ewald treatment produce structures which fit the restraints very well and are very similar to the original in vacuo NMR structure, except for a significant difference in the average helical twist. Figures of merit for the average structure found in the RMD PME decamer simulations in solution are equivalent to the original in vacuo NMR structure while the figures of merit for the free MD simulations are significantly higher. The free MD simulations with the PME method, however, lead to some sequence-dependent structural features in common with the NMR structures, unlike free MD calculations with earlier force fields and protocols. There is some suggestion that the improved handling of electrostatics by PME improves long-range structural aspects which are not well defined by the short-range nature of NMR restraints.

Deducing polymeric structure from aqueous molecular dynamics simulations of oligosaccharides: predictions from simulations of hyaluronan tetrasaccharides compared with hydrodynamic and X-ray fibre diffraction data

Molecular dynamics simulations of the two hyaluronan tetrasaccharides in water predict that over a period of 500 ps, their central linkages populate a single primary minima. Over the same period the peripheral linkages explore this minima, but also a secondary minima. Structures constructed using the primary minima were found to be extended left-handed helices of axial rise per disaccharide (h) 0.8 to 1.0 nm and 2.8 to 4.5 disaccharides per turn (n), in good agreement with n=3 and n=4 helices found by X-ray fibre diffraction studies. We have used the predicted average conformation from molecular dynamics to calculate the translational diffusion coefficients of the oligosaccharide series up to decasaccharide, and compared these with experimental measurements obtained using the method of capillary dispersion. Our calculated values are found to be in good agreement with experiment beyond the size of a tetrasaccharide. A partial digest of hyaluronan in the molecular mass range 10 to 100 kDa was fractionated by gel chromatography. Molecular weights were determined by in-line laser light-scattering measurements, and the translational diffusion coefficients of selected fractions were determined by dynamic laser light-scattering. A similar experiment was performed on hyaluronan with a molecular mass greater than 1MDa. The data suggest a change from rod-like to stiff coil behaviour beyond a molecular weight of 10 kDa. We have also examined the conformations available using the secondary minima, found at the peripheral linkages. In contrast to the extended structures previously described we have found left and right-handed helices with high values of n (5–10) and low values of h. Although there is no experimental evidence for these structures, they are of interest as, over short stretches, they would introduce folds, loops, and turns into the hyaluronan molecule. Such shapes may play an important role in the hydrodynamics of hyaluronan and its interaction with lipids and proteins.

Effect of periodic box size on aqueous molecular dynamics simulation of a DNA dodecamer with particle-mesh Ewald method

The particle-mesh Ewald (PME) method is considered to be both efficient and accurate for the evaluation of long-range electrostatic interactions in large macromolecular systems being studied by molecular dynamics simulations. This method assumes "infinite" periodic boundary conditions resembling the symmetry of a crystal environment. Can such a "solid-state" method accurately portray a macromolecular solute such as DNA in solution? To address this issue, we have performed three 1500-ps PME molecular dynamics (MD) simulations, each with a different box size, on the d(CGCGA6CG)-(CGT6CGCG) DNA dodecamer. The smallest box had the DNA solvated by a layer of water molecules of at least 5 A along each orthogonal direction. The intermediate size box and the largest box had the DNA solvated by a layer of water molecules of at least 10 A and 15 A, respectively, along each orthogonal direction. The intermediate size box in the present study is similar to the box size currently chosen by most workers in the field. Based on a comparison of RMSDs and curvature for this single DNA dodecamer sequence, the larger two box sizes do not appear to afford any extra benefit over the smallest box. The implications of this finding are briefly discussed.

Non-peptidic inhibitors of human leukocyte elastase. 4. Design, synthesis, and in vitro and in vivo activity of a series of beta-carbolinone-containing trifluoromethyl ketones.

A novel series of human leukocyte elastase (HLE) inhibitors containing the beta-carbolinone ring system are reported. The design of these trifluoromethyl ketone-based inhibitors used a combination of structural information obtained from X-ray crystallography and molecular modeling investigations. The beta-carbolinone ring in these compounds serves as a highly efficient peptidiomimetic for the P2-P3 region of peptidyl trifluoromethyl ketone inhibitors of HLE. Several of the beta-carbolinones exhibit significant in vitro potency, with Ki values in the nanomolar range. Using aqueous molecular dynamics simulations, realistic models for the molecular recognition of these inhibitors by HLE have been obtained and are discussed. This series of compounds are found to have excellent selectivity for HLE over a number of other proteolytic enzymes, including closely related enzymes such as porcine pancreatic elastase.