New software for finding transition states by probing accessible,or ergodic,regions

An automated, iterative approach to finding the lowest energy, ionic diffusion paths through a periodic structure has been developed within our new code (written in FORTRAN 77 and named Bubble). The approach is quite general in that it can be applied to find, at a chosen temperature, the accessible (ergodic) regions of a hyper-surface, which is defined across a uniform grid [1 Schön, J.C., Putz, H. and Jansen, M. 1996. Studying the energy hypersurface of continuous system—The threshold algorithm. J. Phys.-Conden. Matt., 8: 143[Crossref] , [Google Scholar]]. We describe both our implementation within the Bubble code and its application to locating the approximate transition states for Mg interstitial diffusion in forsterite, which can then be refined using standard transition state searching [2 Banerjee, A., Adams, N., Simons, J. and Shepard, R. 1985. Search for stationary points on surfaces. J. Phys. Chem., 89: 52[Crossref], [Web of Science ®] , [Google Scholar]].     

Alpha7 Helix Plays an Important Role in the Conformational Stability of PTP1B

Abstract

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix (α7), which forms an allosteric binding site 20 Å away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)- α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Molecular Mechanics Calculations of Proteins. Comparison of Different Energy Minimization Strategies

Abstract

A general strategy for performing energy minimization of proteins using the SYBYL molecular modelling program has been developed. The influence of several variables including energy minimization procedure, solvation, dielectric function and dielectric constant have been investigated in order to develop a general method, which is capable of producing high quality protein structures. Avian pancreatic polypeptide (APP) and bovine pancreatic phospholipase A2 (BP PLA2) were selected for the calculations, because high quality X-ray structures exist and because all classes of secondary structure are represented in the structures. The energy minimized structures were evaluated relative to the corresponding X-ray structures. The overall similarity was checked by calculating RMS distances for all atom positions. Backbone conformation was checked by Ramachandran plots and secondary structure elements evaluated by the length on hydrogen bonds. The dimensions of active site in BP PLA2 is very dependent on electrostatic interactions, due to the presence of the positively charged calcium ion. Thus, the distances between calcium and the calcium-coordinating groups were used as a quality index for this protein. Energy minimized structures of the trimeric PLA2 from Indian cobra (N.n.n. PLA2) were used for assessing the impact of protein-protein interactions. Based on the above mentioned criteria, it could be concluded that using the following conditions: Dielectric constant ? = 4 or 20; a distance dependent dielectric function and stepwise energy minimization, it is possible to reproduce X-ray structures very accurately without including explicit solvent molecules.     

An N-terminal, 830 residues intrinsically disordered region of the cytoskeleton-regulatory protein supervillin contains Myosin II- and F-actin-binding sites

Supervillin, the largest member of the villin/gelsolin family, is a cytoskeleton regulating, peripheral membrane protein. Supervillin increases cell motility and promotes invasive activity in tumors. Major cytoskeletal interactors, including filamentous actin and myosin II, bind within the unique supervillin amino terminus, amino acids 1–830. The structural features of this key region of the supervillin polypeptide are unknown. Here, we utilize circular dichroism and bioinformatics sequence analysis to demonstrate that the N-terminal part of supervillin forms an extended intrinsically disordered region (IDR). Our combined data indicate that the N-terminus of human and bovine supervillin sequences (positions 1–830) represents an IDR, which is the largest IDR known to date in the villin/gelsolin family. Moreover, this result suggests a potentially novel mechanism of regulation of myosin II and F-actin via the intrinsically disordered N-terminal region of hub protein supervillin.     

Structural analysis of the envelope gp120 V3 loop for some HIV-1 variants circulating in the countries of Eastern Europe

The V3 loop on gp120 from human immunodeficiency virus type 1 (HIV-1) is a focus of many research groups involved in anti-AIDS drug development because this region of the protein is a principal target for neutralizing antibodies and a major determinant for cell tropism and syncytium formation. In this study, the nucleotide sequences of the env gene region coding the V3 loop were determined by DNA sequencing methods for four novel HIV-1 strains that circulate in the countries of Eastern Europe, such as Russia, Belarus, Ukraine, etc. Based on the empirical data obtained, the 3D structures of the V3 loops associated with these viral modifications were generated by computer modeling and then compared to discover similarities in the spatial arrangement of this functionally important site of gp120. Despite the HIV-1 genetic variety, several regions of the V3 loop that contain residues critical for cell tropism were shown to be structurally invariant, which may explain its exceptional role in a co-receptor usage. These data together with those on the biological activity of the V3 individual residues clearly show that these conserved structural motifs of gp120 represent potential HIV-1 weak points most suitable for therapeutic intervention.     

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800?K, for 25?ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500?K has several conformations reaching into the protein core, whereas at 600?K some of the alpha-helices were lost. At 800?K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.     

Binding Sites of the Polyamines Putrescine,Cadaverine, Spermidine and Spermine on A- and B-DNA Located by Simulated Annealing

Abstract

Molecular dynamics simulations with simulated annealing are performed on polyamine-DNA systems in order to determine the binding sites of putrescine, cadaverine, spermidine and spermine on A- and B-DNA. The simulations either contain no additional counterions or sufficient Na⁺ ions, together with the charge on the polyamine, to provide 73% neutralisation of the charges on the DNA phosphates. The stabilisation energies of the complexes indicate that all four polyamines should stabilise A-DNA in preference to B-DNA, which is in agreement with experiment in the case of spermine and spermidine, but not in the case of putrescine or cadaverine. The major groove is the preferred binding site on A-DNA of all the polyamines. Putrescine and cadaverine tend to bind to the sugar-phosphate backbone of B-DNA, whereas spermidine and spermine occupy more varied sites, including binding along the backbone and bridging both the major and minor grooves.     

Genetic structure of Barbus spp. populations in the Marches Region of central Italy and its relevance to conservation actions

A genetic survey of Barbus spp. populations in the Marches Region (Adriatic River basins), central Italy, was carried out using mitochondrial and nuclear markers (partial D‐loop, cyt b sequences and microsatellite loci) in order to ascertain their systematic position and to address their genetic structure which is key to conservation action planning. Analyses were conducted on sequences obtained from 91 individuals collected from eight sampling sites in five different rivers, from two specimens provided by the Ichthyological Centre of Rome and mitochondrial sequences of Barbus spp. retrieved from GenBank. Presumptive classification based on external morphological characters was not confirmed by genetic analysis, by means of which all specimens collected in the Marches Region were ascribed to Barbus plebejus. Genetic diversity values (h and π) of sampling groups were all different from 0 except the one sample collected from the upper reaches of the River Tenna, above a hydroelectric dam. Population connectivity and colonization patterns of the studied area were inferred from an analysis of molecular variance distribution and evolutionary relationships among haplotypes. The results point to different levels of isolation among sampling groups due to ecological and anthropogenic factors and the effect of an artificial barrier on genetic variability and conservation status of the population. Finally, this study confirms the uncertainty associated with systematic classification of Barbus spp. based on morphological characters due to the phenotypic plasticity of the species.     

Involvement of serotonin transporter extracellular loop 1 in serotonin binding and transport

Residues Tyr-110 through Gly-115 of serotonin transporter were replaced, one at a time, with cysteine. Of these mutants, only G113C retained full activity for transport, Q111C and N112C retained partial activity, but Y110C, G114C and G115C were inactive. Poor surface expression was at least partly responsible for the lack of transport by G114C and G115C. In membrane preparations, Y110C through G113C all bound a high affinity cocaine analog similarly to the wild type. Treatment with methanethiosulfonate reagents increased the transport activity of Q111C and N112C to essentially wild-type levels but had no measurable effect on the other mutants. The decreased activity of Q111C and N112C resulted from an increase in the K_M for serotonin that was not accompanied by a decrease in serotonin binding affinity. Superfusion experiments indicated a defect in 5-HT exchange. Modification of the inserted cysteine residues reversed the increase in K_M and the poor exchange, also with no effect on serotonin affinity. The results suggest that Gln-111 and Asn-112 are not required for substrate binding but participate in subsequent steps in the transport cycle.