Mechanisms of the acute ischemia-induced arrhythmogenesis--a simulation study

The underlying ionic mechanisms of ischemic-induced arrhythmia were studied by the computer simulation method. To approximate the real situation, ischemic cells were simulated by considering the three major component conditions of acute ischemia (elevated extracellular K(+) concentration, acidosis and anoxia) at the level of ionic currents and ionic concentrations, and a round ischemic zone was introduced into a homogeneous healthy sheet to avoid sharp angle of the ischemic tissue. The constructed models were solved using the operator splitting and adaptive time step methods, and the perturbation finite difference (PFD) scheme was first used to integrate the partial differential equations (PDEs) in the model. The numerical experiments showed that the action potential durations (APDs) of ischemic cells did not exhibited rate adaptation characteristic, resulting in flattening of the APD restitution curve. With reduction of sodium channel availability and long recovery of excitability, refractory period of the ischemic tissue was significantly prolonged, and could no longer be considered as same as APD. Slope of the conduction velocity (CV) restitution curve increased both in normal and ischemic region when pacing cycle length (PCL) was short, and refractory period dispersion increased with shortening of PCL as well. Therefore, dynamic changes of CV and dispersion of refractory period rather than APD were suggested to be the fundamental mechanisms of arrhythmia in regional ischemic myocardium.     

The putative propeptide of MycP1 in mycobacterial type VII secretion system does not inhibit protease activity but improves protein stability

Mycosin-1 protease (MycP1) is a serine protease anchored to the inner membrane of Mycobacterium tuberculosis, and is essential in virulence factor secretion through the ESX-1 type VII secretion system (T7SS). Bacterial physiology studies demonstrated that MycP1 plays a dual role in the regulation of ESX-1 secretion and virulence, primarily through cleavage of its secretion substrate EspB. MycP1 contains a putative N-terminal inhibitory propeptide and a catalytic triad of Asp-His-Ser, classic hallmarks of a subtilase family serine protease. The MycP1 propeptide was previously reported to be initially inactive and activated after prolonged incubation. In this study, we have determined crystal structures of MycP1 with (MycP1^24-422) and without (MycP1^63-422) the propeptide, and conducted EspB cleavage assays using the two proteins. Very high structural similarity was observed in the two crystal structures. Interestingly, protease assays demonstrated positive EspB cleavage for both proteins, indicating that the putative propeptide does not inhibit protease activity. Molecular dynamic simulations showed higher rigidity in regions guarding the entrance to the catalytic site in MycP1^24-422 than in MycP1^63-422, suggesting that the putative propeptide might contribute to the conformational stability of the active site cleft and surrounding regions.     

Identification of the active site of human mitochondrial malonyl‐coenzyme a decarboxylase: A combined computational study

Malonyl‐CoA decarboxylase (MCD) can control the level of malonyl‐CoA in cell through the decarboxylation of malonyl‐CoA to acetyl‐CoA, and plays an essential role in regulating fatty acid metabolism, thus it is a potential target for drug discovery. However, the interactions of MCD with CoA derivatives are not well understood owing to unavailable crystal structure with a complete occupancy in the active site. To identify the active site of MCD, molecular docking and molecular dynamics simulations were performed to explore the interactions of human mitochondrial MCD (HmMCD) and CoA derivatives. The findings reveal that the active site of HmMCD indeed resides in the prominent groove which resembles that of CurA. However, the binding modes are slightly different from the one observed in CurA due to the occupancy of the side chain of Lys183 from the N‐terminal helical domain instead of the adenine ring of CoA. The residues 300 ? 305 play an essential role in maintaining the stability of complex mainly through hydrogen bond interactions with the pyrophosphate moiety of acetyl‐CoA. Principle component analysis elucidates the conformational distribution and dominant concerted motions of HmMCD. MM_PBSA calculations present the crucial residues and the major driving force responsible for the binding of acetyl‐CoA. These results provide useful information for understanding the interactions of HmMCD with CoA derivatives. Proteins 2016; 84:792–802. © 2016 Wiley Periodicals, Inc.     

Atomistic mechanisms of huntingtin N‐terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q_N) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N‐terminal 17‐amino‐acid sequence (htt^NT) positioned just before the Q_N region is important for the binding of huntingtin to membranes. Through all‐atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt^NTQ_N fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps—approach, reorganization, anchoring, and insertion—that are very diverse at the atomic level. On the membrane, the htt^NT peptide forms a stable α‐helix essentially parallel to the membrane with its nonpolar side‐chains—mainly Leu‐4, Leu‐7, Phe‐11 and Leu‐14—positioned in the hydrophobic core of the membrane. Salt‐bridges involving Glu‐5, Glu‐12, Lys‐6, and Lys‐15, as well as hydrogen bonds involving Thr‐3 and Ser‐13 with the phospholipids also stabilize the structure and orientation of the htt^NT peptide. These observations do not significantly change upon adding the Q_N region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt^NT region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q_N region could also play a key role for the oligomerization of htt^NTQ_N on phospholipid membranes. Proteins 2014; 82:1409–1427. © 2014 Wiley Periodicals, Inc.     

Viability of the Endangered Little Bustard <Emphasis Type="Italic">Tetrax tetrax</Emphasis> Population of Western France

Stochastic computer simulations are used to evaluate the sensitivity of Little bustard population parameters, estimating the survival probabilities of the seven endangered Little bustard populations of central-western France for which conservation actions are currentlybeing or have been implemented. Different scenarios of parameter compensation for those nuclei to establish parameter levels assuring population viability are discussed. Adult survival, productivity per female, initial population size and carrying capacity were the most sensitive parameters in a hypothetical, isolated population. Juvenile survival also affected population survival, although its sensitivity was lower. Sex ratio did not have a linear effect on population survival, but probability of extinction increased for extreme values. Productivity per female and initial population size, varied strongly among the populations studied, determining their average time of extinction and growth rate. When a metapopulation scenario was simulated, the survival probabilities of each population and the metapopulations stayed close to 1.0 if no mortality was associated to migration. When mortality during migration was included in the simulations, the metapopulation's probability of survival significantly decreased under 90%. This approach may help managers to correctly address conservation measures and design effective strategies, which should be directed mainly to improve productivity, enhance female survival, and minimise mortality during migration (e.g. promoting insect-rich nesting substrates, avoiding female killing and nest destruction at harvesting, reducing the risk of collision with powerlines, or controlling poaching).     

All-atom Monte Carlo simulation of GCAA RNA folding

We report a detailed all-atom simulation of the folding of the GCAA RNA tetraloop. The GCAA tetraloop motif is a very common and thermodynamically stable secondary structure in natural RNAs. We use our simulation methods to study the folding behavior of a 12-base GCAA tetraloop structure with a four-base helix adjacent to the tetraloop proper. We implement an all-atom Monte Carlo (MC) simulation of RNA structural dynamics using a Go potential. Molecular dynamics (MD) simulation of RNA and protein has realistic energetics and sterics, but is extremely expensive in terms of computational time. By coarsely treating non-covalent energetics, but retaining all-atom sterics and entropic effects, all-atom MC techniques are a useful method for the study of protein and now RNA. We observe a sharp folding transition for this structure, and in simulations at room temperature the state histogram shows three distinct minima: an unfolded state (U), a more narrow intermediated state (I), and a narrow folded state (F). The intermediate consists primarily of structures with the GCAA loop and some helix hydrogen bonds formed. Repeated kinetic folding simulations reveal that the number of helix base-pairs forms a simple 1D reaction coordinate for the I-->N transition.     

A Robust Water Potential Parameterisation

Abstract

We compare molecular dynamics simulation results for the properties of liquid water predicted by four novel water potential models. These models are designed as a combination of parameters taken from the dedicated but brittle TIP3P water potential, and the more flexible but less accurate parameterisations such as the Dreiding and Universal force fields. We find that a hybrid of Dreiding and TIP3P delivers the best results, yielding a density, diffusion coefficient and radial distribution function in good agreement with experiment, performing in some respects even better than the dedicated reference TIP3P model. Another Dreiding based force field predicts semi-quantitative results for the water structure and dynamics while the Universal force field based models are incapable of simulating a condensed phase of water at all, continuing to expand indefinitely. These observations are useful for selecting and designing robust water force field parameterisations that can be used for general simulation purposes.     

Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin’s exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate–aromatic interactions including CH–π and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to β-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.     

Modeling of Branching Patterns in Plants

A major determinant of plant architecture is the arrangement of branches around the stem, known as phyllotaxis. However, the specific form of branching conditions is not known. Here we discuss this question and suggest a branching model which seems to be in agreement with biological observations. Recently, a number of models connected with the genetic network or molecular biology regulation of the processes of pattern formation appeared. Most of these models consider the plant hormone, auxin, transport and distribution in the apical meristem as the main factors for pattern formation and phyllotaxis. However, all these models do not take into consideration the whole plant morphogenesis, concentrating on the events in the shoot or root apex. On the other hand, other approaches for modeling phyllotaxis, where the whole plant is considered, usually are mostly phenomenological, and due to it, do not describe the details of plant growth and branching mechanism. In this work, we develop a mathematical model and study pattern formation of the whole, though simplified, plant organism where the main physiological factors of plant growth and development are taken into consideration. We model a growing plant as a system of intervals, which we will consider as branches. We assume that the number and location of the branches are not given a priori, but appear and grow according to certain rules, elucidated by the application of mathematical modeling. Four variables are included in our model: concentrations of the plant hormones auxin and cytokinin, proliferation and growth factor, and nutrients—we observe a wide variety of plant forms and study more specifically the involvement of each variable in the branching process. Analysis of the numerical simulations shows that the process of pattern formation in plants depends on the interaction of all these variables. While concentrations of auxin and cytokinin determine the appearance of a new bud, its growth is determined by the concentrations of nutrients and proliferation factors. Possible mechanisms of apical domination in the frame of our model are discussed.     

Dependence of folding dynamics and structural stability on the location of a hydrophobic pair in beta-hairpins

We study the dependence of folding time, nucleation site, and stability of a model beta-hairpin on the location of a cross-strand hydrophobic pair, using a coarse-grained off-lattice model with the aid of Monte Carlo simulations. Our simulations have produced 6500 independent folding trajectories dynamically, forming the basis for extensive statistical analysis. Four folding pathways, zipping-out, middle-out, zipping-in, and reptation, have been closely monitored and discussed in all seven sequences studied. A hydrophobic pair placed near the beta-turn or in the middle section effectively speed up folding; a hydrophobic pair placed close to the terminal ends or next to the beta-turn encourages stability of the entire chain.