GPU.proton.DOCK: Genuine Protein Ultrafast proton equilibria consistent DOCKing

GPU.proton.DOCK (Genuine Protein Ultrafast proton equilibria consistent DOCKing) is a state of the art service for in silico prediction of protein–protein interactions via rigorous and ultrafast docking code. It is unique in providing stringent account of electrostatic interactions self-consistency and proton equilibria mutual effects of docking partners. GPU.proton.DOCK is the first server offering such a crucial supplement to protein docking algorithms—a step toward more reliable and high accuracy docking results. The code (especially the Fast Fourier Transform bottleneck and electrostatic fields computation) is parallelized to run on a GPU supercomputer. The high performance will be of use for large-scale structural bioinformatics and systems biology projects, thus bridging physics of the interactions with analysis of molecular networks. We propose workflows for exploring in silico charge mutagenesis effects. Special emphasis is given to the interface-intuitive and user-friendly. The input is comprised of the atomic coordinate files in PDB format. The advanced user is provided with a special input section for addition of non-polypeptide charges, extra ionogenic groups with intrinsic pKa values or fixed ions. The output is comprised of docked complexes in PDB format as well as interactive visualization in a molecular viewer. GPU.proton.DOCK server can be accessed at http://gpudock.orgchm.bas.bg/

Driving calmodulin protein towards conformational shift by changing ionization states of select residues

Proteins are complex systems made up of many conformational sub-states which are mainly determined by the folded structure. External factors such as solvent type, temperature, pH and ionic strength play a very important role in the conformations sampled by proteins. Here we study the conformational multiplicity of calmodulin (CaM) which is a protein that plays an important role in calcium signaling pathways in the eukaryotic cells. CaM can bind to a variety of other proteins or small organic compounds, and mediates different physiological processes by activating various enzymes. Binding of calcium ions and proteins or small organic molecules to CaM induces large conformational changes that are distinct to each interacting partner. In particular, we discuss the effect of pH variation on the conformations of CaM. By using the pKa values of the charged residues as a basis to assign protonation states, the conformational changes induced in CaM by reducing the pH are studied by molecular dynamics simulations. Our current view suggests that at high pH, barrier crossing to the compact form is prevented by repulsive electrostatic interactions between the two lobes. At reduced pH, not only is barrier crossing facilitated by protonation of residues, but also conformations which are on average more compact are attained. The latter are in accordance with the fluorescence resonance energy transfer experiment results of other workers. The key events leading to the conformational change from the open to the compact conformation are (i) formation of a salt bridge between the N-lobe and the linker, stabilizing their relative motions, (ii) bending of the C-lobe towards the N-lobe, leading to a lowering of the interaction energy between the two-lobes, (iii) formation of a hydrophobic patch between the two lobes, further stabilizing the bent conformation by reducing the entropic cost of the compact form, (iv) sharing of a Ca+2 ion between the two lobes

PHEMTO: protein pH-dependent electric moment tools

PHEMTO (protein pH-dependent electric moment tools) is released in response to the high demand in protein science community for evaluation of electrostatic characteristics in relations to molecular recognition. PHEMTO will serve protein scientists with new advanced features for analysis of protein molecular interactions: Electric/dipole moments, their pH-dependence and in silico charge mutagenesis effects on these properties as well as alternative algorithms for electric/dipole moment computation—Singular value decomposition of electrostatic potential (EP) to account for reaction field. The implementation is based on long-term experience—PHEI mean field electrostatics and PHEPS server for evaluation of global and local pH-dependent properties. However, PHEMTO is not just an update of our PHEPS server. Besides standard electrostatics, we offer new, advanced and useful features for analysis of protein molecular interactions. In addition our algorithms are very fast. Special emphasis is given to the interface—intuitive and user-friendly. The input is comprised of the atomic coordinate file in Protein Data Bank format. The advanced user is provided with a special input section for addition of non-polypeptide charges. The output covers actually full electrostatic characteristics but special emphasis is given to electric/dipole moments and their interactive visualization. PHEMTO server can be accessed at http://phemto.orgchm.bas.bg/

Protonation States of Remote Residues Affect Binding-Release Dynamics of the Ligand but not the Conformation of apo Ferric Binding Protein

We have studied the apo (Fe3+ free) form of periplasmic ferric binding protein (FbpA) under different conditions and we have monitored the changes in the binding and release dynamics of H2PO4- that acts as a synergistic anion in the presence of Fe3+. Our simulations predict a dissociation constant of 2.2$\pm$0.2 mM which is in remarkable agreement with the experimentally measured value of 2.3$\pm$0.3 mM under the same ionization strength and pH conditions. We apply perturbations relevant for changes in environmental conditions as (i) different values of ionic strength (IS), and (ii) protonation of a group of residues to mimic a different pH environment. Local perturbations are also studied by protonation or mutation of a site distal to the binding region that is known to mechanically manipulate the hinge-like motions of FbpA. We find that while the average conformation of the protein is intact in all simulations, the H2PO4- dynamics may be substantially altered by the changing conditions. In particular, the bound fraction which is 20$\%$ for the wild type system is increased to 50$\%$ with a D52A mutation/protonation and further to over 90$\%$ at the protonation conditions mimicking those at pH 5.5. The change in the dynamics is traced to the altered electrostatic distribution on the surface of the protein which in turn affects hydrogen bonding patterns at the active site. The observations are quantified by rigorous free energy calculations. Our results lend clues as to how the environment versus single residue perturbations may be utilized for regulation of binding modes in hFbpA systems in the absence of conformational changes.Comment: 26 pages, 4 figure

First report on autochthonous urease-positive Trichophyton rubrum (T-raubitschekii) from South-east Europe

Trichophyton raubitschekii is a dermatophyte belonging to the T. rubrum complex and is differentiated principally by its positive urease activity and production of profuse macroconidia and microconidia in culture. It is classically isolated from African, South-east Asian and Australian aboriginal patients with tinea corporis or tinea cruris. This study was undertaken to screen Greek and Bulgarian clinical isolates identified as T. rubrum for T. raubitschekii and to delineate these strains by two molecular methods used for the first time in T. rubrum epidemiological studies. Ninety-five Greek and 10 Bulgarian strains, originating from various body sites, initially identified as T. rubrum, were screened for urease activity. The biochemical properties and morphology of the urease-positive strains were determined. Strains were delineated with polymerase chain reaction (PCR)-ribotyping amplifying repeat elements of the intergenic spacer region and by PCR fingerprinting. Five Greek and one Bulgarian T. raubitschekii strains were identified comprising isolates from patients with tinea manuum (one), tinea corporis (one), tinea cruris (one) and tinea unguium (three). Only one strain had the classical T. raubitschekii microscopic morphology, whereas the remaining five presented a dominant arthroconidial phenotype. Both typing methods clustered all T. raubitschekii and T. rubrum isolates together in the same group, indicating strain homogeneity in the genetic regions examined. The reported isolation of T. raubitschekii in the Balkan and South-eastern Mediterranean regions extends the geographical distribution of this species. As the more primitive T. raubitschekii probably represents the parental population of T. rubrum, the Greek and Bulgarian T. raubitschekii strains could represent a remnant of the T. rubrum spread that took place after the First World War, rather than being a recent epidemiological event

PHEPS: web-based pH-dependent Protein

PHEPS (pH-dependent Protein Electrostatics Server) is a web service for fast prediction and experiment planning support, as well as for correlation and analysis of experimentally obtained results, reflecting charge-dependent phenomena in globular proteins. Its implementation is based on long-term experience (PHEI package) and the need to explain measured physicochemical characteristics at the level of protein atomic structure. The approach is semi-empirical and based on a mean field scheme for description and evaluation of global and local pH-dependent electrostatic properties: protein proton binding; ionic sites proton population; free energy electrostatic term; ionic groups proton affinities (pKa,i) and their Coulomb interaction with whole charge multipole; electrostatic potential of whole molecule at fixed pH and pH-dependent local electrostatic potentials at user-defined set of points. The speed of calculation is based on fast determination of distance-dependent pair charge-charge interactions as empirical three exponential function that covers charge–charge, charge–dipole and dipole–dipole contributions. After atomic coordinates input, all standard parameters are used as defaults to facilitate non-experienced users. Special attention was given to interactive addition of non-polypeptide charges, extra ionizable groups with intrinsic pKas or fixed ions. The output information is given as plain-text, readable by ‘RasMol’, ‘Origin ’ and the like. The PHEPS server is accessible a