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.     

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.     

Use of Implantable Loop Recorders to Unravel the Cause of Unexplained Syncope

Syncope is a symptom of many underlying disease states, which range from the relatively benign to the life threatening. There are numerous investigations done for patients with recurrent unexplained syncope which may have very low yield when it comes to making a definitive diagnosis. Recently, the implantable loop recorder (ILR) for continuous monitoring of the cardiac rhythm has been launched in India. This review will briefly discuss these current availabel strategies and focus on the usefulness of an ILR in the definitive diagnosis and treatment of patients with a recurrent unexplained syncope.     

The Heterogeneity of the 7S Soybean Protein by Sepharose Gel Chromatography and Disc Gel Electrophoresis

Two kinds of N-acetylmuramidase, M-1 and M-2 enzymes, that were isolated from the cultural broth of Stm. globisporus 1829, were remarkably different in amino acid composition, immunological properties and modes of lytic action from each other. The M-1 enzyme was composed of 186 amino acid residues of which two moles were of half cystine, while the M-2 enzyme was composed of 99 amino acid residues with no cysteine. The hydrolyzing action of the M-2 enzyme was suppressed by the presence of an O-acetyl group on muramic acid residues in the peptidoglycan moiety, while that of the M-l enzyme was independent of the presence of O-acetyl groups. However, the hydrolyzing activity of both enzymes was enhanced when some muramic acid residues were substituted with stem peptides containing alanine, isoglutamine and lysine.     

Development of a loop‐mediated isothermal amplification assay for rapid and sensitive detection of Sporisorium scitamineum in sugarcane

Smut disease caused by Sporisorium scitamineum is one of the most destructive sugarcane diseases worldwide. The pathogen spreads primarily through infected sugarcane setts, and hence, the use of disease‐free planting materials is essential for preventing disease development in the field. In this study, a species‐specific loop‐mediated isothermal amplification (LAMP) assay was developed for rapid and accurate detection of S. scitamineum. Based on the differences in internal transcribed spacer (ITS) sequences of S. scitamineum, a set of four species‐specific primers, F3, B3, FIP and BIP, were designed by using a panel of fungal and bacterial species as controls. After optimization of the reaction conditions, the detection limit of LAMP assay was about 2 fg of the S. scitamineum genomic DNA in 25 µL reaction solution, 100‐fold lower than that of conventional polymerase chain reaction. The assay showed high specificity to discriminate all S. scitamineum isolates from nine other fungal and bacterial pathogens. The LAMP assay also detected smut infection from young sugarcane leaves with no visible smut‐disease symptoms. The findings from this study provide a simple, highly sensitive, rapid and reliable technique for early detection of S. scitamineum, which may be useful for sugarcane quarantine and production of smut‐free seedcanes. This is the first report of LAMP‐based assay for the detection of S. scitamineum in sugarcane.