The conformation of the Congo-red ligand bound to amyloid fibrils HET-s(218–289): a solid-state NMR study

We have previously shown that Congo red (CR) binds site specifically to amyloid fibrils formed by HET-s(218–289) with the long axis of the CR molecule almost parallel to the fibril axis. HADDOCK docking studies indicated that CR adopts a roughly planar conformation with the torsion angle ? characterizing the relative orientation of the two phenyl rings being a few degrees. In this study, we experimentally determine the torsion angle ? at the center of the CR molecule when bound to HET-s(218–289) amyloid fibrils using solid-state NMR tensor-correlation experiments. The method described here relies on the site-specific ¹³C labeling of CR and on the analysis of the two-dimensional magic-angle spinning tensor-correlation spectrum of ¹³C₂-CR. We determined the torsion angle ? to be 19°.     

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Insights into the structural determinants for selective inhibition of nitric oxide synthase isoforms

Selective inhibition of the nitric oxide synthase isoforms (NOS) is a promising approach for the treatment of various disorders. However, given the high active site conservation among all NOS isoforms, the design of selective inhibitors is a challenging task. Analysis of the X-ray crystal structures of the NOS isoforms complexed with known inhibitors most often gives no clues about the structural determinants behind the selective inhibition since the inhibitors share the same binding conformation. Aimed at a better understanding of the structural factors responsible for selective inhibition of NOS isoforms we have performed MD simulations for iNOS, nNOS and eNOS complexed with N^ω-NO₂-L-Arg (1), and with the aminopyridine derivatives 2 and 3. The slightly better selectivity of 1 for nNOS may be assigned to the presence of extra charge–charge interactions due to its “extended” conformation. While the high affinity of 2 for iNOS can be explained by the formation of an iNOS-specific subpocket upon binding, the lack of affinity for eNOS is associated to a conformational change in Glu363. The strong van der Waals and electrostatic interactions between 3 and the active site of nNOS are most likely responsible for its higher affinity for this isoform. Owing to the elongated and narrow binding pocket of iNOS, the correct positioning of 3 over the heme group is difficult, which may account for its lower affinity toward this isoform. Brought together, our results might help to rationalize the design of selective NOS inhibitors.

Computational design of a full-length model of HIV-1 integrase: modeling of new inhibitors and comparison of their calculated binding energies with those previously studied

A full-length model of integrase (IN) of the human immunodeficiency virus type 1 (HIV-1) was constructed based on the distinctly resolved X-ray crystal structures of its three domains, named N-terminal, catalytic core and C-terminal. Thirty-one already known inhibitors with varieties of structural differences as well as nine newly tested ones were docked into the catalytic core. The molecular dynamic (MD) and binding properties of these complexes were obtained by MD calculations. The binding energies calculated by molecular mechanic/Poisson Boltzmann solvation area were significantly correlationed with available IC₅₀. Four inhibitors including two newly designed were also docked into the full-length model and their MD behaviors and binding properties were calculated. It was found that one of the newly designed compounds forms a better complex with HIV-1 IN compared to the rest including raltegravir. MD calculations were performed with AMBER suite of programs using ff99SB force field for the proteins and the general Amber force field for the ligands. In conclusion, the results have produced a promising standpoint not only in the construction of the full-length model but also in development of new drugs against it. However, the role of multimer formation and the involvement of DNAs, and their subsequent effect on the complexation and inhibition, are required to arrive at a conclusive decision.

Amyloid associated with elastin-staining laminar aggregates in the lungs of patients diagnosed with acute respiratory distress syndrome

Background

The heterogeneity of conditions underlying respiratory distress, whether classified clinically as acute lung injury (ALI) or the more severe acute respiratory distress syndrome (ARDS), has hampered efforts to identify and more successfully treat these patients. Examination of postmortem lungs among cases clinically diagnosed as ARDS identified a cohort that showed a consistent morphology at the light and electron microscope levels, and featured pathognomonic structures which we termed elastin-staining laminar structures (ELS).

Methods

Postmortem tissues were stained using the Verhoeff-Van Gieson procedure for elastic fibers, and with Congo red for examination under a polarizing microscope. Similar samples were examined by transmission EM.

Results

The pathognomonic ELS presented as ordered molecular aggregates when stained using the Verhoeff-van Gieson technique for elastic fibers. In several postmortem lungs, the ELS also displayed apple-green birefringence after staining with Congo red, suggesting the presence of amyloid. Remarkably, most of the postmortem lungs with ELS exhibited no significant acute inflammatory cellular response such as neutrophilic reaction, and little evidence of widespread edema except for focal intra-alveolar hemorrhage.

Conclusions

Postmortem lungs that exhibit the ELS constitute a morphologically-identifiable subgroup of ARDS cases. The ordered nature of the ELS, as indicated by both elastin and amyloid stains, together with little morphological evidence of inflammation or edema, suggests that this cohort of ARDS may represent another form of conformational disease. If this hypothesis is confirmed, it will require a new approach in the diagnosis and treatment of patients who exhibit this form of acute lung injury.     

Molecular dynamics and free energy studies of chirality specificity effects on aminobenzo[a]quinolizine inhibitors binding to DPP-IV

The aminobenzo[a]quinolizines were investigated as a novel class of DPP-IV inhibitors. The stereochemistry of this class plays an important role in the bioactivity. In this study, the mechanisms of how different configuration of three chiral centers of this class influences the binding affinity were investigated by molecular dynamics simulations, free energy decomposition analysis. The S configuration for chiral center 3* is decisive for isomers to maintain high bioactivity; the chirality effect of chiral center 2* on the binding affinity is largely dependent, while the S configuration for chiral center 2* is preferable to R configuration for the bioactivity gain; the effect of chiral center 11b* on the binding affinity is insignificant. The chirality specificity for three chiral centers is responsible for distinction of two van der Waals contacts with Tyr547 and Phe357, and of H-bonding interactions with Arg125 and Glu206. Particularly, the Arg125 to act as a bridge in the H-bonding network contributes to stable H-bonding interactions of isomer in DPP-IV active site.

Molecular dynamics analysis of a series of 22 potential farnesyltransferase substrates containing a CaaX-motif

Protein farnesyltransferase (FTase) is an important target in many research fields, more markedly so in cancer investigation since several proteins known to be involved in human cancer development are thought to serve as substrates for FTase and to require farnesylation for proper biological activity. Several FTase inhibitors (FTIs) have advanced into clinical testing. Nevertheless, despite the progress in the field several functional and mechanistic doubts on the FTase catalytic activity have persisted. This work provides some crucial information on this important enzyme by describing the application of molecular dynamics simulations using specifically designed molecular mechanical parameters for a variety of 22 CaaX peptides known to work as natural substrates or inhibitors for this enzyme. The study involves a comparative analysis of several important molecular aspects, at the mechanistic level, of the behavior of substrates and inhibitors at the dynamic level, including the behavior of the enzyme and peptides, as well as their interaction, together with the effect of the solvent. Properties evaluated include the radial distribution function of the water molecules around the catalytically important zinc metal atom and cysteine sulfur of CaaX, the conformations of the substrate and inhibitor and the corresponding RMSF values, critical hydrogen bonds, and several catalytically relevant distances. These results are discussed in light of recent experimental and computational evidence that provides new insights into the activity of this enzyme.

Discovery of potent inhibitors for interleukin-2-inducible T-cell kinase: structure-based virtual screening and molecular dynamics simulation approaches

In our study, a structure-based virtual screening study was conducted to identify potent ITK inhibitors, as ITK is considered to play an important role in the treatment of inflammatory diseases. We developed a structure-based pharmacophore model using the crystal structure (PDB ID: 3MJ2) of ITK complexed with BMS-50944. The most predictive model, SB-Hypo1, consisted of six features: three hydrogen-bond acceptors (HBA), one hydrogen-bond donor (HBD), one ring aromatic (RA), and one hydrophobic (HY). The statistical significance of SB-Hypo1 was validated using wide range of test set molecules and a decoy set. The resulting well-validated model could then be confidently used as a 3D query to screen for drug-like molecules in a database, in order to retrieve new chemical scaffolds that may be potent ITK inhibitors. The hits retrieved from this search were filtered based on the maximum fit value, drug-likeness, and ADMET properties, and the hits that were retained were used in a molecular docking study to find the binding mode and molecular interactions with crucial residues at the active site of the protein. These hits were then fed into a molecular dynamics simulation to study the flexibility of the activation loop of ITK upon ligand binding. This combination of methodologies is a valuable tool for identifying structurally diverse molecules with desired biological activities, and for designing new classes of selective ITK inhibitors.

Rational design of the survivin/CDK4 complex by combining protein–protein docking and molecular dynamics simulations

Survivin, the smallest inhibitor of apoptosis protein (IAP), is a valid target for cancer research. It mediates both the apoptosis pathway and the cell cycle and has been proposed to form a complex with the cyclin-dependent kinase protein CDK4. The resulting complex transports CDK4 from the cytosol to the nucleus, where CDK4 participates in cell division. Survivin has been recognized as a node protein that interacts with several partners; disruption of the formed complexes can lead to new anticancer compounds. We propose a rational model of the survivin/CDK4 complex that fulfills the experimental evidence and that can be used for structure-based design of inhibitors modifying its interface recognition. In particular, the suggested complex involves the alpha helical domain of survivin and resembles the mode of binding of survivin in the survivin/borealin X-ray structure. The proposed model has been obtained by combining protein–protein docking, fractal-based shape complementarity, electrostatics studies and extensive molecular dynamics simulations.

DFT study of the Lewis acid mediated synthesis of 3-acyltetramic acids

The recent emergence and re-emergence of alphaviruses, in particular the chikungunya virus (CHIKV), in numerous countries has invoked a worldwide threat to human health, while simultaneously generating an economic burden on affected countries. There are currently no vaccines or effective drugs available for the treatment of the CHIKV, and with few lead compounds reported, the vital medicinal chemistry is significantly more challenging. This study reports on the discovery of potential inhibitors for the nsP3 macro domain of CHIKV using molecular docking, virtual screening, and molecular dynamics simulations, as well as work done to evaluate and confirm the active site of nsP3. Virtual screening was carried out based on blind docking as well as focused docking, using the database of 1541 compounds from NCI Diversity Set II, to identify hit compounds for nsP3. The top hit compounds were further subjected to molecular dynamic simulations, yielding a greater understanding of the dynamic behavior of nsP3 and its complexes with various ligands, concurrently confirming the outcomes of docking, and establishing in silico lead compounds which target the CHIKV nsP3 enzyme.

Cytotoxic effect and molecular docking of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide—a novel topoisomerase II inhibitor

The preliminary cytotoxic effect of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide hydrochloride (1)—a potent topoisomerase II inhibitor—was measured using a MTT assay. It was found that the compound decreased the number of viable cells in both estrogen receptor-positive MCF-7 and estrogen receptor-negative MDA-MB-231breast cancer cells, with IC₅₀ values of 146?±?2 and 132?±?2 μM, respectively. To clarify the molecular basis of the inhibitory action of 1, molecular docking studies were carried out. The results suggest that 1 targets the ATP binding pocket.

Binding mode of Thioflavin T and other molecular probes in the context of amyloid fibrils—current status

Because understanding amyloid fibrillation in molecular detail is essential for development of strategies to control amyloid formation and overcome neurodegenerative disorders, increased understanding of present molecular probes as well as development of new probes are of utmost importance. To date, the binding modes of these molecular probes to amyloid fibrils are by no means adequately described or understood, and the large number of studies on Thioflavin T (ThT) and Congo Red (CR) binding have resulted in models that are incomplete and conflicting. Different types of binding sites are likely to be present in amyloid fibrils with differences in binding modes. ThT may bind in channels running parallel to the long axis of the fibril. In the channels, ThT may bind in either a monomeric or dimeric form of which the molecular conformation is likely to be planar. CR may bind in grooves formed along the β-sheets as a planar molecule in either a monomeric or supramolecular form.     

Effects of protein binding on a lipid bilayer containing local anesthetic articaine, and the potential of mean force calculation: a molecular dynamics simulation approach

Articaine, as a local anesthetic drug has been simulated in neutral and charged forms, and its interaction with the dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane is investigated by molecular dynamics simulation using GROMACS software. In order to obtain the optimum location of the drug molecules, as they penetrate into the membrane, umbrella sampling is applied and the free energy is calculated. The effect of protein binding to DMPC membrane on the process of drug diffusion through the membrane is considered. Five simulation systems are designed and by applying the potential of mean force, the molecular dynamics simulation on the system is performed. In light of the obtained results, the electrostatic potential, variation of lipid bilayer’s order parameter and the diffusion coefficient of drug are discussed.

Retrospective molecular docking study of WY-25105 ligand to β-secretase and bias of the three-dimensional structure flexibility

β-Secretase (BACE) is a very promising target in the search for a treatment for Alzheimer’s disease using a protein–ligand inhibition approach. Given the many published X-ray structures of BACE protein, structure-based drug design has been used extensively to support new inhibitor discovery programs. Due to the high flexibility and large catalytic site of this protein, sampling of the huge conformational space of the binding site is the big challenge to overcome and is the main limitation of the most widely used docking programs. Incorrect treatment of these pitfalls can introduce bias into ligand docking and could affect the results. This is especially the case with the WY-25105 compound reported by the Wyeth Corporation as a BACE ligand that did not fit into any of the known crystal structures. In the present retrospective study, a set of available X-ray enzyme structures was selected and molecular dynamics simulations were conducted to generate more diverse representative BACE protein conformations. These conformations were then used for a docking study of the WY-25105 compound. The results confirmed the need to use an ensemble of structures in protein–ligand docking for identification of new binding modes in structure-based drug design of BACE inhibitors.

Investigation of the binding network of IGF-I on the cavity surface of IGFBP4

Insulin-like growth factor-binding proteins (IGFBPs) control bioactivity and distribution of insulin-like growth factors (IGFs) through high-affinity complex of IGFBP and IGF. To get more insight into the binding interaction of IGF system, the site-directed mutagenesis and force-driving desorption methods were employed to study the interaction mechanism of IGFBP4 and IGF-I by molecular dynamics (MD) simulation. In IGF-I, residues Gly7 to Asp12 were found to be the hot spots and they mainly anchored on the N-domain of IGFBP4. The contact area, the shape and size of protein, the surroundings of the binding site, the hydrophobic and electrostatic interaction between the two proteins worked as a complex network to regulate the protein-protein interaction. It was also found that the unfolding of the helix was not inevitable in the mutant, and it could be regulated by careful selection of the substituted amino acid.

Effect of novel synthetic evodiamine analogue on sexual behavior in male rats

Currently phosphodiestrase5 (PDE5) inhibitors are the first-line treatment for erectile dysfunction. Drugs such as sildenafil and tadalafil are available as PDE5 inhibitors which are potent and reversible but lack selectivity with side effects such as headache, facial flushing, dyspepsia, and visual disturbances. We herein report for the first time novel condensed thienopyrimidines as evodiamine analogue and their effect on sexual behavior in male rats hitherto unreported. Novel synthetic evodiamine significantly showed improvement in male rat copulatory behavior. The test compound MKAC9 could be of promising importance in the treatment of sexual disorders like desire disorder or erectile dysfunction.

A plausible explanation for enhanced bioavailability of P-gp substrates in presence of piperine: simulation for next generation of P-gp inhibitors

P-glycoprotein (P-gp) has a major role to play in drug pharmacokinetics and pharmacodynamics, since it effluxes many cytotoxic hydrophobic anticancer drugs from gastrointestinal tract, brain, liver and kidney. Piperine is known to enhance the bioavailability of curcumin, as a substrate of P-gp by at least 2000 %. Besides these at least 50 other substrates and inhibitors of P-gp have been reported so far. All P-gp inhibitors have diverse structures. Although little is known about binding of some flavonoids and steroids at the NBD (nucleotide binding domain) of P-gp in the vicinity of ATP binding site inhibiting its hydrolysis, a valid explanation of how P-gp accommodates such a diverse set of inhibitors is still awaited. In the present study, piperine up to 100 μM has not shown observable cytotoxic effect on MDCK cell line, and it has been shown to accumulate rhodamine by fluorescence microscopy and fluorescent activated cell sorter in MDCK cells. Computational simulation for piperine and some first and second generation P-gp inhibitors has shown that these dock at the NBD site of P-gp. A comparative simulation study has been carried out regarding their docking and binding energies. Binding conformation of P-gp co-crystallized complexes with ADP, AMP-PNP (Adenylyl-imidodiphosphate), and ATP were compared with piperine. The receptor based E-pharmacophore of docked piperine has been simulated to find common features amongst P-gp inhibitors. Finally it has been concluded that piperine could be utilized as base molecule for design and development of safe non-toxic inhibitor of P-gp in order to enhance the bioavailability of most of its substrates.

Determination of key receptor–ligand interactions of dopaminergic arylpiperazines and the dopamine D2 receptor homology model

Interest in structure-based G-protein-coupled receptor (GPCR) ligand discovery is huge, given that almost 30 % of all approved drugs belong to this category of active compounds. The GPCR family includes the dopamine receptor subtype D2 (D2DR), but unfortunately—as is true of most GPCRs—no experimental structures are available for these receptors. In this publication, we present the molecular model of D2DR based on the previously published crystal structure of the dopamine D3 receptor (D3DR). A molecular modeling study using homology modeling and docking simulation provided a rational explanation for the behavior of the arylpiperazine ligand. The observed binding modes and receptor–ligand interactions provided us with fresh clues about how to optimize selectivity for D2DR receptors.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

Theoretical study of the solvatochromism of a donor-acceptor bithiophene

The solvation and the solvatochromic behavior of the 5-(methylthio)-5′-nitro-2,2′-bithiophene 1 in diethyl ether, dichloromethane, acetonitrile, methanol and formamide was theoretically investigated with an iterative molecular and quantum mechanics (QM/MM) approach. Calculated longest-wavelength solvatochromic absorption band of 1, obtained as averages of statistically uncorrelated configurations, including the solute and explicit solvent molecules of the first and second solvation layer, were in excellent agreement with the experimental results.