Remarkable disparity in mechanical response among the extracellular domains of type I and II cadherins

Cadherins, a large family of calcium-dependent adhesion molecules, are critical for intercellular adhesion. While crystallographic structures for several cadherins show clear structural similarities, their relevant adhesive strengths vary and their mechanisms of adhesion between types I and II cadherin subfamilies are still unclear. Here, stretching of cadherins was explored experimentally by atomic force microscopy and computationally by steered molecular dynamics (SMD) simulations, where partial unfolding of the E-cadherin ectodomains was observed. The SMD simulations on strand-swapping cadherin dimers displayed similarity in binding strength, suggesting contributions of other mechanisms to explain the strength differences of cell adhesion in vivo. Systematic simulations on the unfolding of the extracellular domains of type I and II cadherins revealed diverse pathways. However, at the earliest stage, a remarkable similarity in unfolding was observed for the various type I cadherins that was distinct from that for type II cadherins. This likely correlates positively with their distinct adhesive properties, suggesting that the initial forced deformation in type I cadherins may be involved in cadherin-mediated adhesion.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:25     

Unravelling novel congeners from acetyllysine mimicking ligand targeting a lysine acetyltransferase PCAF bromodomain

p300/CBP Associated Factor (PCAF) bromodomain (BRD), a lysine acetyltransferases, has emerged as a promising drug target as its dysfunction is linked to onset and progression of several diseases like cancer, diabetes, AIDS, etc. In this study, a three featured E-Pharmacophore (ARR) was generated based on acetyllysine mimicking inhibitor of PCAF BRD which is available as co-crystal structure (PDB ID: 5FDZ). It was used for filtering small molecule databases followed by molecular docking and consequently validated using enrichment calculation. The resulted hits were found to be congeners which show the predictive power of E-Pharmacophore hypothesis. Further, Induced Fit Docking method, Binding energy calculation, ADME prediction, Single Point Energy calculation and Molecular Dynamics simulation were performed to find better hits against PCAF BRD. Based on the results, it was concluded that Asn803, Tyr809 and Tyr802 along with a water molecule (HOH1001) plays crucial role in binding with inhibitor. It is also proposed that four hits from Life Chemicals database namely, F2276-0099, F2276-0008, F2276-0104 and F2276-0106 could act as potent drug molecules for PCAF BRD. Thus, the present study is strongly believed to have bright impact on rational drug design of potent and novel congeners of PCAF BRD inhibitors.     

Mechanism of flexibility control for ATP access of hepatitis C virus NS3 helicase

Hepatitis C virus (HCV) NS3 helicase couples adenosine triphosphate (ATP) binding and hydrolysis to polynucleotide unwinding. Understanding the regulation mechanism of ATP binding will facilitate targeting of the ATP-binding site for potential therapeutic development for hepatitis C. T324, an amino acid residue connecting domains 1 and 2 of NS3 helicase, has been suggested as part of a flexible hinge for opening of the ATP-binding cleft, although the detailed mechanism remains largely unclear. We used computational simulation to examine the mutational effect of T324 on the dynamics of the ATP-binding site. A mutant model, T324A, of the NS3 helicase apo structure was created and energy was minimized. Molecular dynamics simulation was conducted for both wild type and the T324A mutant apo structures to compare their differences. For the mutant structure, histogram analysis of pairwise distances between residues in domains 1 and 2 (E291-Q460, K210-R464 and R467-T212) showed that separation between the two domains was reduced by ～10% and the standard deviation by ～33%. Root mean square fluctuation (RMSF) analysis demonstrated that residues in close proximity to residue 324 have at least 30% RMSF value reductions in the mutant structure. Solvent RMSF analysis showed that more water molecules were trapped near D290 and H293 in domain 1 to form an extensive interaction network constraining cleft opening. We also demonstrated that the T324A mutation established a new atomic interaction with V331, revealing that an atomic interaction cascade from T324 to residues in domains 1 and 2 controls the flexibility of the ATP-binding cleft.     

Comparative study of the binding mode between cytochrome P450 17A1 and prostate cancer drugs in the absence of haem iron

Abstract

According to the X-ray crystal structures of CYP17A1 (including its complexes with inhibitors), it is shown that a hydrogen bond exists between CYP17A1 and its inhibitors (such as abiraterone and TOK-001). Previous short MD simulations (50?ns) suggested that the binding of abiraterone to CYP17A1 is stronger than that of TOK-001. In this work, by carrying out long atomistic MD simulations (200?ns) of CYP17A1 and its complexes with abiraterone and TOK-001, we observed a binding mode between CYP17A1 and abiraterone, which is different from the binding mode between CYP17A1 and TOK-001. In the case of abiraterone binding, the unfilled volume in the active site cavity increases the freedom of movement of abiraterone within CYP17A1, leading to the collective motions of the helices G and B′ as well as the breaking of hydrogen bond existing between the 3β-OH group of abiraterone and N202 of CYP17A1. However, the unfilled volume in the active site cavity can be occupied by the benzimidazole ring of TOK-001, restraining the motion of TOK-001. By pulling the two inhibitors (abiraterone and TOK-001) out of the binding pocket in CYP17A1, we discovered that abiraterone and TOK-001 were moved from their binding sites to the surface of protein similarly through the channels formed by the helices G and B′. In addition, based on the free energy calculations, one can see that it is energetically favorable for the two inhibitors (abiraterone and TOK-001) to enter into the binding pocket in CYP17A1.     

Characterization of the Structure and Melting Behavior of the Loop I fragment of ColE1 RNA I

Abstract

We have synthesized two RNA fragments: a 42-mer corresponding to the full loop I sequence of the loop I region of ColE1 antisense RNA (RNA I), plus three additional Gs at the 5′-end, and a 31-mer which has 11 5′-end nucleotides (G(-2)-U9) deleted. The secondary structure of the 42-mer, deduced from one- and two-dimensional NMR spectra, consists of a stem of 11 base-pairs which contains a U-U base-pair and a bulged C base, a 7 nucleotide loop, and a single-stranded 5′ end of 12 nucleotides. The UV-melting study of the 42-mer further revealed a multi-step melting behavior with transition temperatures 32°C and 71°C clearly discernible. In conjunction with NMR melting study the major transition at 71°C is assigned to the overall melting of the stem region and the 32°C transition is assigned to the opening of the loop region. The deduced secondary structure agrees with that proposed for the intact RNA I and provides structural bases for understanding the specificity of RNase E.     

Exploring the effect of D61G mutation on SHP2 cause gain of function activity by a molecular dynamics study

Noonan syndrome (NS) is a common autosomal dominant congenital disorder which could cause the congenital cardiopathy and cancer predisposition. Previous studies reported that the knock-in mouse models of the mutant D61G of SHP2 exhibited the major features of NS, which demonstrated that the mutation D61G of SHP2 could cause NS. To explore the effect of D61G mutation on SHP2 and explain the high activity of the mutant, molecular dynamic simulations were performed on wild type (WT) of SHP2 and the mutated SHP2-D61G, respectively. The principal component analysis and dynamic cross-correlation mapping, associated with secondary structure, showed that the D61G mutation affected the motions of two regions (residues Asn 58-Thr 59 and Val 460-His 462) in SHP2 from β to turn. Moreover, the residue interaction networks analysis, the hydrogen bond occupancy analysis and the binding free energies were calculated to gain detailed insight into the influence of the mutant D61G on the two regions, revealing that the major differences between SHP2-WT and SHP2-D61G were the different interactions between Gly 61 and Gly 462, Gly 61 and Ala 461, Gln 506 and Ile 463, Gly 61 and Asn 58, Ile 463 and Thr 466, Gly 462 and Cys 459. Consequently, our findings here may provide knowledge to understand the increased activity of SHP2 caused by the mutant D61G.     

Electron Microscopic and Physico-Chemical Studies of DNA Complexes with Synthetic Oligopeptides: Binding Specificity and DNA Compact Structures

Abstract

Binding to DNA of two synthetic peptides, Val-Thr-Thr-Val-Val-NH-NH-Dns and Thr-Val- Thr-Lys-Val-Gly-Thr-Lsy-Val-Gly-Thr-Val-Val-NH-NH-Dns (where Dns is a residue of 5- dimethylaminonaphthalene-l-sulfonic acid), has been studied by circular dichroism, electron microscopy and fluorescence methods. It has been found that these two peptides can self- associate in aqueous solution as follows from the fact that concentration-dependent changes are observed in the UV absorbance and fluorescence spectra. The two peptides can bind to DNA both in self-associated and monomeric forms. The pentapeptide in the β-associated form binds more strongly to poly(dG) · poly(dC) than to poly[d(A-C)] · poly[d(G-T)] and poly(dA) · poly(dT) whereas the tridecapeptide exhibits an opposite order of preferences binding more strongly to poly[d(A-C)] · poly[d(G-T)] and poly(dA) · poly(dT) than to poly(dG) · poly(dC).

Binding is a cooperative process which is accompanied by the DNA compaction at peptide/DNA base pair ratios greater than l. At the initial stage of the compaction process, the coalescence of DNA segments covered by bound peptide molecules leads to the formation of DNA loops stabilized by the interaction between peptide molecules bound to different DNA segments. Further increase in the peptide/DNA ratio leads to the formation of rod-like structures each consisting of two or more double-stranded DNA segments. The final stage of the compaction process involves folding of fibrillar macromolecular complexes into a globular structure containing only one DNA molecule.     

Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10?ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.     

Multiscale simulations of protein folding: application to formation of secondary structures

A multiscale simulation method of protein folding is proposed, using atomic representation of protein and solvent, combing genetic algorithms to determine the key protein structures from a global view, with molecular dynamic simulations to reveal the local folding pathways, thus providing an integrated landscape of protein folding. The method is found to be superior to previously investigated global search algorithms or dynamic simulations alone. For secondary structure formation of a selected peptide, RN24, the structures and dynamics produced by this method agree well with corresponding experimental results. Three most populated conformations are observed, including hairpin, β-sheet and α-helix. The energetic barriers separating these three structures are comparable to the kinetic energy of the atoms of the peptide, implying that the transition between these states can be easily triggered by kinetic perturbations, mainly through electrostatic interactions between charged atoms. Transitions between α-helix and β-sheet should jump over at least two energy barriers and may stay in the energetic trap of hairpin. It is proposed that the structure of proteins should be jointly governed by thermodynamic and dynamic factors; free energy is not the exclusive dominant for stability of proteins.