Interactions between Inhibitors and 5-Lipoxygenase: Insights from Gaussian Accelerated Molecular Dynamics and Markov State Models

Inflammation is a protective stress response triggered by external stimuli, with 5-lipoxygenase (5LOX) playing a pivotal role as a potent mediator of the leukotriene (Lts) inflammatory pathway. Nordihydroguaiaretic acid (NDGA) functions as a natural orthosteric inhibitor of 5LOX, while 3-acetyl-11-keto-β-boswellic acid (AKBA) acts as a natural allosteric inhibitor targeting 5LOX. However, the precise mechanisms of inhibition have remained unclear. In this study, Gaussian accelerated molecular dynamics (GaMD) simulation was employed to elucidate the inhibitory mechanisms of NDGA and AKBA on 5LOX. It was found that the orthosteric inhibitor NDGA was tightly bound in the protein's active pocket, occupying the active site and inhibiting the catalytic activity of the 5LOX enzyme through competitive inhibition. The binding of the allosteric inhibitor AKBA induced significant changes at the distal active site, leading to a conformational shift of residues 168-173 from a loop to an α-helix and significant negative correlated motions between residues 285-290 and 375-400, reducing the distance between these segments. In the simulation, the volume of the active cavity in the stable conformation of the protein was reduced, hindering the substrate's entry into the active cavity and, thereby, inhibiting protein activity through allosteric effects. Ultimately, Markov state models (MSM) were used to identify and classify the metastable states of proteins, revealing the transition times between different conformational states. In summary, this study provides theoretical insights into the inhibition mechanisms of 5LOX by AKBA and NDGA, offering new perspectives for the development of novel inhibitors specifically targeting 5LOX, with potential implications for anti-inflammatory drug development.

Revealing the Role of the Arg and Lys in Shifting Paradigm from BTK Selective Inhibition to the BTK/HCK Dual Inhibition - Delving into the Inhibitory Activity of KIN-8194 against BTK, and HCK in the Treatment of Mutated BTKCys481 Waldenström Macroglobulinemia: A Computational Approach

BACKGROUND

Despite the early success of Bruton's tyrosine kinase (BTK) inhibitors in the treatment of Waldenström macroglobulinemia (WM), these single-target drug therapies have limitations in their clinical applications, such as drug resistance. Several alternative strategies have been developed, including the use of dual inhibitors, to maximize the therapeutic potential of these drugs.

OBJECTIVE

Recently, the pharmacological activity of KIN-8194 was repurposed to serve as a 'dual-target' inhibitor of BTK and Hematopoietic Cell Kinase (HCK). However, the structural dual inhibitory mechanism remains unexplored, hence the aim of this study.

METHODS

Conducting predictive pharmacokinetic profiling of KIN-8194, as well as demonstrating a comparative structural mechanism of inhibition against the above-mentioned enzymes.

RESULTS

Our results revealed favourable binding affinities of -20.17 kcal/mol, and -35.82 kcal/mol for KIN-8194 towards HCK and BTK, respectively. Catalytic residues Arg137/174 and Lys42/170 in BTK and Arg303 and Lys75/173/244/247 in HCK were identified as crucial mediators of the dual binding mechanism of KIN-8194, corroborated by high per-residue energy contributions and consistent high-affinity interactions of these residues. Prediction of the pharmacokinetics and physicochemical properties of KIN-8194 further established its inhibitory potential, evidenced by the favourable absorption, metabolism, excretion, and minimal toxicity properties. Structurally, KIN-8194 impacted the stability, flexibility, solvent-accessible surface area, and rigidity of BTK and HCK, characterized by various alterations observed in the bound and unbound structures, which proved enough to disrupt their biological function.

CONCLUSION

These structural insights provided a baseline for the understanding of the dual inhibitory activity of KIN- 8194. Establishing the cruciality of the interactions between the KIN-8194 and Arg and Lys residues could guide the structure-based design of novel dual BTK/HCK inhibitors with improved therapeutic activities.

Defective ORF8 dimerization in SARS-CoV-2 delta variant leads to a better adaptive immune response due to abrogation of ORF8-MHC1 interaction

In India, during the second wave of the COVID-19 pandemic, the breakthrough infections were mainly caused by the SARS-COV-2 delta variant (B.1.617.2). It was reported that, among majority of the infections due to the delta variant, only 9.8% percent cases required hospitalization, whereas only 0.4% fatality was observed. Sudden dropdown in COVID-19 infections cases were observed within a short timeframe, suggesting better host adaptation with evolved delta variant. Downregulation of host immune response against SARS-CoV-2 by ORF8 induced MHC-I degradation has been reported earlier. The Delta variant carried mutations (deletion) at Asp119 and Phe120 amino acids which are critical for ORF8 dimerization. The deletions of amino acids Asp119 and Phe120 in ORF8 of delta variant resulted in structural instability of ORF8 dimer caused by disruption of hydrogen bonds and salt bridges as revealed by structural analysis and MD simulation studies. Further, flexible docking of wild type and mutant ORF8 dimer revealed reduced interaction of mutant ORF8 dimer with MHC-I as compared to wild-type ORF8 dimer with MHC-1, thus implicating its possible role in MHC-I expression and host immune response against SARS-CoV-2. We thus propose that mutant ORF8 of SARS-CoV-2 delta variant may not be hindering the MHC-I expression thereby resulting in a better immune response against the SARS-CoV-2 delta variant, which partly explains the possible reason for sudden drop of SARS-CoV-2 infection rate in the second wave of SARS-CoV-2 predominated by delta variant in India.

Slower diffusion and anomalous association of R453W lamin A protein alter nuclear architecture in AD-EDMD

Lamins maintain the shape and rigidity of the nucleus in the form of a proteinaceous scaffold underneath the inner nuclear membrane (INM) and provide anchorage to chromatin and other nuclear proteins. Mutations in the human LMNA gene encoding lamin A/C cause about 16 different diseases with distinct phenotypes collectively termed as laminopathies which affect primarily the muscle tissues as well as adipose tissues, neuromuscular junctions and multiple other organs in progeroid syndromes. Lamins contain several domains of which Ig-fold is one of the well characterized and structured domains that harbours many mutations leading to deleterious interactions with other nuclear proteins. In this work, we have elucidated the effects of 3 such mutations namely R453W, W498C and W498R on the dynamics and flexibility of the Ig-fold domain and the consequent effect on the assembly into lamina by live cell imaging, fluorescence correlation spectroscopy (FCS) and molecular dynamics (MD) simulations. From our simulation studies, we concluded that R453W exhibits the highest fluctuation at the residues 475 and 525 in the Ig fold domain compared to the wild type and other mutants. This resulted in pronounced random self-association which could be corroborated by lower diffusivity values obtained from FCS. This is the first report where such an alteration in the full length has been documented by gross changes in diffusional properties as a sequel to a mutation in the Ig fold domain.

Evaluation of immune evasion in SARS-CoV-2 Delta and Omicron variants

Emerging SARS-CoV-2 variants with higher transmissibility and immune escape remain a persistent threat across the globe. This is evident from the recent outbreaks of the Delta (B.1.617.2) and Omicron variants. These variants have originated from different continents and spread across the globe. In this study, we explored the genomic and structural basis of these variants for their lineage defining mutations of the spike protein through computational analysis, protein modeling, and molecular dynamic (MD) simulations. We further experimentally validated the importance of these deletion mutants for their immune escape using a pseudovirus-based neutralization assay, and an antibody (4A8) that binds directly to the spike protein's NTD. Delta variant with the deletion and mutations in the NTD revealed a better rigidity and reduced flexibility as compared to the wild-type spike protein (Wuhan isolate). Furthermore, computational studies of 4A8 monoclonal antibody (mAb) revealed a reduced binding of Delta variant compared to the wild-type strain. Similarly, the MD simulation data and virus neutralization assays revealed that the Omicron also exhibits immune escape, as antigenic beta-sheets appear to be disrupted. The results of the present study demonstrate the higher possibility of immune escape and thereby achieved better fitness advantages by the Delta and Omicron variants, which warrants further demonstrations through experimental evidences. Our study, based on in-silico computational modelling, simulations, and pseudovirus-based neutralization assay, highlighted and identified the probable mechanism through which the Delta and Omicron variants are more pathogenically evolved with higher transmissibility as compared to the wild-type strain.

Co-Binding of JQ1 and Venetoclax Exhibited Synergetic Inhibitory Effect for Cancer Therapy; Potential Line of Treatment for the Waldenström Macroglobulinemia Lymphoma

In recent times, the development of combination therapy has been a focal point in drug discovery. This article explores the potential synergistic effect of co-administration of Bcl2 inhibitor Venetoclax and BET inhibitor JQ1. We envisioned that the 'dual-site'-binding of Bcl2 has significant prospects and paves the way for the next round of rational design of potent Waldenström macroglobulinemia (WM) therapy. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. This study conducted a systematic procedure to explore the preferred binding modes and the structural effects of co-binding at each catalytic active site. Interestingly, a mutual enhanced binding effect was observed - Venetoclax increased the binding affinity of JQ1 by 11.5 %, while JQ1 boosted the binding affinity of Venetoclax by 16.3 % when compared with individual inhibition of each drug. This synergistic binding effect has significantly increased protein stability, with substantial correlated movements and multiple van der Waals interactions. The structural and thermodynamic insights unveiled in this report would assist the future design of improved combined therapy against WM.

Multi-catalytic Sites Inhibition of Bcl2 Induces Expanding of Hydrophobic Groove: A New Avenue Towards Waldenström Macroglobulinemia Therapy

B-cell lymphoma 2 (Bcl2) is a key protein regulator of apoptosis. The hydrophobic groove in Bcl2 is a unique structural feature to this class of enzymes and found to have a profound impact on protein overall structure, function, and dynamics. Dynamics of the hydrophobic groove is an essential determinant of the catalytic activity of Bcl2, an implicated protein in Waldenström macroglobulinemia (WM). The mobility of α3-α4 helices around the catalytic site of the protein remains crucial to its activity. The preferential binding mechanisms of the multi-catalytic sites of the Bcl2 enzyme have been a subject of debate in the literature. In addition to our previous report on the same protein, herein, we further investigate the preferential binding modes and the conformational implications of Venetoclax-JQ1 dual drug binding at both catalytic active sites of Bcl2. Structural analysis revealed asymmetric α3-α4 helices movement with the expansion of the distance between the α3 and α4 helix in Venetoclax-JQ1 dual inhibition by 15.2% and 26.3%, respectively when compared to JQ1 and Venetoclax individual drug inhibition-resulting in remarkable widening of hydrophobic groove. Moreso, a reciprocal enhanced binding effect was observed: Venetoclax increased the binding affinity of JQ1 by 11.5%, while the JQ1 fostered the binding affinity of Venetoclax by 16.3% compared with individual inhibition of each drug. This divergence has also resulted in higher protein stability, and prominent correlated motions were observed with the least fluctuations and multiple van der Waals interactions. Findings offer vital conformational dynamics and structural mechanisms of enzyme-single ligand and enzyme-dual ligand interactions, which could potentially shift the current therapeutic protocol of Waldenström macroglobulinemia.

Impact of single amino acid substitution on the structure and function of TANK-binding kinase-1

Tank-binding kinase 1 (TBK1) is a serine/threonine protein kinase involved in various signaling pathways and subsequently regulates cell proliferation, apoptosis, autophagy, antiviral and antitumor immunity. Dysfunction of TBK1 can cause many complex diseases, including autoimmunity, neurodegeneration, and cancer. This dysfunction of TBK1 may result from single amino acid substitutions and subsequent structural alterations. This study analyzed the effect of substituting amino acids on TBK1 structure, function, and subsequent disease using advanced computational methods and various tools. In the initial assessment, a total of 467 mutations were found to be deleterious. After that, in detailed structural and sequential analyses, 13 mutations were found to be pathogenic. Finally, based on the functional importance, two variants (K38D and S172A) of the TBK1 kinase domain were selected and studied in detail by utilizing all-atom molecular dynamics (MD) simulation for 200 ns. MD simulation, including correlation matrix and principal component analysis, helps to get deeper insights into the TBK1 structure at the atomic level. We observed a substantial change in variants' conformation, which may be possible for structural alteration and subsequent TBK1 dysfunction. However, substitution S172A shows a significant conformational change in TBK1 structure as compared to K38D. Thus, this study provides a structural basis to understand the effect of mutations on the kinase domain of TBK1 and its function associated with disease progression.

Impact of non-synonymous mutations on the structure and function of telomeric repeat binding factor 1

Telomeric repeat binding factor 1 (TRF1) is one of the major components of the shelterin complex. It directly binds to the telomere and controls its function by regulating the telomerase acting on it. Several variations are reported in the TRF1 gene; some are associated with variety of diseases. Here, we have studied the structural and functional significance of these variations in the TRFH domain of TRF1. We have used cutting-edge computational methods such as SIFT, PolyPhen-2, PROVEAN, Mutation Assessor, mCSM, SDM, STRUM, MAESTRO, and DUET to predict the effects of 124 mutations in the TRFH domain of TRF1. Out of 124 mutations, we have identified 12 deleterious mutations with high confidence based on their prediction. To see the impact of the finally selected mutations on the structure and stability of TRF1, all-atom molecular dynamics (MD) simulations on TRF1-Wild type (WT), L79R and P150R mutants for 200 ns were carried out. A significant conformational change in the structure of the P150R mutant was observed. Our integrated computational study provides a comprehensive understanding of structural changes in TRF1 incurred due to the mutations and subsequent function, leading to the progression of many diseases.Communicated by Ramaswamy H. Sarma.

ABBV-744 as a potential inhibitor of SARS-CoV-2 main protease enzyme against COVID-19

A new pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide and become pandemic with thousands new deaths and infected cases globally. To address coronavirus disease (COVID-19), currently no effective drug or vaccine is available. This necessity motivated us to explore potential lead compounds by considering drug repurposing approach targeting main protease (M^pro) enzyme of SARS-CoV-2. This enzyme considered to be an attractive drug target as it contributes significantly in mediating viral replication and transcription. Herein, comprehensive computational investigations were performed to identify potential inhibitors of SARS-CoV-2 M^pro enzyme. The structure-based pharmacophore modeling was developed based on the co-crystallized structure of the enzyme with its biological active inhibitor. The generated hypotheses were applied for virtual screening based PhaseScore. Docking based virtual screening workflow was used to generate hit compounds using HTVS, SP and XP based Glide GScore. The pharmacological and physicochemical properties of the selected lead compounds were characterized using ADMET. Molecular dynamics simulations were performed to explore the binding affinities of the considered lead compounds. Binding energies revealed that compound ABBV-744 binds to the M^pro with strong affinity (ΔG_bind -45.43 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. Our study classified three best compounds which could be considered as promising inhibitors against main protease SARS-CoV-2 virus.

Probing inhibition mechanisms of adenosine deaminase by using molecular dynamics simulations

Adenosine deaminase (ADA) catalyzes the deamination of adenosine, which is important in purine metabolism. ADA is ubiquitous to almost all human tissues, and ADA abnormalities have been reported in various diseases, including rheumatoid arthritis. ADA can be divided into two conformations based on the inhibitor that it binds to: open and closed forms. Here, we chose three ligands, namely, FR117016 (FR0), FR221647 (FR2) (open form), and HDPR (PRH, closed form), to investigate the inhibition mechanism of ADA and its effect on ADA through molecular dynamics simulations. In open forms, Egap and electrostatic potential (ESP) indicated that electron transfer might occur more easily in FR0 than in FR2. Binding free energy and hydrogen bond occupation revealed that the ADA-FR0 complex had a more stable structure than ADA-FR2. The probability of residues Pro159 to Lys171 of ADA-FR0 and ADA-FR2 to form a helix moderately increased compared with that in nonligated ADA. In comparison with FR0 and FR2 PRH could maintain ADA in a closed form to inhibit the function of ADA. The α7 helix (residues Thr57 to Ala73) of ADA in the closed form was mostly unfastened because of the effect of PRH. The number of H bonds and the relative superiority of the binding free energy indicated that the binding strength of PRH to ADA was significantly lower than that of an open inhibitor, thereby supporting the comparison of the inhibitory activities of the three ligands. Alanine scanning results showed that His17, Gly184, Asp295, and Asp296 exerted the greatest effects on protein energy, suggesting that they played crucial roles in binding to inhibitors. This study served as a theoretical basis for the development of new ADA inhibitors.

Molecular mechanism of substrate selectivity of the arginine-agmatine Antiporter AdiC

The arginine-agmatine antiporter (AdiC) is a component of an acid resistance system developed by enteric bacteria to resist gastric acidity. In order to avoid neutral proton antiport, the monovalent form of arginine, about as abundant as its divalent form under acidic conditions, should be selectively bound by AdiC for transport into the cytosol. In this study, we shed light on the mechanism through which AdiC distinguishes Arg⁺ from Arg²⁺ of arginine by investigating the binding of both forms in addition to that of divalent agmatine, using a combination of molecular dynamics simulations with molecular and quantum mechanics calculations. We show that AdiC indeed preferentially binds Arg⁺. The weaker binding of divalent compounds results mostly from their greater tendency to remain hydrated than Arg⁺. Our data suggests that the binding of Arg⁺ promotes the deprotonation of Glu208, a gating residue, which in turn reinforces its interactions with AdiC, leading to longer residence times of Arg⁺ in the binding site. Although the total electric charge of the ligand appears to be the determinant factor in the discrimination process, two local interactions formed with Trp293, another gating residue of the binding site, also contribute to the selection mechanism: a cation-π interaction with the guanidinium group of Arg⁺ and an anion-π interaction involving Glu208.

Understanding the interactions of different substrates with wild-type and mutant acylaminoacyl peptidase using molecular dynamics simulations

Acylaminoacylpeptidase (AAP) belongs to peptidase protein family, which can degrade amyloid β-peptide forms in the brains of patients, and hence leads to Alzheimer's disease. And so, AAP is considered to be a novel target in the design of drugs against Alzheimer's disease. In this investigation, six molecular dynamics simulations were used to find that the interaction between the wild-type and R526V AAP with two different substrates (p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide). Our results were as follows: firstly, Ac-Leu-p-nitroanilide bound to R526V AAP to form a more disordered loop (residues 552-562) in the α/β-hydrolase fold like of AAP, which caused an open and inactive AAP domain form, secondly, binding p-nitrophenylcaprylate and Ac-Leu-p-nitroanilide to AAP can decrease the flexibility of residues 225-250, 260-270, and 425-450, in which the ordered secondary structures may contain the suitable geometrical structure and so it is useful to serine attack. Our theoretical results showed that the binding of the two substrates can induce specific conformational changes responsible for the diverse AAP catalytic specificity. These theoretical substrate-induced structural diversities can help explain the abilities of AAPs to recognize and hydrolyze extremely different substrates.

Exploring the molecular basis of RNA recognition by the dimeric RNA-binding protein via molecular simulation methods

RNA-binding protein with multiple splicing (RBPMS) is critical for axon guidance, smooth muscle plasticity, and regulation of cancer cell proliferation and migration. Recently, different states of the RNA-recognition motif (RRM) of RBPMS, one in its free form and another in complex with CAC-containing RNA, were determined by X-ray crystallography. In this article, the free RRM domain, its wild type complex and 2 mutant complex systems are studied by molecular dynamics (MD) simulations. Through comparison of free RRM domain and complex systems, it's found that the RNA binding facilitates stabilizing the RNA-binding interface of RRM domain, especially the C-terminal loop. Although both R38Q and T103A/K104A mutations reduce the binding affinity of RRM domain and RNA, the underlining mechanisms are different. Principal component analysis (PCA) and Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) methods were used to explore the dynamical and recognition mechanisms of RRM domain and RNA. R38Q mutation is positioned on the homodimerization interface and mainly induces the large fluctuations of RRM domains. This mutation does not directly act on the RNA-binding interface, but some interfacial hydrogen bonds are weakened. In contrast, T103A/K104A mutations are located on the RNA-binding interface of RRM domain. These mutations obviously break most of high occupancy hydrogen bonds in the RNA-binding interface. Meanwhile, the key interfacial residues lose their favorable energy contributions upon RNA binding. The ranking of calculated binding energies in 3 complex systems is well consistent with that of experimental binding affinities. These results will be helpful in understanding the RNA recognition mechanisms of RRM domain.

Unveiling the Mechanism of Arginine Transport through AdiC with Molecular Dynamics Simulations: The Guiding Role of Aromatic Residues

Commensal and pathogenic enteric bacteria have developed several systems to adapt to proton leakage into the cytoplasm resulting from extreme acidic conditions. One such system involves arginine uptake followed by export of the decarboxylated product agmatine, carried out by the arginine/agmatine antiporter (AdiC), which thus works as a virtual proton pump. Here, using classical and targeted molecular dynamics, we investigated at the atomic level the mechanism of arginine transport through AdiC of E. coli. Overall, our MD simulation data clearly demonstrate that global rearrangements of several transmembrane segments are necessary but not sufficient for achieving transitions between structural states along the arginine translocation pathway. In particular, local structural changes, namely rotameric conversions of two aromatic residues, are needed to regulate access to both the outward- and inward-facing states. Our simulations have also enabled identification of a few residues, overwhelmingly aromatic, which are essential to guiding arginine in the course of its translocation. Most of them belong to gating elements whose coordinated motions contribute to the alternating access mechanism. Their conservation in all known E. coli acid resistance antiporters suggests that the transport mechanisms of these systems share common features. Last but not least, knowledge of the functional properties of AdiC can advance our understanding of the members of the amino acid-carbocation-polyamine superfamily, notably in eukaryotic cells.

Exploring the different ligand escape pathways in acylaminoacyl peptidase by random acceleration and steered molecular dynamics simulations

Acylaminoacyl peptidase (APH, EC 3.4.19.1) is a novel class of serine-type protease belonging to the prolyl oligopeptidase (POP) family.

Molecular dynamics simulations of wild type and mutants of botulinum neurotoxin A complexed with synaptic vesicle protein 2C

Using molecular dynamics simulations, we investigate the relationship between the conformational changes of BoNT/A-RBD:SV2C-LD and the interfacial interactions.

Structural and functional implications of the yeast high-affinity tryptophan permease Tat2

Tryptophan is hydrophobic, bulky, and the rarest amino acid found in nutrients. Accordingly, the import machinery can be specialized evolutionarily. Our previous study in Saccharomyces cerevisiae demonstrated that tryptophan import by the high-affinity tryptophan permease Tat2 is accompanied by a large volume increase during substrate import. Nevertheless, the mechanisms by which the permease mediates tryptophan recognition and permeation remain to be elucidated. Here we determined amino acid residues essential for Tat2-mediated tryptophan import. By means of random mutagenesis in combination with site-directed mutagenesis based on crystallographic studies of the Escherichia coli arginine/agmatine antiporter AdiC, we identified 15 amino acid residues in the Tat2 transmembrane domains (TMDs) 1, -3, -5, -8, and -10, which are responsible for tryptophan uptake. T98, Y167, and E286 were assumed to form the central cavity in Tat2. G97/T98 and E286 were located within the putative α-helix break in TMD1 and TMD6, respectively, which are highly conserved among yeast amino acid permeases and bacterial solute transporters. Given the conformational change in AdiC upon substrate binding, G97/T98 and E286 of Tat2 were assumed to mediate a structural shift from an outward-open to a tryptophan-bound-occluded structure upon tryptophan binding, and T320, V322, and F324 became stabilized in TMD7. Such dynamic structural changes may account for the large volume increase associated with tryptophan import occurring concomitantly with a movement of water molecules from the tryptophan binding site. We also propose the working hypothesis that E286 mediates the proton influx that is coupled to tryptophan import.

Network models reveal stability and structural rearrangement of signal recognition particle

The signal recognition particle (SRP) and its receptors (SR) mediate the cotranslational targeting of the membrane and secretory proteins in all cells. In Escherichia coli, SRP is composed of the Ffh protein and the 4.5S SRP RNA. Ffh is a multidomain protein comprising a methionine-rich (M) domain, a helical N domain, and a Ras-like guanine triphosphatase (GTPase) (G) domain. The N and G domains are commonly referred to as one structural unit, the NG domain. In this article, the complex structure of SRP and SR is investigated with the Gaussian network model (GNM) and anisotropic network model (ANM). GNM provides the information of structure stability. It is found that the intermolecular interactions between SRP and SR can obviously decrease the fluctuation of NG domains. Nevertheless, the large structural rearrangement will take place during the cotranslational protein targeting cycle. Hence, the moving directions of fluctuation regions are further ascertained by using cross-correlation analysis and the ANM. The NG domain of Ffh undergoes a clockwise rotation around the GM linker and the M domain of Ffh shows an opposite direction to the NG domain. These functional movements will facilitate the SRP structure to transform into the free form and the sequence-bound form. These simple coarse-grained analyses can be used as a general and quick method for the mechanism studies of protein assembly and supramolecular systems.

Computational model for protein unfolding simulation

The protein folding problem is one of the fundamental and important questions in molecular biology. However, the all-atom molecular dynamics studies of protein folding and unfolding are still computationally expensive and severely limited by the time scale of simulation. In this paper, a simple and fast protein unfolding method is proposed based on the conformational stability analyses and structure modeling. In this method, two structure-based conditions are considered to identify the unstable regions of proteins during the unfolding processes. The protein unfolding trajectories are mimicked through iterative structure modeling according to conformational stability analyses. Two proteins, chymotrypsin inhibitor 2 (CI2) and α -spectrin SH3 domain (SH3) were simulated by this method. Their unfolding pathways are consistent with the previous molecular dynamics simulations. Furthermore, the transition states of the two proteins were identified in unfolding processes and the theoretical Φ values of these transition states showed significant correlations with the experimental data (the correlation coefficients are >0.8). The results indicate that this method is effective in studying protein unfolding. Moreover, we analyzed and discussed the influence of parameters on the unfolding simulation. This simple coarse-grained model may provide a general and fast approach for the mechanism studies of protein folding.