Free Energy Landscape Analysis Research Articles

How Does the Mono-Triazole Derivative Modulate Aβ42 Aggregation and Disrupt a Protofibril Structure: Insights from Molecular Dynamics Simulations.

Clinical studies have identified that abnormal self-assembly of amyloid-β (Aβ) peptide into toxic fibrillar aggregates is associated with the pathology of Alzheimer’s disease (AD). The most acceptable therapeutic approach to stop the progression of AD is to inhibit the formation of β-sheet-rich structures. Recently, we designed and evaluated a series of novel mono-triazole derivatives 4(a–x), where compound 4v was identified as the most potent inhibitor of Aβ42 aggregation and disaggregates preformed Aβ42 fibrils significantly. Moreover, 4v strongly averts the Cu2+-induced Aβ42 aggregation and disaggregates the preformed Cu2+-induced Aβ42 fibrils, halts the generation of reactive oxygen species, and shows neuroprotective effects in SH-SY5Y cells. However, the underlying molecular mechanism of inhibition of Aβ42 aggregation by 4v and disaggregation of preformed Aβ42 fibrils remains obscure. In this work, molecular dynamics (MD) simulations have been performed to explore the conformational ensemble of the Aβ42 monomer and a pentameric protofibril structure of Aβ42 in the presence of 4v. The MD simulations highlighted that 4v binds preferentially at the central hydrophobic core region of the Aβ42 monomer and chains D and E of the Aβ42 protofibril. The dictionary of secondary structure of proteins analysis indicated that 4v retards the conformational conversion of the helix-rich structure of the Aβ42 monomer into the aggregation-prone β-sheet conformation. The binding free energy calculated by the molecular mechanics Poisson–Boltzmann surface area method revealed an energetically favorable process with ΔGbinding = −44.9 ± 3.3 kcal/mol for the Aβ42 monomer–4v complex. The free energy landscape analysis highlighted that the Aβ42 monomer–4v complex sampled conformations with significantly higher helical contents (35 and 49%) as compared to the Aβ42 monomer alone (17%). Compound 4v displayed hydrogen bonding with Gly37 (chain E) and π–π interactions with Phe19 (chain D) of the Aβ42 protofibril. Further, the per-residue binding free energy analysis also highlighted that Phe19 (chain D) and Gly37 (chain E) of the Aβ42 protofibril showed the maximum contribution in the binding free energy. The decreased binding affinity and residue–residue contacts between chains D and E of the Aβ42 protofibril in the presence of 4v indicate destabilization of the Aβ42 protofibril structure. Overall, the structural information obtained through MD simulations indicated that 4v stabilizes the native helical conformation of the Aβ42 monomer and persuades a destabilization in the protofibril structure of Aβ42. The results of the study will be useful in the rational design of potent inhibitors against amyloid aggregation.

Bioinformatic and computational analysis for predominant mutations of the Nrf2/Keap1 complex in pediatric leukemia

The levels of reactive oxygen species (ROS) are tightly controlled and regulated by Nuclear Factor Erythroid-2-Like 2 (Nrf2) transcription factor, which is the main regulator of antioxidant responses and its suppressor protein Kelch-like ECH-associated protein 1 (Keap1). Our previous study has identified six novel changes in Nrf2/Keap1 pathway in pediatric ALL, which were described for the first time. These changes in the pathway are likely to alter the evolutionary process of amino acids and cause structural changes in the final products of genes. In this study, we aimed to compare the pathogenicity of eight determined mutations reported in our previous study by utilizing different programs with different algorithms and molecular dynamics simulation. Since it is too difficult to handle each existing mutation in a wet laboratory, in silico methods may give suggestion to choose the important mutations for further analysis and to establish the appropriate patient population and conduct wet laboratory studies. For this purpose, four different algorithms were used to evaluate the effects of single amino acid mutation. In addition, root-mean-square deviation, root-mean-square fluctuation and free-energy landscape analyses were performed to observe stability, flexibility and energetically favorable conformations, respectively, for each amino acid mutation. As a result, our study emphasizes the importance of Keap1 mutations in pediatric ALL Nrf2/Keap1 pathway, a total of eight mutations, two of which were shown for the first time in our study. Especially the mutations in the Keap1 Broad-Complex, Tramtrack and Bric-à-brac domain are worthy of attention. Communicated by Ramaswamy H. Sarma.

Mechanistic insights into the deleterious roles of Nasu-Hakola disease associated TREM2 variants

Recently, the critical roles played by genetic variants of TREM2 (Triggering Receptor Expressed on Myeloid cells 2) in Alzheimer’s disease have been aggressively highlighted. However, few studies have focused on the deleterious roles of Nasu-Hakola disease (NHD) associated TREM2 variants. In order to get insights into the contributions made by these variants to neurodegeneration, we investigated the influences of four NHD associated TREM2 mutations (Y38C, W50C, T66M, and V126G) on loss-of-function, and followed this with in silico prediction and conventional molecular dynamics simulation. NHD mutations were predicted to be highly deleterious by eight different in silico bioinformatics tools and found to induce conformational changes by molecular dynamics simulation. As compared with the wild-type, the four variants produced substantial differences in the collective motions of loop regions, which not only promoted structural remodeling in the CDR2 (complementarity-determining region 2) loop but also in the CDR1 loop, by changing inter- and intra-loop hydrogen bonding networks. In addition, structural studies in a free energy landscape analysis showed that Y38, T66, and V126 are crucial for maintaining the structural features of CDR1 and CDR2 loops, and that mutations in these positions produced steric clashes and loss of ligand binding. These results showed the presence of mutations in the TREM2 ectodomain induced flexibility and caused structural alterations. Dynamical scenarios, as provided by the present study, may be critical to our understanding of the roles of these TREM2 mutations in neurodegenerative diseases.

Effect of Piedmont mutation (L34V) on the structure, dynamics, and aggregation of Alzheimer’s Aβ40 peptide

The amyloid-β (Aβ) aggregation in the brain has been associated with the development of Alzheimer’s disease (AD). The previous studies have reported that Piedmont mutation (L34V) increases the rate of Aβ40 aggregation. However, the underlying molecular mechanism of the effect of L34V mutation on Aβ40 structure, dynamics, and aggregation remains largely unclear. In the present study, molecular dynamics (MD) simulations were performed to elucidate the effect of L34V mutation on the structural changes and conformational dynamics of Aβ40. The secondary structure analysis highlight that L34V mutation enhances Aβ40 self-assembly due to the formation of aggregation-prone β–sheet structure at the C-terminus of Aβ40 monomeric structure. The higher probability of Asp23–Lys28 salt bridge interaction in Aβ40(L34V) leads to aggregation prone β–sheet conformations, which has the potential to increase the fibril formation rate. The free energy landscape (FEL) analysis depict a sampling of coil conformation in the free energy minima of Aβ40, whereas the aggregation-prone β–sheet conformation was observed at the C-terminal region of Aβ40(L34V) in the minimum energy conformations extracted from FEL of Aβ40(L34V). MD simulations, in agreement with experiment, highlight that L34V mutation increases Aβ40 aggregation as the sampling of the aggregation–prone β–sheet conformation substantially increased. Overall, MD simulations provided atomic level details into the increased fibril formation tendency upon L34V mutation and physical insights into the L34V-mediated conformational as well as structural changes in Aβ40.

Delineating the folding perturbations and molecular mechanisms of Thr-Ala 642 mutation in Rab-GTPase activating protein Akt substrate of 160kDa and its impact on the aetiology of diabetes

Diabetes Mellitus is a complex metabolic disorder with one of the highest prevalence rates in the world. The present study probes into the Thr-Ala 642 mutation of Akt substrate of 160 kDa (AS160) which has been implicated in diabetes by the dysregulation of glucose transported vesicle 4 (GLUT4) translocation. Our study provides a possible evidence on structural basis dysfunction of AS160 and how the association of phosphorylated AS160 with 14-3-3, a downstream binding partner regulating GLUT4 translocation got disrupted due to T642A mutation. We initially derived the disease-causing mutation (Thr642Ala) among others through in-silico based statistical analysis. Subsequently, we interpreted the perturbation induced in the structural arrangement and their impaired interaction in core regions. Due to mutation, the key interfacial interactions between AS160-14-3-3 were changing from Thr642-Asp756, Thr642-Asp757, and Thr642-Lys659 for phosphorylated form, to Ala642-Val681 for mutant. Further, for phosphorylated AS160 the hotspot residues observed were Glu-629, Gln-635, His-641, Lys-653 and Arg-842 which changed to Arg-637, His-641 for mutant. Eventually, the molecular dynamics analysis revealed that local region for phosphorylation site of AS160 is reducing the flexibility, whereas mutation is making the region more flexible. Principal component analysis and Free energy landscape analysis together reveals phosphorylated AS160 is occupying less phase space with more stable landscape when compared to mutant. Our study strongly confers the destabilizing effect of the point mutation at a conserved site providing novel insights right down to the residual level of the conformational dysregulation of AS160 implicating deadly disease.

Investigation of conformational dynamics of Tyr89Cys mutation in protection of telomeres 1 gene associated with familial melanoma

Protection of telomeres 1 (POT1) is a component of the shelterin complex which is crucial for the regulation of telomere length and maintenance. Many naturally occurring mutations in the POT1 gene have been found to be associated with cardiac angiosarcoma, glioma, familial melanoma, and chronic lymphocytic leukemia. In particular, Y89C is a naturally occurring mutation of POT1, responsible for familial melanoma, and the molecular basis of this mutation is unexplored. In this study, we have extensively analyzed the structure of WT and Y89C mutant of POT1 to see the change in the conformational dynamics, free energy landscape, molecular motions and configurational frustration using molecular dynamics (MD) and other bioinformatics approaches. Y89C mutation shows a significant change in the backbone orientation, compactness, residual fluctuation, solvent accessibility, and hydrogen bonding, suggesting an overall destabilization of the protein structure. Besides, essential dynamics, conformation, magnitude, direction of motion and frustration analysis further suggesting the structural loss in POT1 due to Y89C mutation. Free energy landscape analysis also indicates the presence of a single well-defined free-energy minima in case of WT compared to multiple wells defined free energy minima observed in Y89C, clearly suggesting that this mutation leads to reduce the stability of POT1. This study possibly provides a valuable path to understand the molecular basis of Y89C-mediated development of familial melanoma.Communicated by Ramaswamy H. Sarma

Delineating the conformational dynamics of intermediate structures on the unfolding pathway of β-lactoglobulin in aqueous urea and dimethyl sulfoxide

The funnel shaped energy landscape model of the protein folding suggests that progression of folding proceeds through multiple pathways, having the multiple intermediates which leads to multidimensional free-energy surface. Herein, we applied all-atom MD simulation to conduct a comparative study on the structure of β-lactoglobulin (β-LgA) in aqueous mixture of 8 M urea and 8 M dimethyl sulfoxide (DMSO), at different temperatures. The cumulative results of multiple simulations suggest a common unfolding pathway of β-LgA, occurred through the stable and meta-stable intermediates (I), in both urea and DMSO. However, the free-energy landscape (FEL) analyses show that the structural transitions of I-states are energetically different. In urea, FEL shows distinct ensemble of intermediates, I1 and I2, separated by the energy barrier of ∼3.0 kcal mol−1. Similarly, we find the population of two distinct I1 and I2 states in DMSO, however, the I1 appeared transiently around ∼30–35 ns and is short-lived. But, the I2 ensemble is observed structurally compact and long-lived (∼50–150 ns) as compared to unfolding in urea. Furthermore, the I1 and I2 are separated through a high energy barrier of ∼6.0 kcal mol−1. Thus, our results provide the structural insights of intermediates which essentially bear the signature of a different unfolding pathway of β-LgA in urea and DMSO.Abbreviationsβ-LgAβ-lactoglobulinDMSOdimethyl sulfoxideFELfree-energy landscapeGdmClguanidinium chlorideIintermediate stateMGmolten globule statePMEparticle mesh EwaldQfraction of native contactsRMSDroot mean square deviationRMSFroot mean square fluctuationRgradius of gyrationSASAsolvent Accessible Surface AreascSASAthe side chain SASATrptryptophanCommunicated by Ramaswamy H. Sarma

Probing the Conformational and Energy Landscapes of KRAS Membrane Orientation.

Membrane reorientation of oncogenic RAS proteins is emerging as an important modulator of their functions. Previous studies have shown that the most common orientations include those with either the three C-terminal α-helices (OS1) or N-terminal β-strands (OS2) of the catalytic domain facing the membrane. OS1 and OS2 differ by the degree to which the effector-interacting surface is occluded by the membrane. However, the relative stability of these states and the rates of transition between them remained undetermined. How mutations might modulate preferences for specific orientation states is also far from clear. The current work attempted to address these questions through a comprehensive analysis of two 20 μs-long atomistic molecular dynamics simulations. The simulations were conducted on the oncogenic G12D and Q61H KRAS mutants bound to an anionic lipid bilayer. G12D and Q61H are among the most prevalent cancer-causing mutations at the P-loop and switch 2 regions of KRAS, respectively. We found that both mutants fluctuate in a similar manner between OS1 and OS2 via an intermediate orientation OS0, and both favor the signaling competent OS1 and OS0 over the occluded OS2. However, they differ in the details, such as in the extent to which they sample OS1. Analysis of the orientation free-energy landscapes estimated from the simulations indicate that OS1 and OS2 are the most stable states. However, the overall free energy surface is rugged, indicating a large diversity of conformations including at least two substates in each orientation state that differ in stability only by about 0.5-1.0 kcal/mol. Reversible transitions between OS1 and OS2 occur via two well-defined pathways that traverse OS0. In the minimum energy path, helix 4 remains close to the membrane as the angle of the catalytic domain from the membrane plane changes, resulting in a barrier of ∼1 kcal/mol for OS1/OS2 interconversions. Estimation of the rates of the various transitions based on survival probabilities yielded two rate constants in the order of 107 and 106 s-1, which we attribute to intrinsic protein conformational dynamics and transient protein-lipid interactions, respectively. The faster process dominates every transition, confirming a previous suggestion that RAS membrane reorientation is driven by conformational fluctuations rather than protein-lipid interactions.

The investigation of the G-quadruplex aptamer selectivity to Pb2+ ion: a joint molecular dynamics simulation and density functional theory study

The aptamers with the ability to form a G-quadruplex structure can be stable in the presence of some ions. Hence, study of the interactions between such aptamers and ions can be beneficial to determine the highest selective aptamer toward an ion. In this article, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations have been applied to investigate the selectivity of the T30695 aptamer toward Pb2+ in comparison with some ions. The Free Energy Landscape (FEL) analysis indicates that Pb2+ has remained inside the aptamer during the MD simulation, while the other ions have left it. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding energies prove that the conformational stability of the aptamer is the highest in the presence of Pb2+. According to the compaction parameters, the greatest compressed ion-aptamer complex, and hence, the highest ion-aptamer interaction have been induced in the presence of Pb2+. The contact maps clarify the closer contacts between the nucleotides of the aptamer in the presence of Pb2+. The density functional theory (DFT) results show that Pb2+ forms the most stable complex with the aptamer, which is consistent with the MD results. The QM calculations reveal that the N-H bonds and the O…H distances are the longest and the shortest, respectively, in the presence of Pb2+. The obtained results verify that the strongest hydrogen bonds (HBs), and hence, the most compressed aptamer structure are induced by Pb2+. Besides, atoms in molecules (AIM) and natural bond orbital (NBO) analyses confirm the results.Communicated by Ramaswamy H. Sarma

Structural Analysis and Conformational Dynamics of STN1 Gene Mutations Involved in Coat Plus Syndrome

The human CST complex (CTC1–STN1–TEN1) is associated with telomere functions including genome stability. We have systemically analyzed the sequence of STN and performed structure analysis to establish its association with the Coat Plus (CP) syndrome. Many deleterious non-synonymous SNPs have been identified and subjected for structure analysis to find their pathogenic association and aggregation propensity. A 100-ns all-atom molecular dynamics simulation of WT, R135T, and D157Y structures revealed significant conformational changes in the case of mutants. Changes in hydrogen bonds, secondary structure, and principal component analysis further support the structural basis of STN1 dysfunction in such mutations. Free energy landscape analysis revealed the presence of multiple energy minima, suggesting that R135T and D157Y mutations destabilize and alter the conformational dynamics of STN1 and thus may be associated with the CP syndrome. Our study provides a valuable direction to understand the molecular basis of CP syndrome and offer a newer therapeutics approach to address CP syndrome.

Identification of AtHsp90.6 involved in early embryogenesis and its structure prediction by molecular dynamics simulations.

Heat-shock protein of 90 kDa (Hsp90) is a key molecular chaperone involved in folding the synthesized protein and controlling protein quality. Conformational dynamics coupled to ATPase activity in N-terminal domain is essential for Hsp90's function. However, the relevant process is still largely unknown in plant Hsp90s, especially those required for plant embryogenesis which is inextricably tied up with human survival. Here, AtHsp90.6, a member of Hsp90 family in Arabidopsis, was firstly identified as a protein essential for embryogenesis. Thus we modelled AtHsp90.6 in its functionally closed ‘lid-down’ and open ‘lid-up’ states, exploring the nucleotide binding mechanism in these two states. Free energy landscape and electrostatic potential analysis revealed the switching mechanism between these two states. Collectively, this study quantitatively analysed the conformational changes of AtHsp90.6 bound to ATP or ADP. This result may help us understand the mechanism of action of AtHsp90.6 in future.

Molecular mechanism of acetoacetyl-CoA enhanced kinetics for increased bioplastic production from Cupriavidus necator 428

Polyhydroxyalkanoates are gaining importance due to their biodegradable nature and close analogy to plastics. Polyhydroxybutyrate (PHB) is the most widely used bioplastic from polyalkanoate family, which is produced by a legion of bacterial species via phbCAB operon encoding β-ketothiolase (PhaA), NADPH-dependent acetoacetyl-coenzyme A (acetoacetyl-CoA) reductase (PhaB) and polyhydroxyalkanoate synthase (PhaC). Augmentation in the activity of these enzymes is promising for increased PHB production which is achieved by enzyme engineering strategies including non-structural and structural approaches. Our study is deployed on directed evolution-based experimentally reported mutants of PhaB enzyme with increased efficiency due to impact on critical structural factors. We have analyzed and compared the native PhaB with two of its variants Q47L and T173S in complex with their cofactor i.e. NADPH as well as the substrate i.e. acetoacetyl-CoA, via long range molecular dynamics simulations. Interaction profile, MMPBSA, essential dynamics, and free energy landscape analysis revealed that the enzyme efficiency is critically affected by cofactor interactions. It was also observed that mutants have higher equilibrium constant with lesser but optimal affinity for substrate and cofactor than the wild type, which might be the reason for increased efficiency of the mutants via enhanced substrate and cofactor exchange rate. Our study provides insights into the cofactor and substrate binding affinities to PhaB enzyme at atomistic level, which will facilitate designing of highly efficient PhaB enzymes for increased PHB production. Communicated by Ramaswamy H. Sarma

3D-QSAR and molecular recognition of Klebsiella pneumoniae NDM-1 inhibitors

ABSTRACTNew Delhi metallo-β-lactamase-1 (NDM-1) as a target for the development of anti-superbug agents, plays an important role in the resistance of β-lactam antibiotics and has received worldwide attention. Sulfhydryl propionic acid derivatives can effectively inhibit the catalytic activity of NDM-1, but the quantitative structure–activity relationship (QSAR) and inhibitor-target recognition mechanism both remain unclear. In this work, CoMFA and CoMSIA models of sulfhydryl propionic acid inhibitors with high predictive ability were obtained, from which the effect of different substituents on the inhibitory activity against NDM-1 were revealed at the molecular level. Then, two 120-nanosecond comparative molecular dynamics (MD) simulations for NDM-1 enzyme and NDM-1-inhibitor complex systems were performed to study the recognition and inhibition mechanism of sulfhydryl propionic acid derivatives. The binding of inhibitors alters the entire H-bond network of the NDM-1 system accompanied by the formation of strong interactions with I35, W93, H120, H122, D124, H189 and H250, further weakens the recognition of NDM-1 by its endogenic substrates. Finally, the results of free energy landscape and conformation cluster analyses show that NDM-1 underwent a significant conformational change after the association with sulfhydryl propionic acid inhibitors. Our findings can provide theoretical support and help for anti-superbug agents design based on the structures of NDM-1 and sulfhydryl propionic acid derivatives.

Insights into the substrate binding specificity of quorum-quenching acylase PvdQ

Quorum sensing is a cell to cell signaling mechanism that enables them to coordinate their behaviors in a density-dependent manner mediated by small diffusible signaling molecules, which can control the virulence and biofilm gene expression in many Gram-negative and positive bacteria. N-acyl homoserine lactone acylase PvdQ from human opportunistic pathogen Pseudomonas aeruginosa is a quorum-quenching enzyme that can hydrolyze the amide bond of the quorum signaling N-acyl homoserine lactones (AHLs) thereby degrading the signaling molecules, turning off the biofilm phenotype and resulting in a reduction of bacterial virulence. Previous studies demonstrated that PvdQ has different preferences for N-acyl substrates with different acyl chain lengths and substituents. However, the substrate binding specificity determinants of the quorum-quenching enzyme PvdQ with the different bacterial ligands are unknown and unintuitive. Further, elucidation of these determinants can lead to mutants with efficiency and broader substrate promiscuity. To investigate this question, a computational study was carried out combining multiple molecular docking methods, molecular dynamics simulations, residue interaction network analysis, and binding free energy calculations. The main findings are: firstly, the results from pKa predictions support that the pKa of the N-terminus of Serβ1 was depressed due to the surrounding residues. Multiple molecular docking studies provide useful information about the detailed binding modes and binding affinities. Secondly, 300 ns molecular dynamics simulations were carried out to analyze the overall molecular motions of substrate-bound and substrate-free PvdQ. The specific interactions between the active site of PvdQ and different ligands revealed the determinants for the preference among the ligands. A systematic comparison and analysis of the protein dynamic fingerprint of each complex demonstrated that binding of the most favorable ligand, C12-homoserine lactone (C12-HSL), reduced the global motions of the complex and maintained the correct arrangement of the catalytic site. Further, the residue interaction network analysis of each system illustrated that there are more communication contacts and pathways between the residues in the C12-HSL complex as compared to complexes with the other ligands. The binding of the C12-HSL ligand facilitates structural communication between the two knobs and the active site. While the binding of the other ligands tend to impair specific communication pathways between the two knobs and the active site, and lead to a catalytically inefficient state. Finally, simulation results from free energy landscape and binding free energy analysis revealed that the C12-HSL ligand has the lowest binding free energy and greater stability than the less favored ligands. Each of the following residues: Serβ1, Hisβ23, Pheβ24, Metβ30, Pheβ32, Leuβ50, Asnβ57, Thrβ69, Valβ70, Trpβ162, Trpβ186, Asnβ269, Argβ297 and Leuα146, play different roles in substrate binding specificity. This is the first computational study that provides molecular information for structure-dynamic-function relationships of PvdQ with different ligands and demonstrates determinants of bacterial substrate binding specificity.

Conformational modifications induced by internal tandem duplications on the FLT3 kinase and juxtamembrane domains.

The aberrant expression of FLT3 tyrosine kinase is associated primarily with acute myeloid leukaemia. This blood malignancy is often related to the onset of internal tandem duplications (ITDs) in the native sequence of the protein. The ITDs occur mainly in the juxtamembrane domain of the protein and alter the normal activity of the enzyme. In this work, we have studied the native form of FLT3 and six mutants by molecular dynamics simulations. The catalytic activity of FLT3 is exerted by the tyrosine kinase domain (KD) and regulated by the juxtamembrane (JM) domain. Analysis of the dynamics of these two domains have shown that the introduction of ITDs in the JM domain alters both structural and dynamic parameters. The presence of ITDs allowed the protein to span a larger portion of the conformational space, particularly in the JM domain and the activation loop. The FLT3 mutants were found to adopt more stable configurations than the native enzyme. This was due to the different arrangements assumed by the JM domain. Larger fluctuations of the activation loop were found in four of the six mutants. In the native FLT3, the key residue Tyr572 is involved in a strong and stable interaction with an ion pair. This interaction, which is thought to keep the JM in place hence regulating the activity of the enzyme, was found to break in all FLT3 mutants.

Insights into the DNA binding induced thermal stabilization of transcription factor FOXP3

The transcription factor FOXP3 is required for the development and function of regulatory T cells. Here, we studied the dynamics of FOXP3 in the presence and absence of DNA target sequence. Multiple molecular dynamics (MD) simulations were employed to investigate the role of DNA in enhancing the stability of FOXP3 via protein–DNA interactions, and to study the structural transition in protein at three different temperatures (300, 350, and 400 K). Results indicate that FOXP3 is stabilized by DNA even though the temperature rises up to 400 K. FOXP3 is found to undergo significant conformational change at a higher temperature, however, DNA provides greater rigidity and lower overall conformational flexibility, resulting in the global stabilization of FOXP3. The conformational restriction of FOXP3-DNA complex is probed with essential dynamics (ED). Free-energy landscape analysis at high temperature shows the presence of metastable intermediates, suggesting the reason behind the observed thermal stability. Secondary structural snapshots clearly indicate the presence of non-native interactions between DNA and protein. This study increases our understanding of the dynamic behaviors and the interaction mechanism of FOXP proteins and DNA at the atomic level and offers a model for studying the structural biology of transcriptional regulation.Communicated by Ramaswamy H. Sarma

Pyrazinamide drug resistance in RpsA mutant (∆438A) of Mycobacterium tuberculosis: Dynamics of essential motions and free-energy landscape analysis.

Pyrazinamide is an essential first-line antitubercular drug which plays pivotal role in tuberculosis treatment. It is a prodrug that requires amide hydrolysis by mycobacterial pyrazinamidase enzyme for conversion into pyrazinoic acid (POA). POA is known to target ribosomal protein S1 (RpsA), aspartate decarboxylase (PanD), and some other mycobacterial proteins. Spontaneous chromosomal mutations in RpsA have been reported for phenotypic resistance against pyrazinamide. We have constructed and validated 3D models of the native and Δ438A mutant form of RpsA protein. RpsA protein variants were then docked to POA and long range molecular dynamics simulations were carried out. Per residue binding free-energy calculations, free-energy landscape analysis, and essential dynamics analysis were performed to outline the mechanism underlying the high-level PZA resistance conferred by the most frequently occurring deletion mutant of RpsA. Our study revealed the conformational modulation of POA binding site due to the disruptive collective modes of motions and increased conformational flexibility in the mutant than the native form. Residue wise MMPBSA decomposition and protein-drug interaction pattern revealed the difference of energetically favorable binding site in the wild-type (WT) protein in comparison with the mutant. Analysis of size and shape of minimal energy landscape area delineated higher stability of the WT complex than the mutant form. Our study provides mechanistic insights into pyrazinamide resistance in Δ438A RpsA mutant, and the results arising out of this study will pave way for design of novel and effective inhibitors targeting the resistant strains of Mycobacterium tuberculosis.

Binding of Cytotoxic Aβ25-35 Peptide to the Dimyristoylphosphatidylcholine Lipid Bilayer.

Aβ25-35 is a short, cytotoxic, and naturally occurring fragment of the Alzheimer's Aβ peptide. To map the molecular mechanism of Aβ25-35 binding to the zwitterionic dimyristoylphosphatidylcholine (DMPC) bilayer, we have performed replica exchange with solute tempering molecular dynamics simulations using all-atom explicit membrane and water models. Consequences of sequence truncation on the binding mechanism have been measured by utilizing as a control our previous simulations probing binding of the longer peptide Aβ10-40 to the same bilayer. The most intriguing feature of Aβ25-35 binding to the DMPC bilayer is a coexistence of two bound states with strikingly different characteristics: a dominant surface-bound state and a less stable inserted state. In the surface-bound state, the peptide samples extended conformations, in which its unbound C-terminal is pointed away from the bilayer. In contrast, in the inserted state, the C-terminal resides deep in the bilayer hydrophobic core. In both states, the N-terminal remains anchored to the bilayer. Free energy landscape analysis reveals that the two states are separated by a moderate barrier, suggesting that Aβ25-35 monomer may frequently interconvert between them. The net effect of Aβ25-35 binding is a minor impact on the bilayer structure, which contrasts with the considerable bilayer perturbations induced by a longer Aβ10-40 peptide penetrating deep into the bilayer core. Therefore, we conclude that the binding mechanisms of Aβ25-35 and Aβ10-40 peptides are different. Potential implications of our results for Aβ25-35 cytotoxicity are discussed. A comparison of experimental data with our findings reveals a good agreement.

Insights into the structure–function relationship of brown plant hopper resistance protein, Bph14 of rice plant: a computational structural biology approach

Brown plant hopper (BPH) is one of the major destructive insect pests of rice, causing severe yield loss. Thirty-two BPH resistance genes have been identified in cultivated and wild species of rice Although, molecular mechanism of rice plant resistance against BPH studied through map-based cloning, due to non-existence of NMR/crystal structures of Bph14 protein, recognition of leucine-rich repeat (LRR) domain and its interaction with different ligands are poorly understood. Thus, in the present study, in silico approach was adopted to predict three-dimensional structure of LRR domain of Bph14 using comparative modelling approach followed by interaction study with jasmonic and salicylic acids. LRR domain along with LRR-jasmonic and salicylic acid complexes were subjected to dynamic simulation using GROMACS, individually, for energy minimisation and refinement of the structure. Final binding energy of jasmonic and salicylic acid with LRR domain was calculated using MM/PBSA. Free-energy landscape analysis revealed that overall stability of LRR domain of Bph14 is not much affected after forming complex with jasmonic and salicylic acid. MM/PBSA analysis revealed that binding affinities of LRR domain towards salicylic acid is higher as compared to jasmonic acid. Interaction study of LRR domain with salicylic acid and jasmonic acid reveals that THR987 of LRR form hydrogen bond with both complexes. Thus, THR987 plays active role in the Bph14 and phytochemical interaction for inducing resistance in rice plant against BPH. In future, Bph14 gene and phytochemicals could be used in BPH management and development of novel resistant varieties for increasing rice yield.

Structural, molecular motions, and free-energy landscape of Leishmania sterol-14α-demethylase wild type and drug resistant mutant: a comparative molecular dynamics study

Sterol-14α-demethylase (CYP51) is an ergosterol pathway enzyme crucial for the survival of infectious Leishmania parasite. Recent high-throughput metabolomics and whole genome sequencing study revealed amphotericin B resistance in Leishmania is indeed due to mutation in CYP51. The residue of mutation (asparagine 176) is conserved across the kinetoplastidae and not in yeast or humans, portraying its functional significance. In order to understand the possible cause for the resistance, knowledge of structural changes due to mutation is of high importance. To shed light on the structural changes of wild and mutant CYP51, we conducted comparative molecular dynamics simulation study. The active site, substrate biding cavity, substrate channel entrance (SCE), and cavity involving the mutated site were studied based on basic parameters and large concerted molecular motions derived from essential dynamics analyses of 100 ns simulation. Results indicated that mutant CYP51 is stable and less compact than the wild type. Correspondingly, the solvent accessible surface area (SASA) of the mutant was found to be increased, especially in active site and cavities not involving the mutation site. Free-energy landscape analysis disclosed mutant to have a rich conformational diversity than wild type, with various free-energy conformations of mutant having SASA greater than wild type with SCE open. More residues were found to interact with the mutant CYP51 upon docking of substrate to both the wild and mutant CYP51. These results indicate that, relative to wild type, the N176I mutation of CYP51 in Leishmania mexicana could possibly favor increased substrate binding efficiency.