Recognition of ribonuclease A by 3'-5'-pyrophosphate-linked dinucleotide inhibitors: a molecular dynamics/continuum electrostatics analysis

Mechanistic studies on substrate inhibition and substrate activation of ribonuclease A: experimental and in silico investigation

Ribonuclease A (RNase A) is an endonuclease that plays a significant role in antimicrobial activity by the cleavage and hydrolysis of ssRNA. In this study, the interactions between RNase A and cCMP have been investigated, through molecular docking and a 200 ns molecular dynamics simulation. The enzyme kinetic properties were analyzed using saturation curve, Eadie-Hofstee, and Hill plots. The docking results indicate that the cCMP-RNase A complexes are stabilized through hydrophobic interaction, hydrogen bonding, and π-π stacking interaction. The most binding site is observed in the catalytic cleft, especially at residue His12 and His119. Enzyme-ligand docking study indicates that cCMP binds to residues located in the internal cavity of the catalytic site and shows three phases of binding modes. The analysis of MD simulations shows that RMSD, Rg, and RMSF reach equilibrium. The ΔGbinding of the cCMP-RNase A was negative (-619.673 kJ/mol), The comparison between the results pre and post MD simulation showed that ΔGbinding after MD simulation causes to more stable the structure than before simulation. Experimental methods such as saturation, Eadie-Hofstee, and Hill plots confirm theoretical data and show three phases of binding modes as well. Two significant events are demonstrated in the interaction between RNase A and cCMP: substrate activation and substrate inhibition are observed in concentrations below and above 0.8 mM, respectively, for cCMP. Choosing the appropriate concentration of cCMP is very important in investigation of ribonuclease A's catalytic behavour, especially for exploration of the effects of some drugs on biological behaviours related to ribonuclease A.Communicated by Ramaswamy H. Sarma.

Design of configuration-restricted triazolylated β-d-ribofuranosides: a unique family of crescent-shaped RNase A inhibitors

Seven new carbohydrate-bistriazole hybrid molecules were designed taking into consideration the crescent shaped active site of ribonuclease A (RNase A). In this case, the β-d-ribofuranose structure was used as the basic building unit; both the C1 and C4 arms protruding out towards the β-face of the tetrahydrofuran moiety of the ribose sugar provided an overall "U" shape to the basic building block. Several combinations of bistriazole moieties were constructed on the two arms of this basic building block. These mono- and/or bis-substituted 1,2,3-triazole units were linked to acidic functional groups because of the overall basic nature of the hydrolytic site of RNase A. All these compounds were efficient competitive inhibitors of RNase A with inhibition constants (Ki) in the micromolar range. In contrast to the carboxylic acid-modified hybrid molecules, molecules carrying sulfonic acids were found to be more potent because of the stronger interactions with the positively charged active site. The most efficient inhibitor of the series was the disulfonic acid-functionalized carbohydrate-bis-triazole hybrid molecule. Docking studies disclosed that the molecule, because of its well defined "U" shape with flexible arms, fits effectively in the active site; moreover, in all cases, besides the acid groups, the triazole and sugar rings also actively participated in creating the hydrogen bonding network in the cavity of the enzyme active site.

The ammonium sulfate inhibition of human angiogenin

In this study, we investigate the inhibition of human angiogenin by ammonium sulfate. The inhibitory potency of ammonium sulfate for human angiogenin (IC50 = 123.5 ± 14.9 mm) is comparable to that previously reported for RNase A (119.0 ± 6.5 mm) and RNase 2 (95.7 ± 9.3 mm). However, analysis of two X-ray crystal structures of human angiogenin in complex with sulfate anions (in acidic and basic pH environments, respectively) indicates an entirely distinct mechanism of inhibition. While ammonium sulfate inhibits the ribonucleolytic activity of RNase A and RNase 2 by binding to the active site of these enzymes, sulfate anions bind only to peripheral substrate anion-binding subsites of human angiogenin, and not to the active site.

Computational protein design with a generalized Born solvent model: application to Asparaginyl-tRNA synthetase

Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization.

Functional and structural analyses of N-acylsulfonamide-linked dinucleoside inhibitors of RNase A

Molecular probes are useful for both studying and controlling the functions of enzymes and other proteins. The most useful probes have high affinity for their target, along with small size and resistance to degradation. Here, we report on new surrogates for nucleic acids that fulfill these criteria. Isosteres in which phosphoryl [R–O–P(O₂⁻)–O–R′] groups are replaced with N-acylsulfonamidyl [R–C(O)–N⁻–S(O₂)–R′] or sulfonimidyl [R–S(O₂)–N⁻–S(O₂)–R′] groups increase the number of nonbridging oxygens from two (phosphoryl) to three (N-acylsulfonamidyl) or four (sulfonimidyl). Six such isosteres were found to be more potent inhibitors of catalysis by bovine pancreatic RNase A than are parent compounds containing phosphoryl groups. The atomic structures of two RNase A·N-acylsulfonamide complexes were determined at high resolution by X-ray crystallography. The N-acylsulfonamidyl groups were observed to form more hydrogen bonds with active site residues than did the phosphoryl groups in analogous complexes. These data encourage the further development and use of N-acylsulfonamides and sulfonimides as antagonists of nucleic acid-binding proteins.

Kinetics, in silico docking, molecular dynamics, and MM-GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP-binding site inhibitors: the role of water molecules examined

With an aim toward glycogenolysis control in Type 2 diabetes, we have investigated via kinetic experiments and computation the potential of indirubin (IC₅₀ > 50 μM), indirubin-3'-oxime (IC₅₀ = 144 nM), KT5720 (K(i) = 18.4 nM) and staurosporine (K(i) = 0.37 nM) as phosphorylase kinase (PhKγtrnc) ATP-binding site inhibitors, with the latter two revealed as potent inhibitors in the low nM range. Because of lack of structural information, we have exploited information from homologous kinase complexes to direct in silico calculations (docking, molecular dynamics, and MMGBSA) to predict the binding characteristics of the four ligands. All inhibitors are predicted to bind in the same active site area as the ATP adenine ring, with binding dominated by hinge region hydrogen bonds to Asp104:O and Met106:O (all four ligands) and also Met106:NH (for the indirubins). The PhKγtrnc-staurosporine complex has the greatest number of receptor-ligand hydrogen bonds, while for the indirubin-3'-oxime and KT5720 complexes there is an important network of interchanging water molecules bridging inhibitor-enzyme contacts. The MM-GBSA results revealed the source of staurosporine's low nM potency to be favorable electrostatic interactions, while KT5720 has strong van der Waals contributions. KT5720 interacts with the greatest number of protein residues either by direct or 1-water bridged hydrogen bond interactions, and the potential for more selective PhK inhibition based on a KT5720 analogue has been established. Including receptor flexibility in Schrödinger induced-fit docking calculations in most cases correctly predicted the binding modes as compared with the molecular dynamics structures; the algorithm was less effective when there were key structural waters bridging receptor-ligand contacts.

Species specificity of the complement inhibitor compstatin investigated by all-atom molecular dynamics simulations

The development of compounds to regulate the activation of the complement system in non-primate species is of profound interest because it can provide models for human diseases. The peptide compstatin inhibits protein C3 in primate mammals and is a potential therapeutic agent against unregulated activation of complement in humans but is inactive against nonprimate species. Here, we elucidate this species specificity of compstatin by molecular dynamics simulations of complexes between the most potent natural compstatin analog and human or rat C3. The results are compared against an experimental conformation of the human complex, determined recently by X-ray diffraction at 2.4-A resolution. The human complex simulations provide information on the relative contributions to stability of specific C3 and compstatin residues. In the rat simulations, the protein undergoes reproducible conformational changes, which eliminate or weaken specific interactions and reduce the complex stability. The simulation insights can be used to design improved compstatin-based inhibitors for human C3 and active inhibitors against lower mammals.