A critical evaluation of models for complex molecular dynamics: application of NMR studies of double- and single-stranded DNA

The nature of internal and overall motions in native (double-stranded) and denatured (single-stranded) DNA fragments 120–160 base pairs (bp) long is examined by molecular-dynamics modeling using ¹³C-nmr spin-relaxation data obtained over the frequency range of 37–125 MHz. The broad range of ¹³C frequencies is required to differentiate among various models. Relatively narrow linewidths, large nuclear Overhauser enhancements (NOEs), and short T₁ values all vary significantly with frequency and indicate the presence of rapid, restricted internal motions on the nanosecond time scale. For double-stranded DNA monomer fragments (147 bp, 24 Å diam at 32°C), the overall motion is that of an axially symmetric cylinder (τ_x = ～10^?6 s;τ_Z = ～1.8 × 10^?8s), which is in good agreement with values calculated from hydrodynamic theory (τ_x = ～1.8 × 10^?6 s; τ_Z = ～2.7 × 10^?8 s). The DNA internal motion can be modeled as restricted amplitude internal diffusion of individual C? H vectors of deoxyribose methine carbons C1′, C3′, and C4′, either with conic boundary conditions (τ_w = ～4 × 10^?9 s, θ_cone = ～21°) or as a bistable jump (τ_A = τ_B = ～2 × 10^?9 s, θ = ～15°). We discuss the critical role in molecular-dynamics modeling played by the angle (β) that individual C? H vectors make with the long axis of the DNA helix. Heat denaturation brings about increases in both the rate and amplitude of the internal motion (described by the wobble model with τ_W = ～0.2 × 10^?9 s, θ_cone = ～50°), and overall motion is affected by becoming essentially isotropic (τ_x = τ_Z = ～5 × 10^?8 s) for the single-stranded molecules. Since ¹³C-nmr data obtained at various DNA concentrations for C2′ of the deoxyribose ring is not described well by the above models, a new model incorporating an additional internal motion is proposed to take into account the rapid, extensive, and weakly coupled motion of C2′.     

3D-Structure of the interior fusion peptide of HGV/GBV-C by 1H NMR, CD and molecular dynamics studies

In this work, we present a structural characterization of the putative fusion peptide E2(279-298) corresponding to the E2 envelope protein of the HGV/GBV-C virus by (1)H NMR, CD and MD studies performed in H(2)O/TFE and in lipid model membranes. The peptide is largely unstructured in water, whereas in H(2)O/TFE and in model membranes it adopts an helical structure (approximately 65-70%). The partitioning free energy DeltaG ranges from -6 to -7.5 kcal mol(-1). OCD measurements on peptide-containing hydrated and oriented lipid multilayers showed that the peptide adopts a predominantly surface orientation. The (1)H NMR data (observed NOEs, deuterium exchange rates, Halpha chemical shift index and vicinal coupling constants) and the molecular dynamics calculations support the conclusions that the peptide adopts a stable helix in the C-terminal 9-18 residues slightly inserted into the lipid bilayer and a major mobility in the amino terminus of the sequence (1-8 residues).     

NMR dynamics and antibody recognition of the meningococcal lipidated outer membrane protein LP2086 in micellar solution

Neisseria meningitidis is a major cause of meningitis. Although protective vaccination is available against some pathogenic serogroups, serogroup B meningococci have been a challenge for vaccinologists. A family of outer membrane lipoproteins, LP2086 (or factor H binding proteins, fHbp), has been shown to elicit bactericidal antibodies and is currently part of a cocktail vaccine candidate. The NMR structure of the variant LP2086-B01 in micellar solution provided insights on the topology of this family of proteins on the biological membrane. Based on flow cytometry experiments on whole meningococcal cells, binding experiments with monoclonal antibodies, and the NMR structure in micellar solution, we previously proposed that LP2086-B01 anchors the outer bacterial membrane through its lipidated N-terminal cysteine, while a flexible 20 residue linker positions the protein above the layer of lipo-oligosaccharides that surrounds the bacteria. This topology was suggested to increase the antigen exposure to the immune system. In the present work, using micellar solution as a membrane mimicking system, we characterized the backbone dynamics of the variant LP2086-B01 in both its lipidated and unlipidated forms. In addition, binding experiments with a Fab fragment derived from the monoclonal MN86-1042-2 were also performed. Our data suggests that due to the length and flexibility of the N-terminal linker, the antigen is not in contact with the micelle, thus making both N- and C-domains highly available to the host immune system. This dynamic model, combined with the binding data obtained with MN86-1042-2, supports our previously proposed arrangement that LP2086-B01 exposes one face to the extracellular space. Binding of MN86-1042-2 antibody shows that the N-domain is the primary target of this monoclonal, providing further indication that this domain is immunologically important for this family of proteins.     

Larvicidal activity of novel anthraquinone analogues and their molecular docking studies

To investigate the larvicidal activities of novel anthraquinones (1a-1k) against Culex quinquefasciatus mosquito larvae. Novel anthraquinones (1a-1k) derivatives were synthesis via condensation method. The compounds were confirmed through FT-IR spectroscopy, ¹H & ¹³C NMR spectrum, and mass spectral studies. The larvicidal activity of compound 1c was highly active LD₅₀ 20.92 µg/mL against Culex quinquefasciatus compared standard permethrin with LD₅₀ 25.49 µg/mL. Molecular docking studies were carried out for compound 1c against Odorant-binding protein of Culex quinquefasciatus. The compound 1c (−9.8 Kcal/mol) was a potent larvicide with more binding energy than control permethrin (−9.7 Kcal/mol). Therefore, compound (1c) may be more significant inhibitors of mosquito larvicidal.     

Structure and dynamics of a fluorescent DNA oligomer containing the EcoRI recognition sequence: fluorescence, molecular dynamics, and NMR studies

The self-complementary DNA decamer duplex d(CTGAATTCAG)2 and its modified counterpart d(CTGA[2AP]TTCAG)2, where the innermost adenine (6-aminopurine) has been replaced with the fluorescent analogue 2-aminopurine (2AP), have been studied by fluorescence and NMR spectroscopy and simulated by molecular dynamics. Both decamers are recognized and cleaved by the EcoRI restriction endonuclease. 2D NMR results show that both decamers have a standard B-type conformation below 20 degrees C, though a disturbance exists to the 5' side of the 2AP site which may originate from increased local mobility. The fluorescence and fluorescence anisotropy decays of both decamers, as well as the one containing 2AP in only one chain, were studied as a function of temperature. The data show that the 2AP base exists in a temperature-dependent distribution of states and shows rapid motions, suggesting interconversion among these states on a time scale of about 10(-10) s. The integrated fluorescence of the decamer with 2AP in both chains shows a large increase around the helix melting temperature whereas the decamer with one 2AP shows only a mild increase, showing that the mixed helix has a different structural transition as sensed by the 2AP base. The data suggest a model of conformational states which have distinct fluorescence decay times. The various states may differ in the degree of base stacking. Fluctuations in the degree of stacking of the A or 2AP base are supported by molecular dynamics simulations, which additionally show that the 2AP-T or A-T base pair hydrogen bonds remain intact during these large motions.     

Conformational analysis of peptide analogues of silkmoth chorion protein segments using CD, NMR and molecular modelling.

Silkmoth proteins secreted from the follicular cells that surround the oocyte form a large extracellular assembly which is important for protecting and sustaining the structure of the oocyte and the developing embryo. These proteins have been classified into two major families (A and B). Sequence analysis showed conservation of a central domain containing long stretches of six amino acid residue repeats in both families, which have been suggested to be organized in beta-sheet structures. In this work NMR and CD spectra, as well as molecular calculations, have been used to investigate the conformational properties of two synthetic peptides (A and B), analogues of parts of the central domain of silkmoth chorion proteins of the A and B families, respectively. These peptides consist of three tandem repeats of the six-residue basic motif. Analysis of CD spectra of the two peptides in aqueous solutions and mixtures with organic solvents revealed beta-sheet and turn structural elements with a percentage higher than 40%. NOESY spectra at low temperatures (263-273 K) show sequential nOe connectivities (i, i + 1), indicative of a relative flexibility. The presence of HNi-HNi+1 cross-peaks and medium Halphai-HNi+1 connectivities, chemical shift deviations and temperature coefficient data provide, for the first time, experimental evidence that local folded structures around Gly residues occur in peptide segments of chorion proteins in solution. Simulated annealing calculations were used to examine the conformational space of the peptides and to probe the initial steps of amyloid fibril formation in the case of chorion proteins.     

Conformation, orientation and dynamics of dodecylphosphocholine in micellar aggregate: a 3.2 ns molecular dynamics simulation study

A dodecylphosphocholine micelle of 86 monomers with 5776 water molecules has been simulated under NPT conditions for 3.2 ns using GROMACS2.0. The micelle was found to be very dynamic. Some of the C-C bonds, independent of their position in the DPC monomer, adopt gauche conformation and the trans <--> gauche transitions are quite frequent. An average of about 11% of the C-C bonds in the micelle are observed to be in the gauche conformation (i.e., |dihedral angle|< 120 degrees). The terminal methyl groups are randomly distributed all over the micelle whereas the nitrogen atom of phosphocholine headgroup atoms is restricted to the interface region. Some of the monomers were found to lie on the surface. The shape of micelle, influenced by the packing considerations, shows deviations from spherical shape. The phosphocholine headgroup is well solvated and there is no water penetration into the micelle core. The overall features of the micelle of 86 DPC monomers conforms to the lattice model of micelle proposed by Dill and Flory [Dill K A, Flory P J (1981) Proc Natl Acad Sci USA 78, 676-680] and is similar to DPC micelles of smaller aggregate sizes except for the positional preference of the C-C bonds for the gauche conformation and the trans<-->gauche transition times [Tieleman D P, van der Spoel D, Berendsen H J C (2000) J Phys Chem B 104, 6380-6388; Wymore T, Gao X F, Wong T C (1999) J Mol Struct (Theochem) 485-486, 195-210]. It appears that packing considerations play a predominant role in determining the shape and dynamics of the micelle.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

Molecular modeling studies of atorvastatin analogues as HMGR inhibitors using 3D-QSAR,molecular docking and molecular dynamics simulations

3-Hydroxy-3-methylglutaryl coenzyme-A reductase (HMGR) is generally regarded as targets for the treatment of hypercholesterolemia. HMGR inhibitors (more commonly known as statins) are discovered as plasma cholesterol lowering molecules. In this work, 120 atorvastatin analogues were studied using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulation. The results show that the best CoMFA (comparative molecular field analysis) model has q² = 0.558 and r² = 0.977, and the best CoMSIA (comparative molecular similarity indices analysis) model has q² = 0.582 and r² = 0.919. Molecular docking and MD simulation explored the binding relationship of the ligand and the receptor protein. The calculation results indicated that the hydrophobic and electrostatic fields play key roles in QSAR model. After MD simulation, we found four vital residues (Lys735, Arg590, Asp690 and Asn686) and three hydrophobic regions in HMGR binding site. The calculation results show that atorvastatin analogues obtained by introduction of F atoms or gem-difluoro groups could obviously improve the inhibitory activity. The new HMGR inhibitor analogues design in this Letter had been submitted which is being currently synthesized by our laboratories.     

NMR and molecular dynamics studies of the interaction of melatonin with calmodulin

Pineal hormone melatonin (N-acetyl-5-methoxytryptamine) is thought to modulate the calcium/calmodulin signaling pathway either by changing intracellular Ca(2+) concentration via activation of its G-protein-coupled membrane receptors, or through a direct interaction with calmodulin (CaM). The present work studies the direct interaction of melatonin with intact calcium-saturated CaM both experimentally, by fluorescence and nuclear magnetic resonance spectroscopies, and theoretically, by molecular dynamics simulations. The analysis of the experimental data shows that the interaction is calcium-dependent. The affinity, as obtained from monitoring (15)N and (1)H chemical shift changes for a melatonin titration, is weak (in the millimolar range) and comparable for the N- and C-terminal domains. Partial replacement of diamagnetic Ca(2+) by paramagnetic Tb(3+) allowed the measurement of interdomain NMR pseudocontact shifts and residual dipolar couplings, indicating that each domain movement in the complex is not correlated with the other one. Molecular dynamics simulations allow us to follow the dynamics of melatonin in the binding pocket of CaM. Overall, this study provides an example of how a combination of experimental and theoretical approaches can shed light on a weakly interacting system of biological and pharmacological significance.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Polypeptide models of elastin: CD and NMR studies on synthetic poly(X-Gly-Gly).

Poly(X-Gly-Gly), simple structural models for the hydrophobic, proline-devoid, regions of elastin, have been synthesized and studied by circular dichroism and NMR spectroscopies. The results gave evidence of type II beta-turns as the only ordered structure present in the polymers. The stability of the turns has been shown to decrease on hydration and to increase in the series Leu less than Ala less than Val less than Ile.     

Intermolecular interaction of voriconazole analogues with model membrane by DSC and NMR,and their antifungal activity using NMR based metabolic profiling

The development of novel antifungal agents with high susceptibility and increased potency can be achieved by increasing their overall lipophilicity. To enhance the lipophilicity of voriconazole, a second generation azole antifungal agent, we have synthesized its carboxylic acid ester analogues, namely p-methoxybenzoate (Vpmb), toluate (Vtol), benzoate (Vbz) and p-nitrobenzoate (Vpnb). The intermolecular interactions of these analogues with model membrane have been investigated using nuclear magnetic resonance (NMR) and differential scanning calorimetric (DSC) techniques. The results indicate varying degree of changes in the membrane bilayer’s structural architecture and physico-chemical characteristics which possibly can be correlated with the antifungal effects via fungal membrane. Rapid metabolite profiling of chemical entities using cell preparations is one of the most important steps in drug discovery. We have evaluated the effect of synthesized analogues on Candida albicans. The method involves real time ¹H NMR measurement of intact cells monitoring NMR signals from fungal metabolites which gives Metabolic End Point (MEP). This is then compared with Minimum Inhibitory Concentration (MIC) determined using conventional methods. Results indicate that one of the synthesized analogues, Vpmb shows reasonably good activity.     

Conformational properties of oxytocin in dimethyl sulfoxide solution: NMR and restrained molecular dynamics studies.

The conformation of oxytocin, the neurohypophyseal nonapeptide hormone, in solution in deuterated dimethyl sulfoxide has been determined by 1H-nmr. The structural determination is based on the experimental data set of nuclear Overhauser effect restraints. Obtained after the restrained molecular dynamics simulation on an initial structure of extended conformation, five resultant structures satisfy the experimental restraints well. These structures resemble that of the crystal structure of deamino-oxytocin, an analogue of oxytocin, in terms of a close correlation observed both at two beta-turn regions of the 20-membered tocin ring and at the tripeptide tail end. Based on this comparison and analysis of restrained molecular dynamics trajectories, we found that, although the turns are stabilized by the formation of hydrogen bonds, the oxytocin molecule possesses a slight twist in DMSO solution relative to the orientation of deamino-oxytocin in the crystalline state. Analyses of oxytocin conformation indicate that the tripeptide tail is more flexible than the tocin ring.     

Fatty acid synthesis in Xylella fastidiosa: correlations between genome studies, 13C NMR data, and molecular models

Xylella fastidiosa was the first plant pathogen to have its complete genome sequence elucidated. Routine database analyses suggested that two enzymes essential for fatty acid synthesis were missing, one of these is the holo-acyl-carrier-protein synthase. However, here we demonstrate, using (13)C NMR spectroscopy, that X. fastidiosa is indeed able to synthesize fatty acids from acetate via an apparently conventional metabolic pathway. We further identify a gene product HetI, an alternative phosphopantetheinyl transferase, which we propose to fill the missing link. Homology modeling of HetI shows conservation of the Coenzyme A binding site suggesting it to be an active enzyme and reveals several interesting structural features when compared with the surfactin synthase-activating enzyme, on which the model was built. These include a simplified topology due to N- and C-terminal deletions and the observation of a novel serine ladder.     

Iron affects the structure of cell membrane molecular models

The effects of Fe(3+) and Fe(2+) on molecular models of biomembranes were investigated. These consisted of bilayers of dimyristoylphosphatidylcholine (DMPC) and of dimyristoylphosphatidylethanolamine (DMPE), classes of phospholipids located in the outer and inner moieties of cell membranes, respectively. X-ray studies showed that very low concentrations of Fe(3+) affected DMPC organization and 10(-3)M induced a total loss of its multilamellar periodic stacking. Experiments carried out with Fe(2+) on DMPC showed weaker effects than those induced by Fe(3+) ions. Similar experiments were performed on DMPE bilayers. Fe(3+) from 10(-7)M up to 10(-4)M had practically no effect on DMPE structure. However, 10(-3)M Fe(3+) induced a deep perturbation of the multilamellar structure of DMPE. However, 10(-3)M Fe(2+) had no effect on DMPE organization practically. Differential scanning calorimetry measurements also revealed different effects of Fe(3+) and Fe(2+) on the phase transition and other thermal properties of the examined lipids. In conclusion, the results obtained indicate that iron ions interact with phospholipid bilayers perturbing their structures. These findings are consistent with the observation that iron ions change cell membrane fluidity and, therefore, affect its functions.     

Conformation of secretin in dimethyl sulfoxide solution. NMR studies and restrained molecular dynamics

The solution conformation of the 27-residue polypeptide hormone secretin in dimethyl sulfoxide has been determined on the basis of 1H-NMR measurements. The experimental data set used in the structure determination consisted of 98 nuclear-Overhauser-enhancement-derived interproton and dihedral angle restraints from coupling constants. The NH-NH and H alpha-NH NOEs were determined from build-up rates, while the remaining distances were classified in a qualitative manner. The structure calculations consisted of two phases. First, dynamical simulated annealing calculations were carried out to find conformations of the peptide which satisfy NOE and phi dihedral restraints. The convergence of ten calculated structures was good except for those regions of the molecule where NOE data were not unambiguous. From the calculated set another initial structure was built which was again minimized in several 5-ps calculations now employing the full empirical energy function. The resulting structures of secretin reveal conformationally well-defined regions, but not a single uniform secondary structure. The structure is different from the calculated structure from trifluoroethanol/water measurements.     

Solution structures of cyclic and dicyclic analogues of growth hormone releasing factor as determined by two-dimensional NMR and CD spectroscopies and constrained molecular dynamics.

Solution structures were determined for a linear analogue of growth hormone releasing factor (GRF), and cyclic and dicyclic analogues in which the side chains of aspartyl and lysyl residues spaced at positions i-(i + 4) were joined to form a lactam. The four analogues were [Ala15]-GRF-(1-29)-NH2 and its cyclo8-12, cyclo21-25, and dicyclo8-12;21-25 derivatives. The peptides were studied in two solvent systems: 75% methanol/25% water at pH 6.0; and 100% water at pH 3.0. CD spectroscopy was used to assess the overall alpha-helical content. Nuclear magnetic resonance spectroscopy was used to determine the structures in more detail. Nearly complete proton resonance assignments were made for each of the peptides, in both solvents. Nuclear Overhauser effects were converted into distance constraints and applied in the molecular dynamics program CHARMM to evaluate the range of low-energy structures that satisfied the nmr data. In 75% methanol, all of the peptides are comprised of a single alpha-helical segment with fraying of one to three residues at each end. The linear analogue has a tendency to kink. In water, the analogues have two helical segments with flexible regions between them and at the termini of the peptides. The linear analogue is helical at residues 7-14 and 21-28. In the cyclo8-12 analogue, the N-terminal helical region extends to include residues 7-19, while the other helical region is slightly shortened. In the cyclo21-25 analogue, the C-terminal helical region is extended to include residues 19-28, while the N-terminal helical region is destabilized. The dicyclic analogue has the largest N-terminal helix, spanning residues 7-20, but its helical segment at residues 21-28 is not well ordered. All of the analogues exhibit substantial biological activity. The cyclic and dicyclic analogues show dramatically increased resistance to degradation during incubation with human plasma. The i-(i + 4) lactam, therefore, appears to be a synthetic means of stabilizing a local alpha-helical conformation, which may be of general use in the design of active, stable peptides.     

NMR and molecular dynamics studies of tachykinins: conformation of substance P fragment 4-11

The conformation of the C-terminal octapeptide fragment of Substance P (SP4-11, Pro-Gln-Gln-Phe-Phe-Gly-Leu-Met-NH2) has been investigated by 2D-NMR and MD methods. The octapeptide exists in a blend of conformations. The molecule seems to shuttle between conformations with gamma-bends either at Phe5 or Gly6 or Gln3 or Leu7 and between a nearly extended structure.