A superposition free method for protein conformational ensemble analyses and local clustering based on a differential geometry representation of backbone.

Here a differential geometry (DG) representation of protein backbone is explored on the analyses of protein conformational ensembles. The protein backbone is described by curvature, κ, and torsion, τ, values per residue and we propose 1) a new dissimilarity and protein flexibility measurement and 2) a local conformational clustering method. The methods were applied to Ubiquitin and c-Myb-KIX protein conformational ensembles and results show that κ\τ metric space allows to properly judge protein flexibility by avoiding the superposition problem. The dmax measurement presents equally good or superior results when compared to RMSF, especially for the intrinsically unstructured protein. The clustering method is unique as it relates protein global to local dynamics by providing a global clustering solutions per residue. The methods proposed can be especially useful to the analyses of highly flexible proteins. The software written for the analyses presented here is available at https://github.com/AMarinhoSN/FleXgeo for academic usage only.

The role of local residue environmental changes in thermostable mutants of the GH11 xylanase from Bacillus subtilis.

A thermostable variant of the mesophilic xylanase A from Bacillus subtilis (BsXynA-G3_4x) contains the four mutations Gln7His, Gly13Arg, Ser22Pro, and Ser179Cys. The crystal structure of the BsXynA-G3_4x has been solved, and the local environments around each of these positions investigated by molecular dynamics (MD) simulations at 328K and 348K. The structural and MD simulation results were correlated with thermodynamic data of the wild-type enzyme, the 4 single mutants and the BsXynA-G3_4x. This analysis suggests that the overall stabilizing effect is entropic, and is consistent with solvation of charged residues and reduction of main-chain flexibility. Furthermore, increased protein-protein hydrogen bonding and hydrophobic interactions also contribute to stabilize the BsXynA-G3_4x. The study revealed that a combination of several factors is responsible for increased thermostability of the BsXynA-G3_4x; (i) introduction of backbone rigidity in regions of high flexibility, (ii) solvation effects and (iii) hydrophobic contacts.

Characterization of suramin binding sites on the human group IIA secreted phospholipase A2 by site-directed mutagenesis and molecular dynamics simulation.

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the hsPLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors.