Conformation of the oligosaccharide receptor for E-selectin

A tetrasaccharide related to the blood group oligosaccharides, known as sialyl Lewis^X, has been proposed as the receptor for the lectin responsible for leukocyte adhesion named alternatively as E-selectin or ELAM-1. The ¹³C- and ¹H-nmr spectra have been completely assigned for a tetrasaccharide model of this receptor, Neu5Ac α-(2 → 3)-Gal β-(1 → 4)-[ Fuc α-(1 → 3)-] GlcNAc β-NHAc. Quantitative nuclear Overhauser data (NOESY) have been recorded and analyzed by a complete spin matrix simulation method. Conformational space was exhaustively searched and all conformational models whose simulated NOESY spectra matched the experiment were found. Molecular mechanics and molecular dynamics calculations were carried out to test whether the experimental conformations are low energy and thus likely to represent true single conformations for the tetrasaccharide. It was concluded that while the Lewis^X trisaccharide portion of the compound adopts a single conformation, there is likely to be some flexibility about the Neu5Ac α-(2 → 3)-linkage. A model featuring fast exchange between two different conformations of this linkage is found to be consistent with both the nmr experiments and the molecular dynamics simulations. © 1994 John Wiley & Sons, Inc.     

A Combined 2D-NMR and Molecular Dynamics Analysis of the Structure of the Actinomycin D: d(ATGCAT)2 Complex

Abstract

We present a comparative analysis of an NMR experiment and molecular and harmonic dynamics simulations of an actinomycin D: d(ATGCAT)₂ complex. A comparison of NOE measurements and 1/R⁶ weighted proton-proton distances confirm the general correctness of the Actinomycin D-DNA model proposed by Sobell. There are, however, some substantial differences between the proton-proton distances inferred from the NOE results and the molecular and harmonic dynamics simulations. The remaining discrepancies could either come from contributions of other conformations to the average properties of the complex or from uncertainties in the NMR distance analysis. An analysis of the molecular dynamics helix properties, sugar puckers, hydrogen bonding, rms fluctuations and torsional properties are qualitatively consistent with those from previous simulations, but the presence of an intercalated drug leads to some new structural and dynamical features.     

Probing the helical forms of Ca2+‐guluronan junction zones in alginate gels by molecular dynamics 1: Duplexes

A molecular dynamics investigation of the helical forms adopted by (1→4)‐α‐L ‐guluronan in explicit water environment was carried out. Single chains and duplexes were modeled at 300 K starting both from 2₁ or 3₂ helical conformations and in the presence of a neutralizing amount of Ca²⁺ ions. All systems were allowed full conformational freedom. The initial perfect helices with integral screw symmetries were lost at the very beginning of simulations and two distinct behaviors were observed: At equilibrium the 2₁ models mostly retained the 2₁ local helical conformations while exploring the 3₂ ones the rest of the time. In duplexes the two chains, which behaved similarly, were well extended and slightly twisted. By contrast, the chains in 3₂ duplex models were dissimilar and explored a much broader conformational space in which 2₁ and 3₂ local helical conformations were dominant and equally represented but the 3₁ and other conformations were also present. The wide variety of conformations revealed in this study is consistent with the general difficulty in obtaining crystals of Ca²⁺‐guluronate with suitable lateral dimensions for crystallographic studies. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 562–571, 2013.     

The binding properties of the H5N1 influenza virus neuraminidase as inferred from molecular modeling

The avian influenza H5N1 virus has emerged as an important pathogen, causing severe disease in humans and posing a pandemic threat. Substrate specificity is crucial for the virus to obtain the ability to spread from avian to human. Therefore, an investigation of the binding properties of ligands at the molecular level is important for understanding the catalytic mechanism of the avian influenza virus neuraminidase and for designing novel and specific inhibitors of H5N1 neuraminidase. Based on the available crystal structure of H5N1, we have characterized the binding properties between sialic acid, methyl 3’sialyllactoside, methyl 6’sialyllactoside and the H5N1 influenza virus neuraminidase using molecular docking and molecular dynamics simulations. Obtained molecular dynamics trajectories were analyzed in terms of ligand conformations, N1-ligand interactions, and in terms of loop flexibility. It was found that in the N1-SA complex the sialic acid ring undergoes a transition from the B _2,5 to the ² C ₅ conformation. However, in the N1-3SL and N1-6SL complexes sialic acid remained in the distorted boat conformation. The obtained results indicate that 3SL has only weak interactions with the 150-loop, whereas the N1-6SL complex shows strong interactions. Most of the differences arise from the various conformations around the glycosidic linkage, between the sialic acid and galactose, which facilitate the above interactions of 6SL with the enzyme, and as a consequence the interactions between the 150- and 430- loops. This finding suggests that the altered flexibility of loops in and around the active site is one of the reasons why the avian N1 preferentially cleaves sialic acid from α-(2-3)-Gal glycoconjugates over α-(2-6)-Gal. These molecular modeling results are consistent with available experimental results on the specificity of N1.     

Conformational transitions of a dipeptide in water: Effects of imposed pathways using umbrella sampling techniques

The free energy difference between two states of a molecular system separated by an energy barrier can generally be computed using the technique of umbrella sampling along a chosen reaction coordinate or pathway. The effect of a particular choice of pathway upon the obtained free energy difference is investigated by molecular dynamics simulation of a model system consisting of a glycine dipeptide in aqueous solution. Two different reaction coordinates connecting the so-called C₅ and C₇ conformations, one involving intramolecular hydrogen bonds and the other involving the peptide ?, ψ angles, are considered. The Gibbs free energy differences ΔG(C₅ – C₇) are small in both cases, 1.5 ± 1 kJ mol^?1 and 2.2 ± 1 kJ mol ^?1, respectively. The two different reaction coordinates yield free energy differences that are identical to within their statistical error. It is found that the exchange of solute–solute, solute–water, and water–water hydrogen bonds involves free energy changes of less than k_BT, which points at the existence of a multitutde of low free energy pathways connecting the C₅ and C₇ dipeptide conformations. © 1994 John Wiley & Sons, Inc.     

Identification of Distinct Conformations of the Angiotensin-II Type 1 Receptor Associated with the Gq/11 Protein Pathway and the β-Arrestin Pathway Using Molecular Dynamics Simulations

Biased signaling represents the ability of G protein-coupled receptors to engage distinct pathways with various efficacies depending on the ligand used or on mutations in the receptor. The angiotensin-II type 1 (AT₁) receptor, a prototypical class A G protein-coupled receptor, can activate various effectors upon stimulation with the endogenous ligand angiotensin-II (AngII), including the G_q/11 protein and β-arrestins. It is believed that the activation of those two pathways can be associated with distinct conformations of the AT₁ receptor. To verify this hypothesis, microseconds of molecular dynamics simulations were computed to explore the conformational landscape sampled by the WT-AT₁ receptor, the N111G-AT₁ receptor (constitutively active and biased for the G_q/11 pathway), and the D74N-AT₁ receptor (biased for the β-arrestin1 and -2 pathways) in their apo-forms and in complex with AngII. The molecular dynamics simulations of the AngII-WT-AT₁, N111G-AT₁, and AngII-N111G-AT₁ receptors revealed specific structural rearrangements compared with the initial and ground state of the receptor. Simulations of the D74N-AT₁ receptor revealed that the mutation stabilizes the receptor in the initial ground state. The presence of AngII further stabilized the ground state of the D74N-AT₁ receptor. The biased agonist [Sar¹,Ile⁸]AngII also showed a preference for the ground state of the WT-AT₁ receptor compared with AngII. These results suggest that activation of the G_q/11 pathway is associated with a specific conformational transition stabilized by the agonist, whereas the activation of the β-arrestin pathway is linked to the stabilization of the ground state of the receptor.     

Structure and hydrogen bonds of cyclohexapeptide RA-VII by molecular dynamics simulations and quantum chemical calculations

We computationally examined the structure of anti-tumour bicyclic hexapeptide RA-VII. This peptide adopts three conformations (confs.), A, B and C, in dimethyl sulfoxide (DMSO). Although it was experimentally reported that the structure of conf. A is important for anti-tumour activity, the dynamics of confs. A, B and C are not well known. We performed quantum chemical calculations and molecular dynamics (MD) simulations of RA-VII in DMSO. The MD simulations indicated two different local stable structures for conf. C: a structure containing a bent 18-membered ring and another structure containing a rotated peptide bond between Tyr⁶ and d-Ala¹. The root-mean-square fluctuation of the 14-membered ring for conf. A was larger than that for confs. B and C. Ala⁴ formed intramolecular hydrogen bonds more often in conf. A than in the other conformations. A large number of hydrogen bonds and large structural fluctuations are important for the anti-tumour activity of RA-VII. Our results for the structural change of conf. C and the analysis of the dynamics for confs. A, B and C may contribute to the design of new analogues of cyclic peptides.     

Assessing Structural Determinants of Zn2+ Binding to Human HV1 via Multiple MD Simulations

Voltage-gated proton channels (H_V1) are essential for various physiological tasks but are strongly inhibited by Zn²⁺ cations. Some determinants of Zn²⁺ binding have been elucidated experimentally and in computational studies. However, the results have always been interpreted under the assumption that Zn²⁺ binds to monomeric H_V1 despite evidence that H_V1 expresses as a dimer and that the dimer has a higher affinity for zinc than the monomer and experimental data that suggest coordination in the dimer interface. The results of former studies are also controversial, e.g., supporting either one single or two binding sites. Some structural determinants of the binding are still elusive. We performed a series of molecular dynamics simulations to address different structures of the human proton channel, the monomer and two plausible dimer conformations, to compare their respective potential to interact with and bind Zn²⁺ via the essential histidines. The series consisted of several copies of the system to generate independent trajectories and increase the significance compared to a single simulation. The amount of time simulated totals 29.9 μs for 126 simulations of systems comprising ∼59,000 to ∼187,000 atoms. Our approach confirms the existence of two binding sites in monomeric and dimeric human H_V1. The dimer interface is more efficient for attracting and binding Zn²⁺ via the essential histidines than the monomer or a dimer with the histidines in the periphery. The higher affinity is due to the residues in the dimer interface that create an attractive electrostatic potential funneling the zinc cations toward the binding sites.     

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.     

Minimum-Energy Path for a U6 RNA Conformational Change Involving Protonation, Base-Pair Rearrangement and Base Flipping

The U6 RNA internal stem-loop (U6 ISL) is a highly conserved domain of the spliceosome that is important for pre-mRNA splicing. The U6 ISL contains an internal loop that is in equilibrium between two conformations controlled by the protonation state of an adenine (pK_a = 6.5). Lower pH favors formation of a protonated C-A⁺ wobble pair and base flipping of the adjacent uracil. Higher pH favors stacking of the uracil and allows an essential metal ion to bind at this position. Here, we define the minimal-energy path for this conformational transition. To do this, we solved the U6 ISL structure at higher pH (8.0) in order to eliminate interference from the low-pH conformer. This structure reveals disruption of the protonated C-A⁺ pair and formation of a new C-U pair, which explains the preference for a stacked uracil at higher pH. Next, we used nudged elastic band molecular dynamics simulations to calculate the minimum-energy path between the two conformations. Our results indicate that the C-U pair is dynamic, which allows formation of the more stable C-A⁺ pair upon adenine protonation. After formation of the C-A⁺ pair, the unpaired uracil follows a minor-groove base-flipping pathway. Molecular dynamics simulations suggest that the extrahelical uracil is stabilized by contacts with the adjacent helix.     

Conformational analysis and molecular dynamics simulation of α-(1 → 2) and α-(1 → 3) linked rhamnose oligosaccharides: Reconciliation with optical rotation and NMR experiments

Molecular mechanics and dynamics calculations were carried out on the disaccharides α-L-Rhap-(1 → 2)-α-L-Rhap-(1 → OMe) (1) and α-L-Rhap-(1 → 3)-α-L-Rhap-(1 OMe) (2), and the trisaccharide α-L-Rhap-(1 → 2)-α-L-Rhap-(1 → 3)-α-L-Rhap-(1 → OMe) (3). The semiflexible conformational behavior of these molecules was characterized by the occupation of a combination of different glycosidic linkage and side-chain conformational positions whose relative occupations were sensitive to dielectric screening. Molecular dynamics simulations of the trisaccharide 3 showed little difference between the linkage conformations in the trisaccharide and the component disaccharides 1 and 2. Experimental optical rotation data of 1 and 2 were obtained as a function of temperature in varying solvents. The molecular models were combined with the semiempirical theory of Stevens and Sathyanarayana to yield calculated optical rotations. Interpretation of the data of both 1 and 2 implied that a combination of conformations, both in glycosidic and side-chain positions, could explain the experimental data. Solvents effects were important in influencing the conformational mix and averaged optical rotation. Three-bond heteronuclear coupling constants ³J_{C, H} were obtained for the glycosidic linkages of 1 and 2 in D₂O and DMSO. Analysis of the coupling constants with a Karplus curve showed that small reductions in the glycosidic torsion angles of the conformations of the models used here of ca. 10°–15° in ϕ and 5°–10° in ψ were required to give better agreement with experiment; a combination of conformations for both 1 and 2 was consistent with the data. There was a negligible influence on the coupling constants of 1 on changing the solvent from D₂O to DMSO. © 1997 John Wiley & Sons, Inc.     

NMR and conformational studies of linear and cyclic oligo-(1→6)-β-d-glucosamines

The conformational behavior of a series of linear and cyclic oligo-(1→6)-β-d-glucosamines and their N-acetylated derivatives, which are related to fragments of natural poly-N-acetylglucosamine, was studied by theoretical molecular modeling and experimental determination of transglycosidic vicinal coupling constants ³J_C,H and ³J_H,H. Molecular dynamics simulations were performed under several types of conditions varying in the consideration of ionization of amino groups, solvent effect, and temperature. Neural network clustering and asphericity calculations were performed on the basis of molecular dynamics data. It was shown that disaccharide fragments in the studied linear oligosaccharides were not rigid, and tended to have several conformers, thus determining the overall twisted shape with helical elements. In addition, it was found that the behavior of C5–C6 bond depended significantly upon the simulation conditions. The cyclic di-, tri-, and tetrasaccharides mostly had symmetrical ring-shaped conformations. The larger cycles tended to adopt more complicated shapes, and the conformational behavior of their disaccharide fragments was close to that in the linear oligosaccharides.     

Probing the U-Shaped Conformation of Caveolin-1 in a Bilayer

Caveolin induces membrane curvature and drives the formation of caveolae that participate in many crucial cell functions such as endocytosis. The central portion of caveolin-1 contains two helices (H1 and H2) connected by a three-residue break with both N- and C-termini exposed to the cytoplasm. Although a U-shaped configuration is assumed based on its inaccessibility by extracellular matrix probes, caveolin structure in a bilayer remains elusive. This work aims to characterize the structure and dynamics of caveolin-1 (D82–S136; Cav1_82–136) in a DMPC bilayer using NMR, fluorescence emission measurements, and molecular dynamics simulations. The secondary structure of Cav1_82–136 from NMR chemical shift indexing analysis serves as a guideline for generating initial structural models. Fifty independent molecular dynamics simulations (100 ns each) are performed to identify its favorable conformation and orientation in the bilayer. A representative configuration was chosen from these multiple simulations and simulated for 1 μs to further explore its stability and dynamics. The results of these simulations mirror those from the tryptophan fluorescence measurements (i.e., Cav1_82–136 insertion depth in the bilayer), corroborate that Cav1_82–136 inserts in the membrane with U-shaped conformations, and show that the angle between H1 and H2 ranges from 35 to 69°, and the tilt angle of Cav1_82–136 is 27 ± 6°. The simulations also reveal that specific faces of H1 and H2 prefer to interact with each other and with lipid molecules, and these interactions stabilize the U-shaped conformation.     

The beta-strand-loop-beta-strand conformation is marginally populated in beta2-microglobulin (20-41) peptide in solution as revealed by replica exchange molecular dynamics simulations

Solid-state NMR study shows that the 22-residue K3 peptide (Ser²⁰-Lys⁴¹) from β₂-microglobulin (β₂m) adopts a β-strand-loop-β-strand conformation in its fibril state. Residue Pro³² has a trans conformation in the fibril state of the peptide, while it adopts a cis conformation in the native state of full-length β₂m. To get insights into the structural properties of the K3 peptide, and determine whether the strand-loop-strand conformation is encoded at the monomeric level, we run all-atom explicit solvent replica exchange molecular dynamics on both the cis and trans variants. Our simulations show that the conformational space of the trans- and cis-K3 peptides is very different, with 1% of the sampled conformations in common at room temperature. In addition, both variants display only 0.3-0.5% of the conformations with β-strand-loop-β-strand character. This finding, compared to results on the Alzheimer's Aβ peptide, suggests that the biases toward aggregation leading to the β-strand-loop-β-strand conformation in fibrils are peptide-dependent.     

Modeling of ligand binding to G protein coupled receptors: cannabinoid CB1, CB2 and adrenergic β2AR

Cannabinoid and adrenergic receptors belong to the class A (similar to rhodopsin) G protein coupled receptors. Docking of agonists and antagonists to CB₁ and CB₂ cannabinoid receptors revealed the importance of a centrally located rotamer toggle switch and its possible participation in the mechanism of agonist/antagonist recognition. The switch is composed of two residues, F3.36 and W6.48, located on opposite transmembrane helices TM3 and TM6 in the central part of the membranous domain of cannabinoid receptors. The CB₁ and CB₂ receptor models were constructed based on the adenosine A_2A receptor template. The two best scored conformations of each receptor were used for the docking procedure. In all poses (ligand-receptor conformations) characterized by the lowest ligand-receptor intermolecular energy and free energy of binding the ligand type matched the state of the rotamer toggle switch: antagonists maintained an inactive state of the switch, whereas agonists changed it. In case of agonists of β₂AR, the (R,R) and (S,S) stereoisomers of fenoterol, the molecular dynamics simulations provided evidence of different binding modes while preserving the same average position of ligands in the binding site. The (S,S) isomer was much more labile in the binding site and only one stable hydrogen bond was created. Such dynamical binding modes may also be valid for ligands of cannabinoid receptors because of the hydrophobic nature of their ligand-receptor interactions. However, only very long molecular dynamics simulations could verify the validity of such binding modes and how they affect the process of activation.     

Threshold Occupancy and Specific Cation Binding Modes in the Hammerhead Ribozyme Active Site are Required for Active Conformation

The relationship between formation of active in-line attack conformations and monovalent (Na⁺) and divalent (Mg²⁺) metal ion binding in hammerhead ribozyme (HHR) has been explored with molecular dynamics simulations. To stabilize repulsions between negatively charged groups, different requirements of the threshold occupancy of metal ions were observed in the reactant and activated precursor states both in the presence and in the absence of a Mg²⁺ in the active site. Specific bridging coordination patterns of the ions are correlated with the formation of active in-line attack conformations and can be accommodated in both cases. Furthermore, simulation results suggest that the HHR folds to form an electronegative recruiting pocket that attracts high local concentrations of positive charge. The present simulations help to reconcile experiments that probe the metal ion sensitivity of HHR catalysis and support the supposition that Mg²⁺, in addition to stabilizing active conformations, plays a specific chemical role in catalysis.     

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α₂ and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α₄-α₅ KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).     

Multiple interconverting conformers of the cyclic tetrapeptide tentoxin, [cyclo-(L-MeAla1-L-Leu2-MePhe[(Z)Δ]3-Gly4)], as seen by two-dimensional 1H-nmr spectroscopy

The conformations of the phytotoxic cyclic tetrapeptide tentoxin [cyclo-(L -MeAla¹-L -Leu²-MePhe[(Z)Δ]³-Gly⁴ )] have been studied in aqueous solution by two-dimensional proton nmr at various temperatures. Contrary to what is observed in chloroform, tentoxin exhibits multiple exchanging conformations in water. Aggregation phenomena were also observed. Four conformations with different proportions (51, 37, 8, and 4%) were observed at ?5°C. Models were constructed from nmr parameters and restrained molecular dynamics simulations. All the models exhibit cis-trans-cis-trans conformation of the amide bond sequence. The conversion from one form to another is accomplished by a conformational peptide flip consisting of a 180° rotation of a nonmethylated peptide bond. © 1995 John Wiley & Sons, Inc.     

Molecular dynamics simulation of amyloid beta dimer formation

Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.