Impact of Mutations on NPAC Structural Dynamics: Mechanistic Insights from MD Simulations

Deficiency of nucleosome-destabilizing factor GLYR1 dampens spermatogenesis in mice

Aberrant sperm morphology hinders sperm motility and causes male subfertility. Spermatogenesis, a complex process in male germ cell development, necessitates precise regulation of numerous developmental genes. However, the regulatory pathways involved in this process remain partially understood. We have observed the widespread expression of Glyr1, the gene encoding a nucleosome-destabilizing factor, in mouse testicular cells. Our study demonstrates that mice experiencing Glyr1 depletion in spermatogenic cells exhibit subfertility characterized by a diminished count and motility of spermatozoa. Furthermore, the rate of sperm malformation significantly increases in the absence of Glyr1, with a predominant occurrence of head and neck malformation in spermatozoa within the cauda epididymis. Additionally, a reduction in spermatocyte numbers across different meiotic stages is observed, accompanied by diminished histone acetylation in spermatogenic cells upon Glyr1 depletion. Our findings underscore the crucial roles of Glyr1 in mouse spermiogenesis and unveil novel insights into the etiology of male reproductive diseases.

The NPAC-LSD2 complex in nucleosome demethylation

NPAC is a transcriptional co-activator widely associated with the H3K36me3 epigenetic marks present in the gene bodies. NPAC plays a fundamental role in RNA polymerase progression, and its depletion downregulates gene transcription. In this chapter, we review the current knowledge on the functional and structural features of this multi-domain protein. NPAC (also named GLYR1 or NP60) contains a PWWP motif, a chromatin binder and epigenetic reader that is proposed to weaken the DNA-histone contacts facilitating polymerase passage through the nucleosomes. The C-terminus of NPAC is a catalytically inactive dehydrogenase domain that forms a stable and rigid tetramer acting as an oligomerization module for the formation of co-transcriptional multimeric complexes. The PWWP and dehydrogenase domains are connected by a long, mostly disordered, linker that comprises putative sites for protein and DNA interactions. A short dodecapeptide sequence (residues 214-225) forms the binding site for LSD2, a flavin-dependent lysine-specific histone demethylase. This stretch of residues binds on the surface of LSD2 and facilitates the capture and processing of the H3 tail in the nucleosome context, thus promoting the H3K4me1/2 epigenetic mark removal. LSD2 is associated with other two chromatin modifiers, G9a and NSD3. The LSD2-G9a-NSD3 complex modifies the pattern of the post translational modifications deposited on histones, thus converting the relaxed chromatin into a transcriptionally refractory state after the RNA polymerase passage. NPAC is a scaffolding factor that organizes and coordinates the epigenetic activities required for optimal transcription elongation.

A Computational Biology Study on the Structure and Dynamics Determinants of Thermal Stability of the Chitosanase from Aspergillus fumigatus

Environmentally friendly and efficient biodegradation with chitosanase for degrading chitosan to oligosaccharide has been gaining more importance. Here, we studied a chitosanase from Aspergillus fumigatus with potential for production, but does not have the ideal thermal stability. The structure predicted by the Alphafold2 model, especially the binding site and two catalytic residues, has been found to have a high similarity with the experimental structure of the chitosanase V-CSN from the same family. The effects of temperature on structure and function were studied by dynamic simulation and the results showed that the binding site had high flexibility. After heating up from 300 K to 350 K, the RMSD and RMSF of the binding site increased significantly, in particular, the downward shift of loop6 closed the binding site, resulting in the spatial hindrance of binding. The time proportions of important hydrogen bonds at the binding site decreased sharply, indicating that serious disruption of hydrogen bonds should be the main interaction factor for conformational changes. The residues contributing energetically to binding were also revealed to be in the highly flexible region, which inevitably leads to the decrease in the activity stability at high temperature. These findings provide directions for the modification of thermal stability and perspectives on the research of proteins without experimental structures.

Computational Epitope Prediction and Design for Antibody Development and Detection

The design of optimized protein antigens is a fundamental step in the development of new vaccine candidates and in the detection of therapeutic antibodies. A fundamental prerequisite is the identification of antigenic regions that are most prone to interact with antibodies, namely, B-cell epitopes. Here, we describe an efficient structure-based computational method for epitope prediction, called MLCE. In this approach, all that is required is the 3D structure of the antigen of interest. MLCE can be applied to glycosylated proteins, facilitating the identification of immunoreactive versus immune-shielding carbohydrates.

Unlocking SGK1 inhibitor potential of bis-[1-N,7-N, pyrazolo tetraethoxyphthalimido{-4-(3,5-Dimethyl-4-(spiro-3-methylpyazolo)-1,7-dihydro-1H-dipyrazolo[3,4-b;4',3'-e]pyridin-8-yl)}]p-disubstituted phenyl compounds: a computational study

SGK1 (Serum and Glucocorticoid Regulated Kinase 1), a serine/threonine kinase that is activated by various stimuli, including serum and glucocorticoids. It controls inflammation, apoptosis, hormone release, neuro-excitability and cell proliferation, all of which play an important role in cancer progression and metastasis. SGK1 was recently proposed as a potential drug target for cancer, diabetes, and neurodegenerative diseases. In this study, molecular docking, physiochemical, toxicological properties and molecular dynamic simulation of the Bis-[1-N,7-N, Pyrazolo tetraethoxyphthalimido{-4-(3,5-Dimethyl-4-(spiro-3-methylpyazolo)-1,7-dihydro-1H-dipyrazolo[3,4-b;4',3'-e]pyridin-8-yl)}]p-disubstituted phenyl compoundsand reference EMD638683 against new SGK1 target protein. Compared to the reference inhibitor EMD638683, we choose the best compounds (series 2-6) based on the binding energy (in the range from -11.0 to -10.6 kcal/mol). With the exception of compounds 2 and 6, none of the compounds posed a risk for AMES toxicity or carcinogenicity due to their toxicological properties. 100 ns MD simulation accompanied by MM/PBSA energy calculations and PCA. According to MD simulation results, the binding of compounds 3, 4 and 5 stabilizes the SGK1 structure and causes febrile conformational changes compared to EMD638683. As a result of this research, the final selected compounds 3, 4 and 5 can be used as scaffolds to develop promising SGK1 inhibitors for the treatment of related diseases such as cancer.Communicated by Ramaswamy H. Sarma.

Human prion protein: exploring the thermodynamic stability and structural dynamics of its pathogenic mutants

Human familial prion diseases are known to be associated with different single-point mutants of the gene coding for prion protein with a primary focus at several locations of the globular domain. We have identified 12 different single-point pathogenic mutants of human prion protein (HuPrP) with the help of extensive perturbations/mutation technique at multiple locations of HuPrP sequence related to potentiality towards conformational disorders. Among these, some of the mutants include pathogenic variants that corroborate well with the literature reported proteins while majority include some unique single-point mutants that are either not explicitly studied early or studied for variants with different residues at the specific position. Primarily, our study sheds light on the unfolding mechanism of the above mentioned mutants in depth. Besides, we could identify some mutants under investigation that demonstrates not only unfolding of the helical structures but also extension and generation of the β-sheet structures and or simultaneously have highly exposed hydrophobic surface which is assumed to be linked with the production of aggregate/fibril structures of the prion protein. Among the identified mutants, Q212E needs special attention due to its maximum exposure of hydrophobic core towards solvent and E200Q is found to be important due to its maximum extent of β-content. We are also able to identify different respective structural conformations of the proteins according to their degree of structural unfolding and those conformations can be extracted and further studied in detail. Communicated by Ramaswamy H. Sarma.

The Mitochondrial HSP90 Paralog TRAP1: Structural Dynamics, Interactome, Role in Metabolic Regulation, and Inhibitors

The HSP90 paralog TRAP1 was discovered more than 20 years ago; yet, a detailed understanding of the function of this mitochondrial molecular chaperone remains elusive. The dispensable nature of TRAP1 in vitro and in vivo further complicates an understanding of its role in mitochondrial biology. TRAP1 is more homologous to the bacterial HSP90, HtpG, than to eukaryotic HSP90. Lacking co-chaperones, the unique structural features of TRAP1 likely regulate its temperature-sensitive ATPase activity and shed light on the alternative mechanisms driving the chaperone’s nucleotide-dependent cycle in a defined environment whose physiological temperature approaches 50 °C. TRAP1 appears to be an important bioregulator of mitochondrial respiration, mediating the balance between oxidative phosphorylation and glycolysis, while at the same time promoting mitochondrial homeostasis and displaying cytoprotective activity. Inactivation/loss of TRAP1 has been observed in several neurodegenerative diseases while TRAP1 expression is reported to be elevated in multiple cancers and, as with HSP90, evidence of addiction to TRAP1 has been observed. In this review, we summarize what is currently known about this unique HSP90 paralog and why a better understanding of TRAP1 structure, function, and regulation is likely to enhance our understanding of the mechanistic basis of mitochondrial homeostasis.

"In silico identification of ethoxy phthalimide pyrazole derivatives as IL-17A and IL-18 targeted gouty arthritis agents"

Two proinflammatory cytokines, IL17A and IL18, are observed to be elevated in the serum of gout patients and they play a crucial role in the development and worsening of inflammation, which has severe effects. In present study, we have combined molecular docking, molecular dynamics studies and MM-PBSA analysis to study the effectiveness of ethoxy phthalimide pyrazole derivatives (series 3a to 3e) as potential inhibitors against cytokines IL17A and IL18 as a druggable targets. The binding energy of the docked series ranges from -13.5 to -10.0 kcal/mol and extensively interacts with the amino acids in the active pocket of IL17A and IL18. Compound 3e had the lowest binding energy with IL17A at -12.6 kcal/mol compared to control allopurinol (3.32 kcal/mol). With IL18, compound 3a seems to have the lowest binding energy of -9.6 kcal/mol compared to control allopurinol (3.18 kcal/mol). In MD simulation studies, compound 3a forms a stable and energetically stabilized complex with the target protein. Depending on properties of the bound IL17A-3a and IL18-3a complexes was compared by means of MM-PBSA analysis. These derivatives can be used as a scaffold to develop promising IL17A and IL18 inhibitors to assess their potential for gouty arthritis and other related diseases.

Transcription factor protein interactomes reveal genetic determinants in heart disease

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

Targeting a cryptic allosteric site of SIRT6 with small-molecule inhibitors that inhibit the migration of pancreatic cancer cells

SIRT6 belongs to the conserved NAD⁺-dependent deacetylase superfamily and mediates multiple biological and pathological processes. Targeting SIRT6 by allosteric modulators represents a novel direction for therapeutics, which can overcome the selectivity problem caused by the structural similarity of orthosteric sites among deacetylases. Here, developing a reversed allosteric strategy AlloReverse, we identified a cryptic allosteric site, Pocket Z, which was only induced by the bi-directional allosteric signal triggered upon orthosteric binding of NAD⁺. Based on Pocket Z, we discovered an SIRT6 allosteric inhibitor named JYQ-42. JYQ-42 selectively targets SIRT6 among other histone deacetylases and effectively inhibits SIRT6 deacetylation, with an IC₅₀ of 2.33 μmol/L. JYQ-42 significantly suppresses SIRT6-mediated cancer cell migration and pro-inflammatory cytokine production. JYQ-42, to our knowledge, is the most potent and selective allosteric SIRT6 inhibitor. This study provides a novel strategy for allosteric drug design and will help in the challenging development of therapeutic agents that can selectively bind SIRT6.

Elucidating the 3D Structure of a Surface Membrane Antigen from Trypanosoma cruzi as a Serodiagnostic Biomarker of Chagas Disease

Chagas disease (CD) is a vector-borne parasitosis, caused by the protozoan parasite Trypanosoma cruzi, that affects millions of people worldwide. Although endemic in South America, CD is emerging throughout the world due to climate change and increased immigratory flux of infected people to non-endemic regions. Containing of the diffusion of CD is challenged by the asymptomatic nature of the disease in early infection stages and by the lack of a rapid and effective diagnostic test. With the aim of designing new serodiagnostic molecules to be implemented in a microarray-based diagnostic set-up for early screening of CD, herein, we report the recombinant production of the extracellular domain of a surface membrane antigen from T. cruzi (TcSMP) and confirm its ability to detect plasma antibodies from infected patients. Moreover, we describe its high-resolution (1.62 Å) crystal structure, to which in silico epitope predictions were applied in order to locate the most immunoreactive regions of TcSMP in order to guide the design of epitopes that may be used as an alternative to the full-length antigen for CD diagnosis. Two putative, linear epitopes, belonging to the same immunogenic region, were synthesized as free peptides, and their immunological properties were tested in vitro. Although both peptides were shown to adopt a structural conformation that allowed their recognition by polyclonal antibodies raised against the recombinant protein, they were not serodiagnostic for T. cruzi infections. Nevertheless, they represent good starting points for further iterative structure-based (re)design cycles.

Structure-function relationships of the disease-linked A218T oxytocin receptor variant

Various single nucleotide polymorphisms (SNPs) in the oxytocin receptor (OXTR) gene have been associated with behavioral traits, autism spectrum disorder (ASD) and other diseases. The non-synonymous SNP rs4686302 results in the OXTR variant A218T and has been linked to core characteristics of ASD, trait empathy and preterm birth. However, the molecular and intracellular mechanisms underlying those associations are still elusive. Here, we uncovered the molecular and intracellular consequences of this mutation that may affect the psychological or behavioral outcome of oxytocin (OXT)-treatment regimens in clinical studies, and provide a mechanistic explanation for an altered receptor function. We created two monoclonal HEK293 cell lines, stably expressing either the wild-type or A218T OXTR. We detected an increased OXTR protein stability, accompanied by a shift in Ca²⁺ dynamics and reduced MAPK pathway activation in the A218T cells. Combined whole-genome and RNA sequencing analyses in OXT-treated cells revealed 7823 differentially regulated genes in A218T compared to wild-type cells, including 429 genes being associated with ASD. Furthermore, computational modeling provided a molecular basis for the observed change in OXTR stability suggesting that the OXTR mutation affects downstream events by altering receptor activation and signaling, in agreement with our in vitro results. In summary, our study provides the cellular mechanism that links the OXTR rs4686302 SNP with genetic dysregulations associated with aspects of ASD.

Molecular dynamics and protein frustration analysis of human fused in Sarcoma protein variants in Amyotrophic Lateral Sclerosis type 6: An In Silico approach

Amyotrophic lateral sclerosis (ALS) is the most frequent adult-onset motor neuron disorder. The disease is characterized by degeneration of upper and lower motor neurons, leading to death usually within five years after the onset of symptoms. While most cases are sporadic, 5%-10% of cases can be associated with familial inheritance, including ALS type 6, which is associated with mutations in the Fused in Sarcoma (FUS) gene. This work aimed to evaluate how the most frequent ALS-related mutations in FUS, R521C, R521H, and P525L affect the protein structure and function. We used prediction algorithms to analyze the effects of the non-synonymous single nucleotide polymorphisms and performed evolutionary conservation analysis, protein frustration analysis, and molecular dynamics simulations. Most of the prediction algorithms classified the three mutations as deleterious. All three mutations were predicted to reduce protein stability, especially the mutation R521C, which was also predicted to increase chaperone binding tendency. The protein frustration analysis showed an increase in frustration in the interactions involving the mutated residue 521C. Evolutionary conservation analysis showed that residues 521 and 525 of human FUS are highly conserved sites. The molecular dynamics results indicate that protein stability could be compromised in all three mutations. They also affected the exposed surface area and protein compactness. The analyzed mutations also displayed high flexibility in most residues in all variants, most notably in the interaction site with the nuclear import protein of FUS.

Discovery of Novel HCV NS5B polymerase inhibitor, 2-(3,4-dimethyl-5,5-dioxidobenzo[e]pyrazolo[4,3-c][1,2]thiazin-2(4H)-yl)-N-(2-fluorobenzyl)acetamide via molecular docking and experimental approach

Hepatitis C Virus (HCV) is a viral infection posing a severe global threat that left untreated progress to end-stage liver disease, including cirrhosis and Hepatocellular Carcinoma (HCC). Moreover, no prophylactic approach exists so far enabling its prevention. The NS5B polymerase holds special significance as the target of intervention against HCV infection. The current study kindles benzothiazine derivatives against HCV NS5B polymerase through in silico and experimental approaches. Following docking, the compound 2-(3,4-dimethyl-5,5-dioxidobenzo[e]pyrazolo[4,3-c][1,2]thiazin-2(4H)-yl)-N-(2-fluorobenzyl)acetamide was revealed to form effective binding interaction in the proposed site of HCV NS5B with a score of -10 kcal/mol and subsequently was deciphered through MD simulation study which indicated interaction of residues TYR_382, VAL_381 and HIS_467 through hydrophobic interaction and two residues such as GLU_202 and LYS_209 contributed in the formation of water bridges. The subsequent in silico pharmacological analysis revealed its safe drug profile. The cytotoxicity activity of compound 6c indicated to be non-toxic in HepG2 cells at concentration ranges from 0.001-1.0 µM with > 80% cell viability and diminished expression of the HCV NS5B to 98% at the dose of 1.0 µM and 90 % at 0.5µM. Thus the hit compound 6c might be a potent NS5B polymerase inhibitor required to be validated further through in vivo and preclinical studies.

The molecular chaperone TRAP1 in cancer: From the basics of biology to pharmacological targeting

TRAP1, the mitochondrial component of the Hsp90 family of molecular chaperones, displays important bioenergetic and proteostatic functions. In tumor cells, TRAP1 contributes to shape metabolism, dynamically tuning it with the changing environmental conditions, and to shield from noxious insults. Hence, TRAP1 activity has profound effects on the capability of neoplastic cells to evolve towards more malignant phenotypes. Here, we discuss our knowledge on the biochemical functions of TRAP1 in the context of a growing tumor mass, and we analyze the possibility of targeting its chaperone functions for developing novel anti-neoplastic approaches.

Computational evidence of new putative allosteric sites in the acetylcholinesterase receptor

Acetylcholinesterase (AChE), with a rigid structure and buried active site at the end of a deep narrow gorge, is interesting enough to solve the paradox between high catalytic activity and unavailability of the active site in treatment of Alzheimer's disease (AD). In this way, the blind docking process is performed on an ensemble of AChE structures created with molecular dynamics (MD) simulations to survey the whole space of AChE to find multiple access pathways to the active site and ranking them based on their affinity scores. Our results show that there are other allosteric binding sites in the protein structure whose inhibition, can affect protein function by disrupting the release of the Acetylcholine (AC) degradation products. In this study, inhibitory activities of Hybride14 and two natural compounds (Papaverine and Palmatine) were evaluated for all possible allosteric sites via docking method. The results confirmed the non-competitive inhibition mechanism. The best binding mode for these inhibitors and efficacy of hydrogen bonds and hydrophobic interactions on inhibitory activities of ligands were also disclosed. Furthermore, our studies provide significant molecular insight for AChE inhibition that could aid in the development of new drugs for AD's treatment.

Unraveling the Allosteric Cross-Talk between the Coactivator Peptide and the Ligand-Binding Site in the Glucocorticoid Receptor

The glucocorticoid receptor (GR) is a nuclear receptor that controls critical biological processes by regulating the transcription of specific genes. There is a known allosteric cross-talk between the ligand and coregulator binding sites within the GR ligand-binding domain that is crucial for the control of the functional response. However, the molecular mechanisms underlying such an allosteric control remain elusive. Here, molecular dynamics (MD) simulations, bioinformatic analysis, and biophysical measurements are integrated to capture the structural and dynamic features of the allosteric cross-talk within the GR. We identified a network of evolutionarily conserved residues that enables the allosteric signal transduction, in agreement with experimental data. MD simulations clarify how such a network is dynamically interconnected and offer a mechanistic explanation of how different peptides affect the intensity of the allosteric signal. This study provides useful insights to elucidate the GR allosteric regulation, ultimately providing a foundation for designing novel drugs.

A Comparative Study to Decipher the Structural and Dynamics Determinants Underlying the Activity and Thermal Stability of GH-11 Xylanases

With the growing need for renewable sources of energy, the interest for enzymes capable of biomass degradation has been increasing. In this paper, we consider two different xylanases from the GH-11 family: the particularly active GH-11 xylanase from Neocallimastix patriciarum, NpXyn11A, and the hyper-thermostable mutant of the environmentally isolated GH-11 xylanase, EvXyn11^TS. Our aim is to identify the molecular determinants underlying the enhanced capacities of these two enzymes to ultimately graft the abilities of one on the other. Molecular dynamics simulations of the respective free-enzymes and enzyme–xylohexaose complexes were carried out at temperatures of 300, 340, and 500 K. An in-depth analysis of these MD simulations showed how differences in dynamics influence the activity and stability of these two enzymes and allowed us to study and understand in greater depth the molecular and structural basis of these two systems. In light of the results presented in this paper, the thumb region and the larger substrate binding cleft of NpXyn11A seem to play a major role on the activity of this enzyme. Its lower thermal stability may instead be caused by the higher flexibility of certain regions located further from the active site. Regions such as the N-ter, the loops located in the fingers region, the palm loop, and the helix loop seem to be less stable than in the hyper-thermostable EvXyn11^TS. By identifying molecular regions that are critical for the stability of these enzymes, this study allowed us to identify promising targets for engineering GH-11 xylanases. Eventually, we identify NpXyn11A as the ideal host for grafting the thermostabilizing traits of EvXyn11^TS.

Targeting the HSP90-CDC37-kinase chaperone cycle: A promising therapeutic strategy for cancer

Heat shock protein 90 (HSP90) is an indispensable molecular chaperone that facilitates the maturation of numerous oncoproteins in cancer cells, including protein kinases, ribonucleoproteins, steroid hormone receptors, and transcription factors. Although over 30 HSP90 inhibitors have steadily entered clinical trials, further clinical advancement has been restricted by their limited efficacy, inevitable heat shock response, and multiple side-effects, likely induced via an ATP inhibition mechanism. Since both ATP and various co-chaperones play essential roles in the HSP90 chaperone cycle to achieve integrated function, optimal therapeutics require an understanding of the dynamic interactions among HSP90, ATP, and cochaperones. To date, continuous research has promoted the exploration of the cochaperone cell division cycle 37 (CDC37) as a kinase-specific recognizer and has shown that the HSP90-CDC37-kinase complex is particularly relevant in cancers. Indeed, disrupting the HSP90-CDC37-kinase complex, rather than totally blocking the ATP function of HSP90, is emerging as an alternative way to avoid the limitations of current inhibitors. In this review, we first briefly introduce the HSP90-CDC37-kinase cycle and present the currently available approaches for inhibitor development targeting this cycle and provide insights into selective regulation of the kinase clients of HSP90 by more directional ways.

Exploiting Folding and Degradation Machineries To Target Undruggable Proteins: What Can a Computational Approach Tell Us?

Advances in genomics and proteomics have unveiled an ever-growing number of key proteins and provided mechanistic insights into the genesis of pathologies. This wealth of data showed that changes in expression levels of specific proteins, mutations, and post-translational modifications can result in (often subtle) perturbations of functional protein-protein interaction networks, which ultimately determine disease phenotypes. Although many such validated pathogenic proteins have emerged as ideal drug targets, there are also several that escape traditional pharmacological regulation; these proteins have thus been labeled "undruggable". The challenges posed by undruggable targets call for new sorts of molecular intervention. One fascinating solution is to perturb a pathogenic protein's expression levels, rather than blocking its activities. In this Concept paper, we shall discuss chemical interventions aimed at recruiting undruggable proteins to the ubiquitin proteasome system, or aimed at disrupting protein-protein interactions in the chaperone-mediated cellular folding machinery: both kinds of intervention lead to a decrease in the amount of active pathogenic protein expressed. Specifically, we shall discuss the role of computational strategies in understanding the molecular determinants characterizing the function of synthetic molecules typically designed for either type of intervention. Finally, we shall provide our perspectives and views on the current limitations and possibilities to expand the scope of rational approaches to the design of chemical regulators of protein levels.

Molecular Dynamics Reveals a DNA-Induced Dynamic Switch Triggering Activation of CRISPR-Cas12a

CRISPR-Cas12a is a genome-editing system, recently also harnessed for nucleic acid detection, which is promising for the diagnosis of the SARS-CoV-2 coronavirus through the DETECTR technology. Here, a collective ensemble of multimicrosecond molecular dynamics characterizes the key dynamic determinants allowing nucleic acid processing in CRISPR-Cas12a. We show that DNA binding induces a switch in the conformational dynamics of Cas12a, which results in the activation of the peripheral REC2 and Nuc domains to enable cleavage of nucleic acids. The simulations reveal that large-amplitude motions of the Nuc domain could favor the conformational activation of the system toward DNA cleavages. In this process, the REC lobe plays a critical role. Accordingly, the joint dynamics of REC and Nuc shows the tendency to prime the conformational transition of the DNA target strand toward the catalytic site. Most notably, the highly coupled dynamics of the REC2 region and Nuc domain suggests that REC2 could act as a regulator of the Nuc function, similar to what was observed previously for the HNH domain in the CRISPR-associated nuclease Cas9. These mutual domain dynamics could be critical for the nonspecific binding of DNA and thereby for the underlying mechanistic functioning of the DETECTR technology. Considering that REC is a key determinant in the system's specificity, our findings provide a rational basis for future biophysical studies aimed at characterizing its function in CRISPR-Cas12a. Overall, our outcomes advance our mechanistic understanding of CRISPR-Cas12a and provide grounds for novel engineering efforts to improve genome editing and viral detection.

Markov state models and NMR uncover an overlooked allosteric loop in p53

The tumor suppressor p53 is the most frequently mutated gene in human cancer, and thus reactivation of mutated p53 is a promising avenue for cancer therapy. Analysis of wildtype p53 and the Y220C cancer mutant long-timescale molecular dynamics simulations with Markov state models and validation by NMR relaxation studies has uncovered the involvement of loop L6 in the slowest motions of the protein. Due to its distant location from the DNA-binding surface, the conformational dynamics of this loop has so far remained largely unexplored. We observe mutation-induced stabilization of alternate L6 conformations, distinct from all experimentally-determined structures, in which the loop is both extended and located further away from the DNA-interacting surface. Additionally, the effect of the L6-adjacent Y220C mutation on the conformational landscape of the functionally-important loop L1 suggests an allosteric role to this dynamic loop and the inactivation mechanism of the mutation. Finally, the simulations reveal a novel Y220C cryptic pocket that can be targeted for p53 rescue efforts. Our approach exemplifies the power of the MSM methodology for uncovering intrinsic dynamic and kinetic differences among distinct protein ensembles, such as for the investigation of mutation effects on protein function.

Wildtype and Y220C L1 and L6 loops conformational landscape, with MSM-identified L6 states highlighted on the right.

Simple Model of Protein Energetics To Identify Ab Initio Folding Transitions from All-Atom MD Simulations of Proteins

A fundamental requirement to predict the native conformation, address questions of sequence design and optimization, and gain insights into the folding mechanisms of proteins lies in the definition of an unbiased reaction coordinate that reports on the folding state without the need to compare it to reference values, which might be unavailable for new (designed) sequences. Here, we introduce such a reaction coordinate, which does not depend on previous structural knowledge of the native state but relies solely on the energy partition within the protein: the spectral gap of the pair nonbonded energy matrix (ENergy Gap, ENG). This quantity can be simply calculated along unbiased MD trajectories. We show that upon folding the gap increases significantly, while its fluctuations are reduced to a minimum. This is consistently observed for a diverse set of systems and trajectories. Our approach allows one to promptly identify residues that belong to the folding core as well as residues involved in non-native contacts that need to be disrupted to guide polypeptides to the folded state. The energy gap and fluctuations criteria are then used to develop an automatic detection system which allows us to extract and analyze folding transitions from a generic MD trajectory. We speculate that our method can be used to detect conformational ensembles in dynamic and intrinsically disordered proteins, revealing potential preorganization for binding.

Probing intrinsic dynamics and conformational transition of HIV gp120 by molecular dynamics simulation

The HIV envelope glycoprotein gp120 has evolved two distinct conformational states to balance viral infection and immune escape. One is a closed state resistant to most neutralization antibodies, and the other is an open state responsible for the binding of the receptor and coreceptors. Although the structures of gp120 in these two conformational states have been determined, a detailed molecular mechanism involving intrinsic dynamics and conformational transition is still elusive. In this study, μs-scale molecular dynamics simulation is performed to probe molecular dynamics and conformational transition away from the open state and approach the closed state. Our results reveal that open gp120 shows a larger structural deviation, higher conformational flexibility, and more conformational diversity than the form in the closed state, providing a structural explanation for receptor or coreceptor affinity at the open state and the neutralization resistance of closed conformation. Seven regions with greatly decreased coupled motions in the open states have been observed by dynamic cross-correlation analysis, indicating that conformational transition can be mainly attributed to the relaxation of intrinsic dynamics. Three conformations characterized by the structural orientations of the V1/V2 region and the V3 loop, suggesting gp120 is intrinsically dynamic from the open state to the closed state. Taken together, these findings shed light on the understanding of the conformational control mechanism of HIV.

The HIV envelope glycoprotein gp120 has evolved two distinct conformational states to balance viral infection and immune escape.

mTOR–mLST8 interaction: hot spot identification through quantum biochemistry calculations

Quantum calculation of mTOR–mLST8 interaction.