Browse
Submit Data
Databases
API
Help

Dataset Information

15 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway.

ABSTRACT: Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold's mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

SUBMITTER: Zhang P

PROVIDER: S-EPMC7732773 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Unfolding simulations reveal the mechanism of extreme unfolding cooperativity in the kinetically stable alpha-lytic protease.

Project description:Kinetically stable proteins, those whose stability is derived from their slow unfolding kinetics and not thermodynamics, are examples of evolution's best attempts at suppressing unfolding. Especially in highly proteolytic environments, both partially and fully unfolded proteins face potential inactivation through degradation and/or aggregation, hence, slowing unfolding can greatly extend a protein's functional lifetime. The prokaryotic serine protease alpha-lytic protease (alphaLP) has done just that, as its unfolding is both very slow (t(1/2) approximately 1 year) and so cooperative that partial unfolding is negligible, providing a functional advantage over its thermodynamically stable homologs, such as trypsin. Previous studies have identified regions of the domain interface as critical to alphaLP unfolding, though a complete description of the unfolding pathway is missing. In order to identify the alphaLP unfolding pathway and the mechanism for its extreme cooperativity, we performed high temperature molecular dynamics unfolding simulations of both alphaLP and trypsin. The simulated alphaLP unfolding pathway produces a robust transition state ensemble consistent with prior biochemical experiments and clearly shows that unfolding proceeds through a preferential disruption of the domain interface. Through a novel method of calculating unfolding cooperativity, we show that alphaLP unfolds extremely cooperatively while trypsin unfolds gradually. Finally, by examining the behavior of both domain interfaces, we propose a model for the differential unfolding cooperativity of alphaLP and trypsin involving three key regions that differ between the kinetically stable and thermodynamically stable classes of serine proteases.

| S-EPMC2829044 | biostudies-literature

Transmembrane Pore Structures of ?-Hairpin Antimicrobial Peptides by All-Atom Simulations.

Project description:Protegrin-1 is an 18-residue ?-hairpin antimicrobial peptide (AMP) that has been suggested to form transmembrane ?-barrels in biological membranes. However, alternative structures have also been proposed. Here, we performed multimicrosecond, all-atom molecular dynamics simulations of various protegrin-1 oligomers on the membrane surface and in transmembrane topologies. The membrane surface simulations indicated that protegrin dimers are stable, whereas trimers and tetramers break down. Tetrameric arcs remained stably inserted in lipid membranes, but the pore water was displaced by lipid molecules. Unsheared protegrin ?-barrels opened into ?-sheets that surrounded stable aqueous pores, whereas tilted barrels with sheared hydrogen bonding patterns were stable in most topologies. A third type of observed pore consisted of multiple small oligomers surrounding a small, partially lipidic pore. We also considered the ?-hairpin AMP tachyplesin, which showed less tendency to oligomerize than protegrin: the octameric bundle resulted in small pores surrounded by six peptides as monomers and dimers, with some peptides returning to the membrane surface. The results imply that multiple configurations of protegrin oligomers may produce aqueous pores and illustrate the relationship between topology and putative steps in protegrin-1's pore formation. However, the long-term stability of these structures needs to be assessed further.

| S-EPMC5686775 | biostudies-literature

All-Atom Simulations Reveal How Single-Point Mutations Promote Serpin Misfolding.

Project description:Protein misfolding is implicated in many diseases, including serpinopathies. For the canonical inhibitory serpin ?1-antitrypsin, mutations can result in protein deficiencies leading to lung disease, and misfolded mutants can accumulate in hepatocytes, leading to liver disease. Using all-atom simulations based on the recently developed bias functional algorithm, we elucidate how wild-type ?1-antitrypsin folds and how the disease-associated S (Glu264Val) and Z (Glu342Lys) mutations lead to misfolding. The deleterious Z mutation disrupts folding at an early stage, whereas the relatively benign S mutant shows late-stage minor misfolding. A number of suppressor mutations ameliorate the effects of the Z mutation, and simulations on these mutants help to elucidate the relative roles of steric clashes and electrostatic interactions in Z misfolding. These results demonstrate a striking correlation between atomistic events and disease severity and shine light on the mechanisms driving chains away from their correct folding routes.

| S-EPMC5961751 | biostudies-literature

All-Atom MD Simulations of the HBV Capsid Complexed with AT130 Reveal Secondary and Tertiary Structural Changes and Mechanisms of Allostery.

Project description:The hepatitis B virus (HBV) capsid is an attractive drug target, relevant to combating viral hepatitis as a major public health concern. Among small molecules known to interfere with capsid assembly, the phenylpropenamides, including AT130, represent an important antiviral paradigm based on disrupting the timing of genome packaging. Here, all-atom molecular dynamics simulations of an intact AT130-bound HBV capsid reveal that the compound increases spike flexibility and improves recovery of helical secondary structure in the spike tips. Regions of the capsid-incorporated dimer that undergo correlated motion correspond to established sub-domains that pivot around the central chassis. AT130 alters patterns of correlated motion and other essential dynamics. A new conformational state of the dimer is identified, which can lead to dramatic opening of the intradimer interface and disruption of communication within the spike tip. A novel salt bridge is also discovered, which can mediate contact between the spike tip and fulcrum even in closed conformations, revealing a mechanism of direct communication across these sub-domains. Altogether, results describe a dynamical connection between the intra- and interdimer interfaces and enable mapping of allostery traversing the entire core protein dimer.

| S-EPMC8065791 | biostudies-literature

All-Atom Simulations Reveal a Key Interaction Network in the HLA-E/NKG2A/CD94 Immune Complex Fine-Tuned by the Nonameric Peptide.

Project description:Natural killer (NK) cells, an important part of the innate immune system, can clear a wide variety of pathological challenges, including tumor, senescent, and virally infected cells. They express various activating and inhibitory receptors on their surface, and the balance of interactions between them and specific ligands displayed on the surface of target cells is critical for NK cell cytolytic function and target cell protection. The CD94/NKG2A heterodimer is one of the inhibitory receptors that interacts with its trimeric ligand consisting of HLA-E, β2m, and a nonameric peptide. Here, multi-microsecond-long all-atom molecular dynamics simulations of eight immune complexes elucidate the subtleties of receptor (NKG2A/CD94)-ligand (HLA-E/β2m/peptide) molecular recognition that mediate the NK cell protection from a geometric and energetic perspective. We identify key differences in the interactions between the receptor and ligand complexes, which are via an entangled network of hydrogen bonds fine-tuned by the ligand-specific nonameric peptide. We further reveal that the receptor protein NKG2A regulates the NK cell activity, while its CD94 partner forms the majority of the energetically important interactions with the ligand. This knowledge rationalizes the atomistic details of the fundamental NK cell protection mechanism and may enable a variety of opportunities in rational-based drug discovery for diverse pathologies including viral infections and cancer and elimination of senescent cells associated with potential treatment of many age-related diseases.

| S-EPMC8389527 | biostudies-literature

Stability tests on known and misfolded structures with discrete and all atom molecular dynamics simulations.

Project description:The brevity of molecular dynamics simulations often limits their utility in developing and evaluating structural models of proteins. The duration of simulations can be increased greatly using discrete molecular dynamics (DMD). However, the trade off is that coarse graining, implicit solvent, and other time-saving procedures reduce the accuracy of DMD simulations. Here we address some of these issues by comparing results of DMD and conventional all atom MD simulations on proteins of known structure and misfolded proteins. DMD simulations were performed at a range of temperatures to identify a 'physiological' temperature for DMD that mimicked molecular motions of conventional MD simulations at 310K. We also compared results obtained with a new implicit solvent model developed here based on Miyazawa-Jernigan interaction pair potential to those obtained with a previously used model based on Kyte-Doolittle hydropathy scale. We compared DMD and all atom molecular dynamics with explicit water by simulating both correctly and incorrectly folded structures, and monomeric and dimeric α β-barrel structures to analyze the ability of these procedures to distinguish between good and bad models. Deviations from the correct structures were substantially greater with DMD, as would be expected from coarse-graining and longer simulation time. Deviations were smallest for β-strands and greatest for coiled loops. Structures of the incorrectly folded models were very poorly preserved during the DMD simulations; but both methods were able to distinguish between the correct and the incorrect structures based on differences in the magnitudes of the root mean squared deviation (RMSD) from the starting conformation.

| S-EPMC3104411 | biostudies-literature

Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics.

Project description:As theory and experiment have shown, protein dehydration is a major contributor to protein folding. Dehydration upon folding can be characterized directly by all-atom simulations of fast pressure drops, which create desolvated pockets inside the nascent hydrophobic core. Here, we study pressure-drop refolding of three ?-repressor fragment (?6-85) mutants computationally and experimentally. The three mutants report on tertiary structure formation via different fluorescent helix-helix contact pairs. All-atom simulations of pressure drops capture refolding and unfolding of all three mutants by a similar mechanism, thus validating the nonperturbative nature of the fluorescent contact probes. Analysis of simulated interprobe distances shows that the ?-helix 1-3 pair distance displays a slower characteristic time scale than the 1-2 or 3-2 pair distance. To see whether slow packing of ?-helices 1 and 3 is reflected in the rate-limiting folding step, fast pressure-drop relaxation experiments captured refolding on a millisecond time scale. These experiments reveal that refolding monitored by 1-3 contact formation indeed is much slower than when monitored by 1-2 or 3-2 contact formation. Unlike the case of the two-state folder [three-?-helix bundle (?3D)], whose drying and core formation proceed in concert, ?6-85 repeatedly dries and rewets different local tertiary contacts before finally forming a solvent-excluded core, explaining the non-two-state behavior observed during refolding in molecular dynamics simulations. This work demonstrates that proteins can explore desolvated pockets and dry globular states numerous times before reaching the native conformation.

| S-EPMC6431218 | biostudies-literature

Structures of asymmetric ClpX hexamers reveal nucleotide-dependent motions in a AAA+ protein-unfolding machine.

Project description:ClpX is a AAA+ machine that uses the energy of ATP binding and hydrolysis to unfold native proteins and translocate unfolded polypeptides into the ClpP peptidase. The crystal structures presented here reveal striking asymmetry in ring hexamers of nucleotide-free and nucleotide-bound ClpX. Asymmetry arises from large changes in rotation between the large and small AAA+ domains of individual subunits. These differences prevent nucleotide binding to two subunits, generate a staggered arrangement of ClpX subunits and pore loops around the hexameric ring, and provide a mechanism for coupling conformational changes caused by ATP binding or hydrolysis in one subunit to flexing motions of the entire ring. Our structures explain numerous solution studies of ClpX function, predict mechanisms for pore elasticity during translocation of irregular polypeptides, and suggest how repetitive conformational changes might be coupled to mechanical work during the ATPase cycle of ClpX and related molecular machines.

| S-EPMC2778613 | biostudies-literature

Project description:Nuclear spin relaxation is a powerful method for studying molecular dynamics at atomic resolution. Recent methods development in biomolecular NMR spectroscopy has enabled detailed investigations of molecular dynamics that are critical for biological function, with prominent examples addressing allostery, enzyme catalysis, and protein folding. Dynamic processes with similar correlation times are often detected in multiple locations of the molecule, raising the question of whether the underlying motions are correlated (corresponding to concerted fluctuations involving many atoms distributed across extended regions of the molecule) or uncorrelated (corresponding to independent fluctuations involving few atoms in localized regions). Here, we have used (13)C(alpha)(i - 1)/(13)C(alpha)(i) differential multiple-quantum spin relaxation to provide direct evidence for correlated dynamics of consecutive amino acid residues in the protein sequence. By monitoring overlapping pairs of residues (i - 1 and i, i and i + 1, etc.), we identified correlated motions that extend through continuous segments of the sequence. We detected significant correlated conformational transitions in the native state of the E140Q mutant of the calmodulin C-terminal domain. Previous work has shown that this domain exchanges between two major conformational states that resemble the functionally relevant open and closed states of the WT protein, with a mean correlation time of approximately 20 micros. The present results reveal that an entire alpha-helix undergoes partial unraveling in a transient and cooperative manner.

| S-EPMC1287973 | biostudies-literature

Million-atom molecular dynamics simulations reveal the interfacial interactions and assembly of plant PSII-LHCII supercomplex.

Project description:Protein-protein interface interactions dictate efficient excitation energy transfer from light-harvesting antennas to the photosystem II (PSII) core. In this work, we construct a 1.2 million atom-scale model of plant C2S2-type PSII-LHCII supercomplex and perform microsecond-scale molecular dynamics (MD) simulations to explore the interactions and assembly mechanisms of the sizeable PSII-LHCII supercomplex. We optimize the nonbonding interactions of the PSII-LHCII cryo-EM structure using microsecond-scale MD simulations. Binding free energy calculations with component decompositions reveal that hydrophobic interactions predominantly drive antenna-core association and the antenna-antenna interactions are relatively weak. Despite the positive electrostatic interaction energies, hydrogen bonds and salt bridges mainly provide directional or anchoring forces for interface binding. Analysis of the roles of small intrinsic subunits of PSII suggests that LHCII and CP26 first interact with small intrinsic subunits and then bind to the core proteins, whereas CP29 adopts a one-step binding process to the PSII core without the assistance of other factors. Our study provides insights into the molecular underpinnings of the self-organization and regulation of plant PSII-LHCII. It lays the framework for deciphering the general assembly principles of photosynthetic supercomplexes and possibly other macromolecular structures. The finding also has implications for repurposing photosynthetic systems to enhance photosynthesis.

| S-EPMC9969236 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data