Project description:Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK a values for each nucleotide.
Project description:The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.
Project description:The interaction of ?-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an ?-helix within a membrane, or that change the probability that the ?-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an ?-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.
Project description:Holliday junctions play a central role in genetic recombination, DNA repair and other cellular processes. We combine simulations and experiments to evaluate the ability of the 3SPN.2 model, a coarse-grained representation designed to mimic B-DNA, to predict the properties of DNA Holliday junctions. The model reproduces many experimentally determined aspects of junction structure and stability, including the temperature dependence of melting on salt concentration, the bias between open and stacked conformations, the relative populations of conformers at high salt concentration, and the inter-duplex angle (IDA) between arms. We also obtain a close correspondence between the junction structure evaluated by all-atom and coarse-grained simulations. We predict that, for salt concentrations at physiological and higher levels, the populations of the stacked conformers are independent of salt concentration, and directly observe proposed tetrahedral intermediate sub-states implicated in conformational transitions. Our findings demonstrate that the 3SPN.2 model captures junction properties that are inaccessible to all-atom studies, opening the possibility to simulate complex aspects of junction behavior.
Project description:The coronavirus disease 19 (COVID-19) pandemic is causing a global health crisis and has already caused a devastating societal and economic burden. The pathogen, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has a high sequence and architecture identity with SARS-CoV, but far more people have been infected by SARS-CoV-2. Here, combining the structural data from cryo-electron microscopy and structure prediction, we constructed bottom-up Martini coarse-grained models of intact SARS-CoV and SARS-CoV-2 envelopes. Microsecond molecular dynamics simulations were performed, allowing us to explore their dynamics and supramolecular organization. Both SARS-CoV and SARS-CoV-2 envelopes present a spherical morphology, with structural proteins forming multiple string-like islands in the membrane and clusters between the heads of spike proteins. Critical differences between the SARS-CoV and SARS-CoV-2 envelopes are the interaction pattern between the spike proteins and the flexibility of the spike proteins. Our models provide structural and dynamic insights into the SARS virus envelopes and could be used for further investigation, such as drug design and membrane fusion and fission processes.
Project description:Many biological cellular processes occur at the micro- or millisecond time scale. With traditional all-atom molecular modeling techniques it is difficult to investigate the dynamics of long time scales or large systems, such as protein aggregation or activation. Coarse graining (CG) can be used to reduce the number of degrees of freedom in such a system, and reduce the computational complexity. In this paper the first version of a coarse grained model for transmembrane proteins is presented. This model differs from other coarse grained protein models due to the introduction of a novel angle potential as well as a hydrogen bonding potential. These new potentials are used to stabilize the backbone. The model has been validated by investigating the adaptation of the hydrophobic mismatch induced by the insertion of WALP-peptides into a lipid membrane, showing that the first step in the adaptation is an increase in the membrane thickness, followed by a tilting of the peptide.
Project description:Deoxyribozymes (DNAzymes) are single-stranded DNA that catalyze nucleic acid biochemistry. Although a number of DNAzymes have been discovered by in vitro selection, the relationship between their tertiary structure and function remains unknown. We focus here on the well-studied 10-23 DNAzyme, which cleaves mRNA with a catalytic efficiency approaching that of RNase A. Using coarse-grained Brownian dynamics simulations, we find that the DNAzyme bends its substrate away from the cleavage point, exposing the reactive site and buckling the DNAzyme catalytic core. This hypothesized transition state provides microscopic insights into experimental observations concerning the size of the DNAzyme/substrate complex, the impact of the recognition arm length, and the sensitivity of the enzymatic activity to point mutations of the catalytic core. Upon cleaving the pertinent backbone bond in the substrate, we find that the catalytic core of the DNAzyme unwinds and the overall complex rapidly extends, in agreement with experiments on the related 8-17 DNAzyme. The results presented here provide a starting point for interpreting experimental data on DNAzyme kinetics, as well as developing more detailed simulation models. The results also demonstrate the limitations of using a simple physical model to understand the role of point mutations.
Project description:An algorithm is described for the interaction of a hierarchy of objects that seeks to circumvent a fundamental problem in coarse-grained modelling which is the loss of fine detail when components become bundled together. A "currants-in-jelly" model is developed that provides a flexible approach in which the contribution of the soft high-level objects (jelly-like) are employed to protect the underlying atomic structure (currants), while still allowing them to interact. Idealised chains were used to establish the parameters to achieve this degree of interaction over a hierarchy spanning four levels and in a more realistic example, the distortion experienced by a protein domain structure during collision was measured and the parameters refined. This model of steric repulsion was then combined with sets of predicted distance constraints, derived from correlated mutation analysis. Firstly, an integral trans-membrane protein was modelled in which the packing of the seven helices was refined but without topological rearrangement. Secondly, an RNA structure was 'folded' under the predicted constraints, starting only from its 2-dimensional secondary structure prediction.
Project description:The structural and mechanical properties of monomeric actin (G-actin), the trimer nucleus, and actin filaments (F-actins) are determined as a function of the conformation of the DNase I-binding loop (DB loop) by using all-atom molecular dynamics simulations and coarse-grained (CG) analysis. Recent x-ray structures of ADP-bound G-actin (G-ADP) by Otterbein et al. [Otterbein, L. R., Graceffa, P. & Dominguez, R. (2001) Science 293, 708-711] and ATP-bound G-actin (G-ATP) by Graceffa and Dominguez [Graceffa, P. & Dominguez, R. (2003) J. Biol. Chem. 278, 34172-34180] indicate that the DB loop of actin does not have a well defined secondary structure in the ATP state but folds into an alpha-helix in the ADP state. MD simulations and CG analysis indicate that such a helical DB loop significantly weakens the intermonomer interactions of actin assemblies and thus leads to a wider, shorter, and more disordered filament. The computed persistence lengths of F-actin composed of G-ATP (16 microm) and of G-ADP (8.5 microm) agree well with the experimental values for the two states. Therefore, the loop-to-helix transition of the DB loop may be one of the factors that lead to the changes in structural and mechanical properties of F-actin after ATP hydrolysis. This result may provide a direct connection between the conformational changes of an actin monomer and the structural and mechanical properties of the cytoskeleton. The information provided by MD simulations also helps to understand the possible origin of the special features of actin dynamics.
Project description:Maculatin 1.1 (M1.1) is a membrane-active antimicrobial peptide (AMP) from an Australian tree frog that forms a kinked amphipathic alpha-helix in the presence of a lipid bilayer or bilayer-mimetic environment. To help elucidate its mechanism of membrane-lytic activity, we performed a total of approximately 8 micros of coarse-grained molecular dynamics (CG-MD) simulations of M1.1 in the presence of zwitterionic phospholipid membranes. Several systems were simulated in which the peptide/lipid ratio was varied. At a low peptide/lipid ratio, M1.1 adopted a kinked, membrane-interfacial location, consistent with experiment. At higher peptide/lipid ratios, we observed spontaneous, cooperative membrane insertion of M1.1 peptide aggregates. The minimum size for formation of a transmembrane (TM) aggregate was just four peptides. The absence of a simple and well-defined central channel, along with the exclusion of lipid headgroups from the aggregates, suggests that a pore-like model is an unlikely explanation for the mechanism of membrane lysis by M1.1. We also performed an extended 1.25 micros simulation of the permeabilization of a complete liposome by multiple peptides. Consistent with the simpler bilayer simulations, formation of monomeric interfacial peptides and TM peptide clusters was observed. In contrast, major structural changes were observed in the vesicle membrane, implicating induced membrane curvature in the mechanism of active antimicrobial peptide lysis. This contrasted with the behavior of the nonpore-forming model peptide WALP23, which inserted into the vesicle to form extended clusters of TM alpha-helices with relatively little perturbation of bilayer properties.