PRIMO/PRIMONA: A coarse‐grained model for proteins and nucleic acids that preserves near‐atomistic accuracy

The new coarse graining model PRIMO/PRIMONA for proteins and nucleic acids is proposed. This model combines one to several heavy atoms into coarse‐grained sites that are chosen to allow an analytical, high‐resolution reconstruction of all‐atom models based on molecular bonding geometry constraints. The accuracy of proposed reconstruction method in terms of structure and energetics is tested and compared with other popular reconstruction methods for a variety of protein and nucleic acid test sets. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

A coarse‐grained potential for fold recognition and molecular dynamics simulations of proteins

A coarse‐grained potential for protein simulations and fold ranking is presented. The potential is based on a two‐point model of individual amino acids and a specific implementation of hydrogen bonding. Parameters are determined for distance dependent pair interactions, pseudo bonds, angles, and torsions. A scaling factor for a hydrogen bonding term is also determined. Iterative sampling for 4867 proteins reproduces distributions of internal coordinates and distances observed in the Protein Data Bank. The adjustment of the potential and resampling are in the spirit of the generalized ensemble approach. No native structure information (e.g., secondary structure) is used in the calculation of the potential or in the simulation of a particular protein. The potential is subject to two tests as follows: (i) simulations of 956 globular proteins in the neighborhood of their native folds (these proteins were not used in the training set) and (ii) discrimination between native and decoy structures for 2470 proteins with 305,000 decoys and the “Decoys ‘R’ Us” dataset. In the first test, 58% of tested proteins stay within 5 Å from the native fold in Molecular Dynamics simulations of more than 20 nanoseconds using the new potential. The potential is also useful in differentiating between correct and approximate folds providing significant signal for structure prediction algorithms. Sampling with the potential consistently regenerates the distribution of distances and internal coordinates it learned. Nevertheless, during Molecular Dynamics simulations structures are found that reproduce the learned distributions but are far from the native fold. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Folding and stability of helical bundle proteins from coarse‐grained models

We develop a coarse‐grained model where solvent is considered implicitly, electrostatics are included as short‐range interactions, and side‐chains are coarse‐grained to a single bead. The model depends on three main parameters: hydrophobic, electrostatic, and side‐chain hydrogen bond strength. The parameters are determined by considering three level of approximations and characterizing the folding for three selected proteins (training set). Nine additional proteins (containing up to 126 residues) as well as mutated versions (test set) are folded with the given parameters. In all folding simulations, the initial state is a random coil configuration. Besides the native state, some proteins fold into an additional state differing in the topology (structure of the helical bundle). We discuss the stability of the native states, and compare the dynamics of our model to all atom molecular dynamics simulations as well as some general properties on the interactions governing folding dynamics. Proteins 2013; 81:1200–1211. © 2013 Wiley Periodicals, Inc.     

Improved coarse‐grained model for studying sequence dependent phase separation of disordered proteins

We present improvements to the hydropathy scale (HPS) coarse‐grained (CG) model for simulating sequence‐specific behavior of intrinsically disordered proteins (IDPs), including their liquid–liquid phase separation (LLPS). The previous model based on an atomistic hydropathy scale by Kapcha and Rossky (KR scale) is not able to capture some well‐known LLPS trends such as reduced phase separation propensity upon mutations (R‐to‐K and Y‐to‐F). Here, we propose to use the Urry hydropathy scale instead, which was derived from the inverse temperature transitions in a model polypeptide with guest residues X. We introduce two free parameters to shift (Δ) and scale (µ) the overall interaction strengths for the new model (HPS‐Urry) and use the experimental radius of gyration for a diverse group of IDPs to find their optimal values. Interestingly, many possible (Δ, µ) combinations can be used for typical IDPs, but the phase behavior of a low‐complexity (LC) sequence FUS is only well described by one of these models, which highlights the need for a careful validation strategy based on multiple proteins. The CG HPS‐Urry model should enable accurate simulations of protein LLPS and provide a microscopically detailed view of molecular interactions.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

Coarse‐ and fine‐grained models for proteins: Evaluation by decoy discrimination

Coarse‐grained models for protein structure are increasingly used in simulations and structural bioinformatics. In this study, we evaluated the effectiveness of three granularities of protein representation based on their ability to discriminate between correctly folded native structures and incorrectly folded decoy structures. The three levels of representation used one bead per amino acid (coarse), two beads per amino acid (medium), and all atoms (fine). Multiple structure features were compared at each representation level including two‐body interactions, three‐body interactions, solvent exposure, contact numbers, and angle bending. In most cases, the all‐atom level was most successful at discriminating decoys, but the two‐bead level provided a good compromise between the number of model parameters which must be estimated and the accuracy achieved. The most effective feature type appeared to be two‐body interactions. Considering three‐body interactions increased accuracy only marginally when all atoms were used and not at all in medium and coarse representations. Though two‐body interactions were most effective for the coarse representations, the accuracy loss for using only solvent exposure or contact number was proportionally less at these levels than in the all‐atom representation. We propose an optimization method capable of selecting bead types of different granularities to create a mixed representation of the protein. We illustrate its behavior on decoy discrimination and discuss implications for data‐driven protein model selection. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Analyses of the folding properties of ferredoxin‐like fold proteins by means of a coarse‐grained Gō model: Relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi‐state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on‐going research issues. In this study, we investigate ferredoxin‐like fold proteins which have a simple topology and multi‐state transitions. We analyze the folding processes by means of a coarse‐grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high‐free‐energy intermediate state, and ADA2h and S6, which fold into the native structure through two‐state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54‐55, are reproduced and compared with experimental observations. We show that the ferredoxin‐like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two‐state mechanism. Proteins 2014; 82:954–965. © 2013 Wiley Periodicals, Inc.     

Toward decrypting the allosteric mechanism of the ryanodine receptor based on coarse‐grained structural and dynamic modeling

The ryanodine receptors (RyRs) are a family of calcium (Ca) channels that regulate Ca release by undergoing a closed‐to‐open gating transition in response to action potential or Ca binding. The allosteric mechanism of RyRs gating, which is activated/regulated by ligand/protein binding >200 Å away from the channel gate, remains elusive for the lack of high‐resolution structures. Recent solution of the closed‐form structures of the RyR1 isoform by cryo‐electron microscopy has paved the way for detailed structure‐driven studies of RyRs functions. Toward elucidating the allosteric mechanism of RyRs gating, we performed coarse‐grained modeling based on the newly solved closed‐form structures of RyR1. Our normal mode analysis captured a key mode of collective motions dominating the observed structural variations in RyR1, which features large outward and downward movements of the peripheral domains with the channel remaining closed, and involves hotspot residues that overlap well with key functional sites and disease mutations. In particular, we found a key interaction between a peripheral domain and the Ca‐binding EF hand domain, which may allow for direct coupling of Ca binding to the collective motions as captured by the above mode. This key mode was robustly reproduced by the normal mode analysis of the other two closed‐form structures of RyR1 solved independently. To elucidate the closed‐to‐open conformational changes in RyR1 with amino‐acid level of details, we flexibly fitted the closed‐form structures of RyR1 into a 10‐Å cryo‐electron microscopy map of the open state. We observed extensive structural changes involving the peripheral domains and the central domains, resulting in the channel pore opening. In sum, our findings have offered unprecedented structural and dynamic insights to the allosteric mechanism of RyR1 via modulation of the key collective motions involved in RyR1 gating. The predicted hotspot residues and open‐form conformation of RyR1 will guide future mutational and functional studies. Proteins 2015; 83:2307–2318. © 2015 Wiley Periodicals, Inc.     

Application of information theory to a three‐body coarse‐grained representation of proteins in the PDB: Insights into the structural and evolutionary roles of residues in protein structure

Knowledge‐based methods for analyzing protein structures, such as statistical potentials, primarily consider the distances between pairs of bodies (atoms or groups of atoms). Considerations of several bodies simultaneously are generally used to characterize bonded structural elements or those in close contact with each other, but historically do not consider atoms that are not in direct contact with each other. In this report, we introduce an information‐theoretic method for detecting and quantifying distance‐dependent through‐space multibody relationships between the sidechains of three residues. The technique introduced is capable of producing convergent and consistent results when applied to a sufficiently large database of randomly chosen, experimentally solved protein structures. The results of our study can be shown to reproduce established physico‐chemical properties of residues as well as more recently discovered properties and interactions. These results offer insight into the numerous roles that residues play in protein structure, as well as relationships between residue function, protein structure, and evolution. The techniques and insights presented in this work should be useful in the future development of novel knowledge‐based tools for the evaluation of protein structure. Proteins 2014; 82:3450–3465. © 2014 Wiley Periodicals, Inc.     

A new reaction useful for chemical cross-linking between nucleic acids and proteins

Cytosine in nucleic acids can be converted into N4-aminocytosine by treatment with a mixture of hydrazine and bisulfite. The hydrazino group thus formed at position 4 of the pyrimidine ring can be linked to a sulhydryl group in proteins by the use of bromopyruvate as a linker. Successful use of this scheme of chemical cross-linking between nucleic acid and protein was demonstrated in the linking of poly(C) with glutathione, and of RNA with protein in the E. coli 30 S ribosomal subunit.     

Exploration of conformational changes in lactose permease upon sugar binding and proton transfer through coarse‐grained simulations

Escherichia coli lactose permease (LacY) actively transports lactose and other galactosides across cell membranes through lactose/H⁺ symport process. Lactose/H⁺ symport is a highly complex process that involves sugar translocation, H⁺ transfer, and large‐scale protein conformational changes. The complete picture of lactose/H⁺ symport is largely unclear due to the complexity and multiscale nature of the process. In this work, we develop the force field for sugar molecules compatible with PACE, a hybrid and coarse‐grained force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid. After validation, we implement the new force field to investigate the binding of a ‐d ‐galactopyranosyl‐1‐thio‐ ‐d ‐galactopyranoside (TDG) molecule to a wild‐type LacY. Results show that the local interactions between TDG and LacY at the binding pocket are consistent with the X‐ray experiment. Transitions from inward‐facing to outward‐facing conformations upon TDG binding and protonation of Glu269 have been achieved from ～5.5 µs simulations. Both the opening of the periplasmic side and closure of the cytoplasmic side of LacY are consistent with double electron–electron resonance and thiol cross‐linking experiments. Our analysis suggests that the conformational changes of LacY are a cumulative consequence of interdomain H‐bonds breaking at the periplasmic side, interdomain salt‐bridge formation at the cytoplasmic side, and the TDG orientational changes during the transition.     

A fine‐grained spatial prediction model for the red‐listed vascular plant Scorzonera humilis

We explored the applied use of distribution modelling as a tool for making spatial predictions of occurrences of the red‐listed vascular plant species Scorzonera humilis in a study area in southeast Norway. Scorzonera is typical of extensively managed semi‐natural grasslands. A Maxent model was trained on all known records of the species, accurately georeferenced and gridded to fine resolution (grid cells of 25×25 m). Model performance was assessed on the training data by data‐splitting (by which some records were set off for evaluation) and on independent evaluation data collected in the field. Of the eight predictor variables used in the modelling, distance to roads and to arable land were most important followed by land‐cover class and altitude. Judged from the area under curve (AUC), the model was good to excellent and a significant, positive relationship was found between relative probabilities of occurrence predicted by the model and true probability of presence provided by the independently collected evaluation data. The model was used together with the evaluation data to estimate presence of Scorzonera humilis in 0.7% of the grid cells in the study area. The grid cells in which the model predicted highest probability for Scorzonera to be present had a true probability of presence of ca 12%, i.e. 17×higher than in an average cell. The present study demonstrates that, even when only simple predictor variables are available, spatial prediction modelling contributes important knowledge about rare species such as prevalence estimates, spatial prediction maps and insights into the species’ autecology. Spatial prediction modelling also makes cost‐efficient monitoring of rare species possible. However, it is pointed out that these benefits require evaluation of the model on independently sampled evaluation data.     

Nucleobase‐caged peptide nucleic acids: PNA/PNA duplex destabilization and light‐triggered PNA/PNA recognition

The 2‐(o‐nitrophenyl)‐propyl (NPP) group is used as caging group to mask the nucleobases adenine and cytosine in N‐(2‐aminoethyl)glycine peptide nucleic acids (aeg‐PNA). The adeninyl and cytosinyl nucleo amino acid building blocks Fmoc‐a^NPP‐aeg‐OH and Fmoc‐c^NPP‐aeg‐OH were synthesized and incorporated into PNA sequences by Fmoc solid phase synthesis relying on high stability of the NPP nucleobase protecting group toward Fmoc‐cleavage, coupling, capping, and resin cleavage conditions. Removal of the nucleobase caging group was achieved by UV‐LED irradiation at 365 nm. The nucleobase caging groups provided sterical crowding effecting the Watson–Crick base pairing, and thereby, the PNA double strand stabilities. Duplex formation can completely be suppressed for complementary PNA containing caging groups in both strands. PNA/PNA recognition can be completely restored by UV light‐triggered release of the photolabile protecting group. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.