Analysis of the conformational features of Watson–Crick duplex fragments by molecular mechanics and quantum mechanics methods

It is generally accepted that important features of the Watson–Crick duplex (WCD) originate from the molecular structure of its subunits. However, it is still unclear what properties of each subunit are responsible for significant features of the WCD structure. Computations of deoxydinucleoside monophosphate (dDMP) complexes with Na ions on the basis of the density functional theory (DFT) have shown that conformational properties of minimal single-stranded fragments of DNA play a pivotal role in the origin of the unique features of the WCD. The directionality of the sugar-phosphate backbone (SPB) and preferable ranges of its torsion angles combined with the difference in geometry between purines and pyrimidines have been found to define the nucleotide sequence dependence of the WCD three-dimensional structure. In this work, density functional theory computations were extended to minimal duplex fragments, that is, complementary dDMP (cdDMP) complexes with Na ions. Using several computational methods and various functionals, energy minima were searched for the BI conformation of cdDMP complexes with different nucleotide sequences. Two sequences were optimized using an ab initio method at the MP2/6-31++G** level of theory. An analysis of the SPB torsion angles, sugar-ring puckering, and mutual base positions in the optimized structures showed that the conformational features of cdDMP complexes with Na ions remained within the BI ranges and become more similar to the corresponding features that WCDs display in a crystal. Qualitatively, the main features of each cdDMP complex were invariant with different computational methods, although the values of certain conformational parameters could vary, but still within the limits that are typical of the corresponding family. Common functionals that are employed in DFT calculations were observed to overestimate the distance between base pairs, while MP2 computations and new complex functionals yielded structures with atom–atom contacts that are too close. Several energy minima that correspond to the BI conformation have been proven to exist for certain cdDMP complexes with Na ions, indicating that the topography of the potential energy surface is complex. This circumstance accounts for the variation of conformational parameters among duplex fragments with the same nucleotide sequence. The common AMBER and CHARMM molecular mechanics force fields reproduce many conformational characteristics of dDMPs and their complementary complexes with Na ions, but fail to reproduce certain details of the nucleotide sequence dependence of the WCD conformation.     

Modelling and predicting the binding mechanics of HIV P1053-0.5β antibody complex

Simulating antigen–antibody interactions is crucial in understanding the mechanics of antigen–antibody binding in medical science. In this study, molecular dynamics simulations are performed to analyse the dissociation of the P1053-0.5β antibody complex structure. The two-dimensional free energy profiles of the complex structure are extracted using the weighted histogram analysis method, and the binding pathway is then predicted using a modified form of the MaxFlux-PRM method. The simulation results suggest that 10 amino residues (i.e. Leu11, Val13, Asp34, Arg112, Thr101, Gly127, Val229, Ser231, Ile235 and Arg236) play a key role in relaxing the antibody structure, thereby facilitating the binding of the 0.5β antibody-P1053 peptide system.     

Investigations of Takeout proteins’ ligand binding and release mechanism using molecular dynamics simulation

Takeout (To) proteins exist in a diverse range of insect species. They are involved in many important processes of insect physiology and behaviors. As the ligand carriers, To proteins can transport the small molecule to the target tissues. However, ligand release mechanism of To proteins is unclear so far. In this contribution, the process and pathway of the ligand binding and release are revealed by conventional molecular dynamics simulation, steered molecular dynamics simulation and umbrella sampling methods. Our results show that the α4-side of the protein is the unique gate for the ligand binding and release. The structural analysis confirms that the internal cavity of the protein has high rigidity, which is in accordance with the recent experimental results. By using the potential of mean force calculations in combination with residue cross correlation calculation, we concluded that the binding between the ligand and To proteins is a process of conformational selection. Furthermore, the conformational changes of To proteins and the hydrophobic interactions both are the key factors for ligand binding and release.     

The molecular assembly of amyloid aβ controls its neurotoxicity and binding to cellular proteins

Accumulation of β-sheet-rich peptide (Aβ) is strongly associated with Alzheimer's disease, characterized by reduction in synapse density, structural alterations of dendritic spines, modification of synaptic protein expression, loss of long-term potentiation and neuronal cell death. Aβ species are potent neurotoxins, however the molecular mechanism responsible for Aβ toxicity is still unknown. Numerous mechanisms of toxicity were proposed, although there is no agreement about their relative importance in disease pathogenesis. Here, the toxicity of Aβ 1-40 and Aβ 1-42 monomers, oligomers or fibrils, was evaluated using the N2a cell line. A structure-function relationship between peptide aggregation state and toxic properties was established. Moreover, we demonstrated that Aβ toxic species cross the plasma membrane, accumulate in cells and bind to a variety of internal proteins, especially on the cytoskeleton and in the endoplasmatic reticulum (ER). Based on these data we suggest that numerous proteins act as Aβ receptors in N2a cells, triggering a multi factorial toxicity.     

Investigation of the binding mode of (−)-meptazinol and bis-meptazinol derivatives on acetylcholinesterase using a molecular docking method

Molecular docking has been performed to investigate the binding mode of (-)-meptazinol (MEP) with acetylcholinesterase (AChE) and to screen bis-meptazinol (bis-MEP) derivatives for preferable synthetic candidates virtually. A reliable and practical docking method for investigation of AChE ligands was established by the comparison of two widely used docking programs, FlexX and GOLD. In our hands, we had more luck using GOLD than FlexX in reproducing the experimental poses of known ligands (RMSD<1.5 A). GOLD fitness values of known ligands were also in good agreement with their activities. In the present GOLD docking protocol, (-)-MEP seemed to bind with the enzyme catalytic site in an open-gate conformation through strong hydrophobic interactions and a hydrogen bond. Virtual screening of a potential candidate compound library suggested that the most promising 15 bis-MEP derivatives on the list were mainly derived from (-)-MEP with conformations of (S,S) and (SR,RS) and with a 2- to 7-carbon linkage. Although there are still no biological results to confirm the predictive power of this method, the current study could provide an alternate tool for structural optimization of (-)-MEP as new AChE inhibitors. [Figure: see text].     

Synthesis of novel inhibitors of β-glucuronidase based on benzothiazole skeleton and study of their binding affinity by molecular docking

Benzothiazole derivatives 1-26 have been synthesized and their in vitro β-glucuronidase potential has been evaluated. Compounds 4 (IC(50)=8.9 ± 0.25 μM), 5 (IC(50)=36.1 ± 1.80 μM), 8 (IC(50)=8.9 ± 0.38 μM), 13 (IC(50)=19.4 ± 1.00 μM), 16 (IC(50)=4.23 ± 0.054 μM), and 18 (IC(50)=2.26 ± 0.06 μM) showed β-glucuronidase activity potent than the standard (d-saccharic acid 1,4-lactone, IC(50)=48.4 ± 1.25 μM). Compound 9 (IC(50)=94.0 ± 4.16 μM) is found to be the least active among the series. All active analogs were also evaluated for cytotoxicity and none of the compounds showed any cytotoxic effect. Furthermore, molecular docking studies were performed using the gold 3.0 program to investigate the binding mode of benzothiazole derivatives. This study identifies a novel class of β-glucuronidase inhibitors.     

Nucleotide binding domain 1 pharmacophore modeling for visualization and analysis of P-glycoprotein–flavonoid molecular interactions

Background

P-glycoprotein (P-gp) is a 170-kDa membrane protein. It provides a barrier function and help to excrete toxins from the body as a transporter. Some bioflavonoids have been shown to block P-gp activity.

Objective

To evaluate the important amino acid residues within nucleotide binding domain 1 (NBD1) of P-gp that play a key role in molecular interactions with flavonoids using structure-based pharmacophore model.

Methods

In the molecular docking with NBD1 models, a putative binding site of flavonoids was proposed and compared with the site for ATP. The binding modes for ligands were achieved using LigandScout to generate the P-gp–flavonoid pharmacophore models.

Results

The binding pocket for flavonoids was investigated and found these inhibitors compete with the ATP for binding site in NBD1 including the NBD1 amino acid residues identified by the in silico techniques to be involved in the hydrogen bonding and van der Waals (hydrophobic) interactions with flavonoids.

Conclusion

These flavonoids occupy with the same binding site of ATP in NBD1 proffering that they may act as an ATP competitive inhibitor.

     

Binding of γ-crystallin substrate prevents the binding of copper and zinc ions to the molecular chaperone α-crystallin

α-Crystallin is a small heat shock protein and molecular chaperone. Binding of Cu2+ and Zn2+ ions to α-crystallin leads to enhanced chaperone function. Sequestration of Cu2+ by α-crystallin prevents metal-ion mediated oxidation. Here we show that binding of human γD-crystallin (HGD, a natural substrate) to human αA-crystallin (HAA) is inversely related to the binding of Cu2+/Zn2+ ions: The higher the amount of bound HGD, the lower the amount of bound metal ions. Thus, in the aging lens, depletion of free HAA will not only lower chaperone capacity but also lower Cu2+ sequestration, thereby promoting oxidation and cataract.     

A molecular dynamics simulation study of polyamine– and sodium–DNA. Interplay between polyamine binding and DNA structure

Four different molecular dynamics (MD) simulations have been performed for infinitely long ordered DNA molecules with different counterions, namely the two natural polyamines spermidine(3+) (Spd³⁺) and putrescine(2+) (Put²⁺), the synthetic polyamine diaminopropane(2+) (DAP²⁺), and the simple monovalent cation Na⁺. All systems comprised a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charge. The simulation setup mimics the DNA state in oriented DNA fibers, previously studied using NMR and other experimental methods. In this paper the interplay between polyamine binding and local DNA structure is analyzed by investigating how and if the minor groove width of DNA depends on the presence and dynamics of the counterions. The results of the MD simulations reveal principal differences in the polyamine–DNA interactions between the natural [spermine(4+), Spd³⁺, Put²⁺] and the synthetic (DAP²⁺) polyamines.Abbreviations DAP diaminopropane - DDD Drew–Dickerson dodecamer - MD molecular dynamics - Put putrescine - RDF radial distribution function - Spd spermidine - Spm spermine     

Probing the “fingers” domain binding pocket of Hepatitis C virus NS5B RdRp and D559G resistance mutation via molecular docking,molecular dynamics simulation and binding free energy calculations

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the “fingers” sub-domain of NS5B. The “fingers” domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the “fingers” sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π?π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN?Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50?ns MD simulation, the π?π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance.

Communicated by Ramaswamy H. Sarma     

Biophysical insights into the binding characteristics of bovine serum albumin with dipyridamole and the influence of molecular interaction with β cyclodextrin

Abstract

The binding characteristic of anti-platelet drug dipyridamole has been investigated with a transport protein, serum albumin. A multi-spectroscopic approach has been employed, and the results were well supported by in silico molecular docking and simulation studies. The fluorescence quenching of serum albumin at three different temperatures revealed that the mechanism involved is static and the binding constant of the interaction was found to be of the order of 10⁴ M^?1. The reaction was found to be spontaneous and involved hydrophobic interactions. Synchronous, 3D fluorescence and CD spectroscopy indicated a change in conformation of bovine serum albumin (BSA) on interaction with DP. Using site-selective markers, the binding site of DP was found to be in subdomain IB. Molecular docking studies further corroborated these results. Molecular dynamic (MD) simulations showed lower RMSD values on interaction, suggesting the existence of a stable complex between DP and BSA. Furthermore, since β-Cyclodextrin (βCD) is used to improve the solubility of DP in ophthalmic solutions, therefore, the effect of (βCD) on the interaction of BSA and DP was also studied, and it was found that in the presence of βCD, the binding constant for BSA-DP interaction decreased. The present study is an attempt to characterize the transport of DP and to improve its bioavailability, consequently helping in dosage design to achieve optimum therapeutic levels.

Communicated by Ramaswamy H. Sarma     

Steered molecular dynamics simulation of the binding of the β2 and β3 regions in domain-swapped human cystatin C dimer

The crystal structure of the human cystatin C (hCC) dimer revealed that a stable twofold-symmetric dimer was formed via 3D domain swapping. Domain swapping with the need for near-complete unfolding has been proposed as a possible route for amyloid fibril initiation. Thus, the interesting interactions that occur between the two molecules may be important for the further aggregation of the protein. In this work, we performed steered molecular dynamics (SMD) simulations to investigate the dissociation of the β2 and β3 strands in the hCC dimer. The energy changes observed during the SMD simulations showed that electrostatic interactions were the dominant interactions involved in stabilizing the two parts of the dimer during the early stages of SMD simulation, whereas van der Waals (VDW) interactions and electrostatic interactions were equally matched during the latter stages. Furthermore, our data indicated that the two parts of the dimer are stabilized by intermolecular hydrogen bonds among the residues Arg51 (β2), Gln48 (β2), Asp65 (β3), and Glu67 (β3), salt bridges among the residues Arg53 (β2), Arg51 (β2), and Asp65 (β3), and VDW interactions among the residues Gln48 (β2), Arg51 (β2), Glu67 (β3), Asp65 (β3), Phe63 (β3), and Asn61 (β3). The residues Gln48 (β2), Arg51 (β2), Asp65 (β3) and Glu67 (β3) appear to be crucial, as they play important roles in both electrostatic and VDW interactions. Thus, the present study determined the key residues involved in the stabilization of the domain-swapped dimer structure, and also provided molecular-level insights into the dissociation process of the hCC dimer.     

GRIMM — an interactive personal-computer graphics interface to molecular mechanics

GRIMM, an interactive personal-computer graphics program designed to be used as an interface to the molecular mechanics program MM2, is described. It facilitates the interactive preparation of input data for MM2 and the visualization of the resultant structure. GRIMM can also be used to construct, modify and manipulate molecular structures shown on a screen.     

Probing the binding mechanisms of α-tocopherol to trypsin and pepsin using isothermal titration calorimetry,spectroscopic, and molecular modeling methods

α-Tocopherol is a required nutrient for a variety of biological functions. In this study, the binding of α-tocopherol to trypsin and pepsin was investigated using isothermal titration calorimetry (ITC), steady-state and time-resolved fluorescence measurements, circular dichroism (CD) spectroscopy, and molecular modeling methods. Thermodynamic investigations reveal that α-tocopherol binds to trypsin/pepsin is synergistically driven by enthalpy and entropy. The fluorescence experimental results indicate that α-tocopherol can quench the fluorescence of trypsin/pepsin through a static quenching mechanism. The binding ability of α-tocopherol with trypsin/pepsin is in the intermediate range, and one molecule of α-tocopherol combines with one molecule of trypsin/pepsin. As shown by circular dichroism (CD) spectroscopy, α-tocopherol may induce conformational changes of trypsin/pepsin. Molecular modeling displays the specific binding site and gives information about binding forces and α-tocopherol-tryptophan (Trp)/tyrosine (Tyr) distances. In addition, the inhibition rate of α-tocopherol on trypsin and pepsin was studied. The study provides a basic data set for clarifying the binding mechanisms of α-tocopherol with trypsin and pepsin and is helpful for understanding its biological activity in vivo.     

Co-crystal structures of PKG Iβ (92-227) with cGMP and cAMP reveal the molecular details of cyclic-nucleotide binding

Background

Cyclic GMP-dependent protein kinases (PKGs) are central mediators of the NO-cGMP signaling pathway and phosphorylate downstream substrates that are crucial for regulating smooth muscle tone, platelet activation, nociception and memory formation. As one of the main receptors for cGMP, PKGs mediate most of the effects of cGMP elevating drugs, such as nitric oxide-releasing agents and phosphodiesterase inhibitors which are used for the treatment of angina pectoris and erectile dysfunction, respectively.

Methodology/Principal Findings

We have investigated the mechanism of cyclic nucleotide binding to PKG by determining crystal structures of the amino-terminal cyclic nucleotide-binding domain (CNBD-A) of human PKG I bound to either cGMP or cAMP. We also determined the structure of CNBD-A in the absence of bound nucleotide. The crystal structures of CNBD-A with bound cAMP or cGMP reveal that cAMP binds in either syn or anti configurations whereas cGMP binds only in a syn configuration, with a conserved threonine residue anchoring both cyclic phosphate and guanine moieties. The structure of CNBD-A in the absence of bound cyclic nucleotide was similar to that of the cyclic nucleotide bound structures. Surprisingly, isothermal titration calorimetry experiments demonstrated that CNBD-A binds both cGMP and cAMP with a relatively high affinity, showing an approximately two-fold preference for cGMP.

Conclusions/Significance

Our findings suggest that CNBD-A binds cGMP in the syn conformation through its interaction with Thr193 and an unusual cis-peptide forming residues Leu172 and Cys173. Although these studies provide the first structural insights into cyclic nucleotide binding to PKG, our ITC results show only a two-fold preference for cGMP, indicating that other domains are required for the previously reported cyclic nucleotide selectivity.     

Prediction of the binding mode between GSK3β and a peptide derived from GSKIP using molecular dynamics simulation

GSK3β plays an important role in many physiological functions; dysregulated GSK3β is involved in human diseases such as diabetes, cancer, and Alzheimer's disease. This study uses MD simulations to determine the interaction between GSK3β and a peptide derived from GSKIP, a novel GSK3β interacting protein. Results show that GSKIPtide is inlaid in a binding pocket consisting of an α-helix and an extended loop near the carboxy-terminal end. This binding pocket is hydrophobic, and is responsible for the protein-protein interaction of two other GSK3β interacting proteins: FRAT and Axin. The GSKIPtide binding mode is closer to that of AxinGID (in the Axin-GSK3-interacting domain). The single-point mutations of V267G and Y288F in GSK3β differentiate the binding modes between GSK3 and GSKIPtide, AxinGID, and FRATide. The V2677G mutation of GSK3β reduces the GSKIPtide binding affinity by 70% and abolishes the binding affinity with AxinGID, but has no effect on FRATide. However, GSK3β Y288F completely abolishes the FRATide binding without affecting GSKIPtide or AxinGID binding. An analysis of the GSK3β-GSKIPtide complex structure and the X-ray crystal structures of GSK3β-FRATide and GSK3β-AxinGID complexes suggests that the hydroxyl group of Y288 is crucial to maintaining a hydrogen bond network in GSK3β-FRATide. The hydrophobic side chain of V267 maintains the integrity of helix-helix ridge-groove hydrophobic interaction for GSK3β-GSKIPtide and GSK3β-AxinGID. This study simulates these two mutant systems to provide atomic-level evidence of the aforementioned experimental results and validate the wild-type complex structure prediction.     

Crystal and molecular structures of substituted α-d-glucopyranosides

The crystal and molecular structures of methyl 2,4,6-tri-O-pivaloyl-α-d-glucopyranoside (1), methyl 4,6-O-(R)-benzylidene-2-O-pivaloyl-α-d-glucopyranoside (2), and methyl 4,6-O-(R)-benzylidene-2,3-di-O-pivaloyl-α-d-glucopyranoside (3) were determined by X-ray analysis. Crystals of 1 are orthorhombic, space group P2₁2₁2₁ with the unit cell a = 13.026(2), b = 16.832, c = 11.929(2) Å, Z = 4. Crystals of 2 are monoclinic, space group P2₁. The unit-cell parameters are a = 6.519(1), b = 14.664(4), c = 10.635(4) Å, β = 93.18(1)°, Z = 2. Crystals of 3 are orthorhombic, space group P2₁2₁2₁ with a = 10.006(3), b = 13.874(3), c = 18.527(5) Å, Z = 4. The structures were solved by MULTAN and refined by a full-matrix procedure to final values of R = 0.084 (1), 0.048 (2), and 0.069 (3). The pyranose ring in each compound adopts the ⁴C₁ conformation. The 1,3-dioxane rings in 2 and 3 show a chair conformation. The molecular packing in 1 is through the hydrogen bonds involving HO-3 and the 6-O-pivaloyl carbonyl group [HO-3 ⋯ O-9, 2.855(8) Å], which connect the molecules into a chain along

. The endocyclic oxygen atom is involved in an intermolecular hydrogen-bond with HO-3 [2.848(4) Å], joining molecules of 2 into the chains along

. There are no free hydroxyl groups in 3 and molecular packing reflects van der Waals interactions only.     

Abnormal actin binding of aberrant β-tropomyosins is a molecular cause of muscle weakness in TPM2-related nemaline and cap myopathy

NM (nemaline myopathy) is a rare genetic muscle disorder defined on the basis of muscle weakness and the presence of structural abnormalities in the muscle fibres, i.e. nemaline bodies. The related disorder cap myopathy is defined by cap-like structures located peripherally in the muscle fibres. Both disorders may be caused by mutations in the TPM2 gene encoding β-Tm (tropomyosin). Tm controls muscle contraction by inhibiting actin-myosin interaction in a calcium-sensitive manner. In the present study, we have investigated the pathogenetic mechanisms underlying five disease-causing mutations in Tm. We show that four of the mutations cause changes in affinity for actin, which may cause muscle weakness in these patients, whereas two show defective Ca2+ activation of contractility. We have also mapped the amino acids altered by the mutation to regions important for actin binding and note that two of the mutations cause altered protein conformation, which could account for impaired actin affinity.     

Enhancement of protonated molecular ions and ammonium·molecular ion complexes in direct-liquid-introduction-mass spectrometry of carbohydrates

Addition of ammonium acetate to the mobile phase in direct-liquid-introduction mass spectrometry enhances the abundance of the protonated molecular ion or ammonium·molecular ion complex for compounds of biological interest. The efficacy of the method was investigated by comparing mass spectra obtained, with and without ammonium acetate, for a variety of underivatized, per-O-acetylated, and per-O-alkylated carbohydrates, and for several underivatized peptides. The mass spectra of the per-O-alkylated carbohydrates obtained by direct-liquid-introduction mass spectrometry with ammonium acetate were also compared to those obtained by thermospray mass spectrometry.