A step toward the prediction of the fluorescence lifetimes of tryptophan residues in proteins based on structural and spectral data

A method is presented that allows the calculation of the lifetimes of tryptophan residues on the basis of spectral and structural data. It is applied to two different proteins. The calcium binding protein from the sarcoplasm of the muscles of the sand worm Nereis diversicolor (NSCP) changes its conformation upon binding of Ca2+ or Mg2+. NSCP contains three tryptophan residues at position 4, 57, and 170, respectively. The fluorescence lifetimes of W57 are investigated in a mutant in which W4 and W170 have been replaced. The time resolved fluorescence properties of W57 are linked to its different microconformations, which were determined by a molecular dynamics simulation map. Together with the determination of the radiative rate constant from the wavelength of maximum intensity of the decay associated spectra, it was possible to determine an exponential relation between the nonradiative rate constant and the distance between the indole CE3 atom and the carbonyl carbon of the peptide bond reflecting a mechanism of electron transfer as the main determinant of the value for the nonradiative rate constant. This result allows the calculation of the fluorescence lifetimes from the protein structure and the spectra. This method was further tested for the tryptophan of Ha-ras p21 (W32) and for W43 of the Tet repressor, which resulted in acceptable values for the predicted lifetimes.     

Dynamics and structural stability effects of germline PTEN mutations associated with cancer versus autism phenotypes

Individuals with germline mutations in the tumor suppressor gene phosphatase and tensin homolog (PTEN), irrespective of clinical presentation, are diagnosed with PTEN hamartoma tumor syndrome (PHTS). PHTS confers a high risk of breast, thyroid, and other cancers or autism spectrum disorder (ASD) with macrocephaly. It remains unclear why mutations in one gene can lead to seemingly disparate phenotypes. Thus, we sought to identify differences in ASD vs. cancer-associated germline PTEN missense mutations by investigating putative structural effects induced by each mutation. We utilized a theoretical computational approach combining in silico structural analysis and molecular dynamics (MD) to interrogate 17 selected mutations from our patient population: six mutations were observed in patients with ASD (only), six mutations in patients with PHTS-associated cancer (only), four mutations shared across both phenotypes, and one mutation with both ASD and cancer. We demonstrate structural stability changes where all six cancer-associated mutations showed a global decrease in structural stability and increased dynamics across the domain interface with a proclivity to unfold, mediating a closed (inactive) active site. In contrast, five of the six ASD-associated mutations showed localized destabilization that contribute to the partial opening of the active site. Our results lend insight into distinctive structural effects of germline PTEN mutations associated with PTEN-ASD vs. those associated with PTEN-cancer, potentially aiding in identification of the shared and separate molecular features that contribute to autism or cancer, thus, providing a deeper understanding of genotype–phenotype relationships for germline PTEN mutations.     

The Bio3D packages for structural bioinformatics

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D‐eddm package supports both experimental and theoretical simulation‐generated structures, is integrated with other methods for dissecting sequence‐structure–function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/ .     

Investigation of the molecular similarity in closely related protein systems: The PrP case study

The amyloid conversion is a massive detrimental modification affecting several proteins upon specific physical or chemical stimuli characterizing a plethora of diseases. In many cases, the amyloidogenic stimuli induce specific structural features to the protein conferring the propensity to misfold and form amyloid deposits. The investigation of mutants, structurally similar to their native isoform but inherently prone to amyloid conversion, may be a viable strategy to elucidate the structural features connected with amyloidogenesis. In this article, we present a computational protocol based on the combination of molecular dynamics (MD) and grid‐based approaches suited for the pairwise comparison of closely related protein structures. This method was applied on the cellular prion protein (PrP^C) as a case study and, in particular, addressed to the quali/quantification of the structural features conferred by either E200K mutations and treatment with CaCl₂, both able to induce the scrapie conversion of PrP. Several schemes of comparison were developed and applied to this case study, and made up suitable of application to other protein systems. At this purpose an in‐house python codes has been implemented that, together with the parallelization of the GRID force fields program, will spread the applicability of the proposed computational procedure. Proteins 2015; 83:1751–1765. © 2015 Wiley Periodicals, Inc.     

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations

Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in β-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.     

A communication network within the cytoplasmic domain of toll-like receptors has remained conserved during evolution

Toll/IL-1R (TIR) domain, that is, the cytoplasmic domain, in toll-like receptors (TLRs) from different species showed high sequence conservation in stretches spread across the surface as well as the core of the domain. To probe the structure–function significance of these residues, especially those coming from the core of TIR domains, we analyzed molecular dynamics trajectories of sequence similarity based models of human TIR domains. This study brought forth that N-terminal of the TIR domain simultaneously interacts with the flanking residues of the BB loop and central β-sheets. At the same time, residues of the central β-strands form favorable contacts with the DD loop and C-terminal, thus forming a two-way circuit between the N- and C-termini. In this work, the array of intradomain interactions is termed as communication network. Importantly, the “hubs” of this communication network were found to be conserved in all human TLRs. Earlier mutagenesis–function correlation work brought forth that certain mutations in the “core” of the TIR domain of TLR4 (e.g. in IFI767–769AAA and L815A) led to almost complete abrogation of signaling and reasoning for this dramatic loss-of-function has remained unclear, since these sites are not surface exposed. Using MD studies, we show here that this communication network gets disrupted in mutants of human TLR4 which were earlier reported to be functionally compromised. Extension of MD studies to heterodimer of TLR1/2 suggested that this evolutionarily conserved communication network senses the interactions formed upon dimerization and relays it to surfaces which are not involved in direct interdomain contacts.     

Solution Conformation of [D-Pen2,D-Pen5]enkephalin in Water: A NMR and Molecular Dynamics Study

The solution conformation of [D -Pen²,D -Pen⁵] enkephalin (DPDPE), a highly potent δ-selective opioid agonist, was examined by means of NMR, molecular mechanics and molecular dynamics methods. The structural information in the solvent water was obtained employing one- and two-dimensional methods of ¹H and ¹³C-NMR spectroscopy. Based on the distance geometry technique using the ROE data as input, 400 conformers were obtained and considered in the structure analysis. Alternatively, about 2000 conformers were stochastically generated and related to the NMR data after energy minimization. The structure analysis provides one conformer in agreement with all NMR data, which belongs to the lowest energy conformation group. This structure may serve as a reference conformer for DPDPE analogues synthesized with the aim of activity increase.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Mutational analysis of FUS gene and its structural and functional role in amyotrophic lateral sclerosis 6

Amyotrophic lateral sclerosis 6 (ALS6) is an autosomal recessive disorder caused by heterozygous mutation in the Fused in Sarcoma (FUS) gene. ALS6 is a neurodegenerative disorder, which affects the upper and lower motor neurons in the brain and spinal cord, resulting in fatal paralysis. ALS6 is caused by the genetic mutation in the proline/tyrosine-nuclear localization signals of the Fused in sarcoma Protein (FUS). FUS gene also known as TLS (Translocated in liposarcoma), which encodes a protein called RNA-binding protein-Fus (FUS), has a molecular weight of 75?kDa. In this analysis, we applied computational approach to filter the most deleterious and neurodegenerative disease of ALS6-associated mutation on FUS protein. We found H517Q as most deleterious and disease associated using PolyPhen 2.0, I-Mutant 3.0, SIFT, SNPs&GO, PhD-SNP, Pmut, and Mutpred tools. Molecular dynamics simulation (MDS) approach was conducted to investigate conformational changes in the mutant protein structure with respect to its native conformation. MDS results showed the flexibility loss in mutant (H517Q) FUS protein. Due to mutation, FUS protein became more rigid in nature and might alter the structural and functional behavior of protein and play a major role in inducing ALS6. The results obtained from this investigation would help in the field of pharmacogenomics to develop a potent drug target against FUS-associated neurodegenerative diseases.     

Theoretical studies on the interaction of guanine riboswitch with guanine and its closest analogues

Experimental studies (M. Mandal, B. Boese, J.E. Barrick, W.C. Winkler and R.R. Breaker, Riboswitches control fundamental biochemical pathways in bacillus subtilis and other bacteria, Cell 113 (2003), pp. 577–586) demonstrated that, besides recognising guanine with high specificity, guanine riboswitch could also bind guanine analogues, but the alteration of every functionalised position on the guanine heterocycle could cause a substantial loss of binding affinity. To investigate the nature of guanine riboswitch recognising metabolites, molecular docking and molecular dynamics simulation were carried out on diverse guanine analogues. The calculation results reveal that (1) most guanine analogues could bind to guanine riboswitch at the same binding pocket, with identical orientations and dissimilar binding energies, which is related to the positions of the functional groups; (2) the two tautomers of xanthine adopt different binding modes, and the enol-tautomer shows similar binding mode and affinity of hypoxanthine, which agrees well with the experimental results and (3) the riboswitch could form stable complexes with guanine analogues by hydrogen bonding contacts with U51 and C74. Particularly, U51 plays an important role in stabilising the complexes.     

Nonspecific yet decisive: Ubiquitination can affect the native‐state dynamics of the modified protein

Ubiquitination is one of the most common post‐translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site‐specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48‐linked tetra‐ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half‐life through proteasomal degradation.     

Congeneric but still distinct: how closely related trypsin ligands exhibit different thermodynamic and structural properties

A congeneric series of benzamidine-type ligands with a central proline moiety and a terminal cycloalkyl group—linked by a secondary amine, ether, or methylene bridge—was synthesized as trypsin inhibitors. This series of inhibitors was investigated by isothermal titration calorimetry, crystal structure analysis in two crystal forms, and molecular dynamics simulations. Even though all of these congeneric ligands exhibited essentially the same affinity for trypsin, their binding profiles at the structural, dynamic, and thermodynamic levels are very distinct. The ligands display a pronounced enthalpy/entropy compensation that results in a nearly unchanged free energy of binding, even though individual enthalpy and entropy terms change significantly across the series. Crystal structures revealed that the secondary amine-linked analogs scatter over two distinct conformational families of binding modes that occupy either the inside or of the outside the protein's S3/S4 specificity pocket. In contrast, the ether-linked and methylene-linked ligands preferentially occupy the hydrophobic specificity pocket. This also explains why the latter ligands could only be crystallized in the conformationally restricting closed crystal form whereas the derivative with the highest residual mobility in the series escaped our attempts to crystallize it in the closed form; instead, a well-resolved structure could only be achieved in the open form with the ligand in disordered orientation. These distinct binding modes are supported by molecular dynamics simulations and correlate with the shifting enthalpic/entropic signatures of ligand binding. The examples demonstrate that, at the molecular level, binding modes and thermodynamic binding signatures can be very different even for closely related ligands. However, deviating binding profiles provide the opportunity to optimally address a given target.     

A tale of two GTPases in cotranslational protein targeting

Guanosine triphosphatases (GTPases) comprise a superfamily of proteins that provide molecular switches to regulate numerous cellular processes. The "GTPase switch" paradigm, in which a GTPase acts as a bimodal switch that is turned "on" and "off" by external regulatory factors, has been used to interpret the regulatory mechanism of many GTPases. Recent work on a pair of GTPases in the signal recognition particle (SRP) pathway has revealed a distinct mode of GTPase regulation. Instead of the classical GTPase switch, the two GTPases in the SRP and SRP receptor undergo a series of conformational changes during their dimerization and reciprocal activation. Each conformational rearrangement provides a point at which these GTPases can communicate with and respond to their upstream and downstream biological cues, thus ensuring the spatial and temporal precision of all the molecular events in the SRP pathway. We suggest that the SRP and SRP receptor represent an emerging class of "multistate" regulatory GTPases uniquely suited to provide exquisite control over complex cellular pathways that require multiple molecular events to occur in a highly coordinated fashion.     

An Effective Solvation Term Based on Atomic Occupancies for Use in Protein Simulations

Abstract

Several approaches to the treatment of solvent effects based on continuum models are reviewed and a new method based on occupied atomic volumes (occupancies) is proposed and tested. The new method describes protein-water interactions in terms of atomic solvation parameters, which represent the solvation free energy per unit of volume. These parameters were determined for six different atoms types, using experimental free energies of solvation. The method was implemented in the GROMOS and PRESTO molecular simulation program suites. Simulations with the solvation term require 20-50% more CPU time than the corresponding vacuum simulations and are approximately 20 times faster than explicit water simulations. The method and parameters were tested by carrying out 200 ps simulations of BPTI in water, in vacuo, and with the solvation term. The performance of the solvation term was assessed by comparing the structures and energies from the solvation simulations with the equivalent quantities derived from several BPTI crystal structures and from the explicit water and vacuum simulations. The model structures were evaluated in terms of exposed total surface, buried and exposed polar surfaces, secondary structure preservation, number of hydrogen bonds, energy contributions, and positional deviations from BPTI crystal structures. Vacuum simulations produced unrealistic structures with respect to all criteria applied. The structures resulting from the simulations with explicit water were closer to the 5PTI crystal structure, although part of the secondary structure dissolved. The simulations with the effective solvation term produce structures that are normal according to all evaluations and in most respects are remarkably similar to the 5PTI crystal structure despite considerable positional fluctuations during the simulations. The segments where the model and crystal structures differ are known to be flexible and the observed difference may be physically realistic. The effective solvation term based on occupancies is not only very efficient in terms of computer time but also results in meaningful structural properties for BPTI. It may therefore be generally useful in molecular dynamics of macromolecules.     

CD4-gp120 interaction interface - a gateway for HIV-1 infection in human: molecular network,modeling and docking studies

The major causative agent for Acquired Immune Deficiency Syndrome (AIDS) is Human Immunodeficiency Virus-1 (HIV-1). HIV-1 is a predominant subtype of HIV which counts on human cellular mechanism virtually in every aspect of its life cycle. Binding of viral envelope glycoprotein-gp120 with human cell surface CD4 receptor triggers the early infection stage of HIV-1. This study focuses on the interaction interface between these two proteins that play a crucial role for viral infectivity. The CD4–gp120 interaction interface has been studied through a comprehensive protein–protein interaction network (PPIN) analysis and highlighted as a useful step towards identifying potential therapeutic drug targets against HIV-1 infection. We prioritized gp41, Nef and Tat proteins of HIV-1 as valuable drug targets at early stage of viral infection. Lack of crystal structure has made it difficult to understand the biological implication of these proteins during disease progression. Here, computational protein modeling techniques and molecular dynamics simulations were performed to generate three-dimensional models of these targets. Besides, molecular docking was initiated to determine the desirability of these target proteins for already available HIV-1 specific drugs which indicates the usefulness of these protein structures to identify an effective drug combination therapy against AIDS.     

Molecular dynamics simulations provide insights into the origin of gleevec’s selectivity toward human tyrosine kinases

Protein kinases are critical drug targets against cancer. Since the discovery of Gleevec, a specific inhibitor of Abl kinase, the capability of this drug to distinguish between Abl and other tyrosine kinases, such as Src, has been intensely investigated but the origin of Gleevec’s selectivity to Abl against Src is less studied. Here, we performed molecular dynamics (MD) simulations, dynamical cross-correlation matrices (DCCM), dynamical network analysis, and binding free energy calculations to explore Gleevec’s selectivity based on the crystal structures of Abl, Src, and their common ancestors (ANC-AS) and the two constructed mutation systems (AS→Abl and AS→Src). MD simulations revealed that the conformation of the phosphate-binding loop (P-loop) was altered significantly in the AS→Abl system. DCCM results unraveled that mutations increased anticorrelated motions in the AS→Abl system. Community network analysis suggested that the P-loop established special contacts in the AS→Abl system that are devoid in the AS→Src system. The binding free energy calculations unveiled that the affinity of Gleevec to AS→Abl increased to near the Abl level, whereas its affinity to AS→Src decreased to near the Src level. Analysis of individual residue contributions showed that the differences were located mainly at the P-loop. This study is valuable for understanding the sensitivity of Gleevec to human tyrosine kinases.

Communicated by Ramaswamy H. Sarma     

Ab initio simulations of Cu binding sites on the N-terminal region of prion protein

The human prion protein binds Cu²⁺ ions in the octarepeat domain of the N-terminal tail up to full occupancy at pH 7.4. Recent experiments have shown that the HGGG octarepeat subdomain is responsible for holding the metal bound in a square-planar configuration. By using first principle ab initio molecular dynamics simulations of the Car–Parrinello type, the coordination of copper to the binding sites of the prion protein octarepeat region is investigated. Simulations are carried out for a number of structured binding sites. Results for the complexes Cu(HGGGW)(wat), Cu(HGGG), and [Cu(HGGG)]₂ are presented. While the presence of a Trp residue and a water molecule does not seem to affect the nature of the copper coordination, high stability of the bond between copper and the amide nitrogen of deprotonated Gly residues is confirmed in all cases. For the more interesting [Cu(HGGG)]₂ complex, a dynamically entangled arrangement of the two domains with exchange of amide nitrogen bonds between the two copper centers emerges, which is consistent with the short Cu–Cu distance observed in experiments at full copper occupancy.     

Molecular dynamics simulation of the substitution of serine for the conserved glycine in the G-loop in the cdc28-srm yeast mutant using the crystal lattice of human CDK2 kinase

Two-nanosecond molecular dynamics simulations of the crystal lattice of an active complex of pT160-CDK2 kinase/cyclin A/ATP-Mg²⁺/substrate were performed. The simulations showed that the structures of the wild-type CDK2 complex and the mutant CDK2 complex involving the substitution G16S-CDK2 corresponding to the yeast substitution G20S-CDC28 differ noticeably and the differences between the structural conformations are most pronounced in the regions that play a key role in the kinase functioning. The results of the computer calculations were used to consider the structural elements that may affect the kinase activity, the regulatory phosphorylation, and the binding of protein kinase with cyclins and substrates.     

Motional timescale predictions by molecular dynamics simulations: Case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin‐lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.