Thermodynamically Important Contacts in Folding of Model Proteins

We introduce a quantity, the entropic susceptibility, that measures the thermodynamic importance-for the folding transition-of the contacts between amino acids in model proteins. Using this quantity, we find that only one equilibrium run of a computer simulation of a model protein is sufficient to select a subset of contacts that give rise to the peak in the specific heat observed at the folding transition. To illustrate the method, we identify thermodynamically important contacts in a model 46-mer. We show that only about 50% of all contacts present in the protein native state are responsible for the sharp peak in the specific heat at the folding transition temperature, while the remaining 50% of contacts do not affect the specific heat.Comment: 5 pages, 5 figures; to be published in PR

The resource potential of the agro-industrial complex of the north-caucasian federal district as a factor of innovative development of the macro-region

Purpose: Innovative development of the region is largely determined by the available resources and the ability of the economic system to reproduce them. The article deals with the problem of increasing the efficiency of the use of the resource potential of the North-Caucasian Federal District (NCFD). Design/Methodology/Approach: The elements of the structure of the resource potential of the agro-region, including natural, material, technical, labor and innovation component are revealed. An important factor in the sustainable socio-economic development of the North-Caucasian Federal District (NCFD) is the production of high-quality products of the agro-industrial complex based on maximum use of local raw materials and own resources. Findings: To modernize the agricultural regional economy, it is necessary to base agro-industrial production mainly on innovative activities that allow the creation of a single production space. Practical implications: The results of the study can be used in the target complex programs of development of regional economy with the aim of strengthening and enhancing the effectiveness of regional institutional and industrial transformations. Originality/Value: The resource potential is proposed to be determined in the context of the availability and limitations of regional production conditions.peer-reviewe

Ab initio RNA folding by discrete molecular dynamics: From structure prediction to folding mechanisms

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 Å deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNAPhe, pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses

Molecular Dynamics Simulation of Amyloid Beta Dimer Formation

Recent experiments with amyloid-beta (Abeta) peptide suggest that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation: (1) we use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations, and (2) we employ all-atom molecular mechanics simulations to estimate the thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts ten different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that (a) all dimer conformations have higher free energies compared to their corresponding monomeric states, and (b) the free energy difference between the Abeta(1-42) and the analogous Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of stable planar beta-strand Abeta dimers.Comment: 32 pages (preprint format), 3 figure

Molecular Origin of Polyglutamine Aggregation in Neurodegenerative Diseases

Expansion of polyglutamine (polyQ) tracts in proteins results in protein aggregation and is associated with cell death in at least nine neurodegenerative diseases. Disease age of onset is correlated with the polyQ insert length above a critical value of 35–40 glutamines. The aggregation kinetics of isolated polyQ peptides in vitro also shows a similar critical-length dependence. While recent experimental work has provided considerable insights into polyQ aggregation, the molecular mechanism of aggregation is not well understood. Here, using computer simulations of isolated polyQ peptides, we show that a mechanism of aggregation is the conformational transition in a single polyQ peptide chain from random coil to a parallel β-helix. This transition occurs selectively in peptides longer than 37 glutamines. In the β-helices observed in simulations, all residues adopt β-strand backbone dihedral angles, and the polypeptide chain coils around a central helical axis with 18.5 ± 2 residues per turn. We also find that mutant polyQ peptides with proline-glycine inserts show formation of antiparallel β-hairpins in their ground state, in agreement with experiments. The lower stability of mutant β-helices explains their lower aggregation rates compared to wild type. Our results provide a molecular mechanism for polyQ-mediated aggregation

Hybrid Dynamics Simulation Engine for Metalloproteins

Quality computational description of metalloproteins is a great challenge due to the vast span of time- and lengthscales characteristic of their existence. We present an efficient new method that allows for robust characterization of metalloproteins. It combines quantum mechanical (QM) description of the metal-containing active site, and extensive dynamics of the protein captured by discrete molecular dynamics (DMD) (QM/DMD). DMD samples the entire protein, including the backbone, and most of the active site, except for the immediate coordination region of the metal. QM operates on the part of the protein of electronic and chemical significance, which may include tens to hundreds of atoms. The breathing quantum-classical boundary provides a continuous mutual feedback between the two machineries. We test QM/DMD using the Fe-containing electron transporter protein, rubredoxin, and its three mutants as a model. QM/DMD can provide a reliable balanced description of metalloproteins’ structure, dynamics, and electronic structure in a reasonable amount of time. As an illustration of QM/DMD capabilities, we then predict the structure of the Ca2+ form of the enzyme catechol O-methyl transferase, which, unlike the native Mg2+ form, is catalytically inactive. The Mg2+ site is ochtahedral, but the Ca2+ is 7-coordinate and features the misalignment of the reacting parts of the system. The change is facilitated by the backbone adjustment. QM/DMD is ideal and fast for providing this level of structural insight

Dynamics Intrinsic to Cystic Fibrosis Transmembrane Conductance Regulator Function and Stability

The cystic fibrosis transmembrane conductance regulator (CFTR) requires dynamic fluctuations between states in its gating cycle for proper channel function, including changes in the interactions between the nucleotide-binding domains (NBDs) and between the intracellular domain (ICD) coupling helices and NBDs. Such motions are also linked with fluctuating phosphorylation-dependent binding of CFTR’s disordered regulatory (R) region to the NBDs and partners. Folding of CFTR is highly inefficient, with the marginally stable NBD1 sampling excited states or folding intermediates that are aggregation-prone. The severe CF–causing F508del mutation exacerbates the folding inefficiency of CFTR and leads to impaired channel regulation and function, partly as a result of perturbed NBD1–ICD interactions and enhanced sampling of these NBD1 excited states. Increased knowledge of the dynamics within CFTR will expand our understanding of the regulated channel gating of the protein as well as of the F508del defects in folding and function

Discrete molecular dynamics simulations of peptide aggregation

We study the aggregation of peptides using the discrete molecular dynamics simulations. At temperatures above the alpha-helix melting temperature of a single peptide, the model peptides aggregate into a multi-layer parallel beta-sheet structure. This structure has an inter-strand distance of 0.48 nm and an inter-sheet distance of 1.0 nm, which agree with experimental observations. In this model, the hydrogen bond interactions give rise to the inter-strand spacing in beta-sheets, while the Go interactions among side chains make beta-strands parallel to each other and allow beta-sheets to pack into layers. The aggregates also contain free edges which may allow for further aggregation of model peptides to form elongated fibrils.Comment: 15 pages, 8 figure