Project description:Prokaryotic CRISPR-Cas adaptive immune systems utilize sequence-specific RNA-guided endonucleases to defend against infection by viruses, bacteriophages, and mobile elements, while these foreign genetic elements evolve diverse anti-CRISPR proteins to overcome the CRISPR-Cas-mediated defense of the host. Recently, AcrIIA2 and AcrIIA4, encoded by Listeria monocytogene prophages, were shown to block the endonuclease activity of type II-A Streptococcus pyogene Cas9 (SpyCas9). We now report the crystal structure of AcrIIA4 in complex with single-guide RNA-bound SpyCas9, thereby establishing that AcrIIA4 preferentially targets critical residues essential for PAM duplex recognition, as well as blocks target DNA access to key catalytic residues lining the RuvC pocket. These structural insights, validated by biochemical assays on key mutants, demonstrate that AcrIIA4 competitively occupies both PAM-interacting and non-target DNA strand cleavage catalytic pockets. Our studies provide insights into anti-CRISPR-mediated suppression mechanisms for inactivating SpyCas9, thereby broadening the applicability of CRISPR-Cas regulatory tools for genome editing.
Project description:Effective data reduction methods are necessary for uncovering the inherent conformational relationships present in large molecular dynamics (MD) trajectories. Clustering algorithms provide a means to interpret the conformational sampling of molecules during simulation by grouping trajectory snapshots into a few subgroups, or clusters, but the relationships between the individual clusters may not be readily understood. Here we show that network analysis can be used to visualize the dominant conformational states explored during simulation as well as the connectivity between them, providing a more coherent description of conformational space than traditional clustering techniques alone. We compare the results of network visualization against 11 clustering algorithms and principal component conformer plots. Several MD simulations of proteins undergoing different conformational changes demonstrate the effectiveness of networks in reaching functional conclusions.
Project description:The temperature dependence of helical propensities for the peptides Ac-ZGG-(KAAAA)(3)X-NH(2) (Z = Y or G, X = A, K, and D-Arg) were studied both experimentally and by MD simulations. Good agreement is observed in both the absolute helical propensities as well as relative helical content along the sequence; the global minimum on the calculated free energy landscape corresponds to a single alpha-helical conformation running from K4 to A18 with some terminal fraying, particularly at the C-terminus. Energy component analysis shows that the single helix state has favorable intramolecular electrostatic energy due to hydrogen bonds, and that less-favorable two-helix globular states have favorable solvation energy. The central lysine residues do not appear to increase helicity; however, both experimental and simulation studies show increasing helicity in the series X = Ala --> Lys --> D-Arg. This C-capping preference was also experimentally confirmed in Ac-(KAAAA)(3)X-GY-NH(2) and (KAAAA)(3)X-GY-NH(2) sequences. The roles of the C-capping groups, and of lysines throughout the sequence, in the MD-derived ensembles are analyzed in detail.
Project description:Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease caused by the misfolding of Cu, Zn superoxide dismutase (SOD1). Several earlier studies have shown that monomeric apo SOD1 undergoes significant local unfolding dynamics and is the predecessor for aggregation. Here, we have employed atomistic molecular dynamics (MD) simulations to study the structure and dynamics of monomeric apo and holo SOD1 in water, aqueous urea and aqueous urea-TMAO (trimethylamine oxide) solutions. Loop IV (zinc-binding loop) and loop VII (electrostatic loop) of holo SOD1 are considered as functionally important loops as they are responsible for the structural stability of holo SOD1. We found larger local unfolding of loop IV and VII of apo SOD1 as compared to holo SOD1 in water. Urea induced more unfolding in holo SOD1 than apo SOD1, whereas the stabilization of both the form of SOD1 was observed in ternary solution (i.e. water/urea/TMAO solution) but the extent of stabilization was higher in holo SOD1 than apo SOD1. The partially unfolded structures of apo SOD1 in water, urea and holo SOD1 in urea were identified by the exposure of the hydrophobic cores, which are highly dynamic and these may be the initial events of aggregation in SOD1. Our simulation studies support the formation of aggregates by means of the local unfolding of monomeric apo SOD1 as compared to holo SOD1 in water.
Project description: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.
Project description:Interactions of amyloid-β (Aβ) with neuronal membrane are associated with the progression of Alzheimer's disease (AD). Ganglioside GM1 has been shown to promote the structural conversion of Aβ and increase the rate of peptide aggregation; but the exact nature of interaction driving theses processes remains to be explored. In this work, we have carried out atomistic-scale computer simulations (totaling 2.65 µs) to investigate the behavior of Aβ monomer and dimers in GM1-containing raft-like membrane. The oligosaccharide head-group of GM1 was observed to act as scaffold for Aβ-binding through sugar-specific interactions. Starting from the initial helical peptide conformation, a β-hairpin motif was formed at the C-terminus of the GM1-bound Aβ-monomer; that didn't appear in absence of GM1 (both in fluid POPC and liquid-ordered cholesterol/POPC bilayers and also in aqueous medium) within the simulation time span. For Aβ-dimers, the β-structure was further enhanced by peptide-peptide interactions, which might influence the propensity of Aβ to aggregate into higher-ordered structures. The salt-bridges and inter-peptide hydrogen bonds were found to account for dimer stability. We observed spontaneous formation of intra-peptide D(23)-K(28) salt-bridge and a turn at V(24)GSN(27) region - long been accepted as characteristic structural-motifs for amyloid self-assembly. Altogether, our results provide atomistic details of Aβ-GM1 and Aβ-Aβ interactions and demonstrate their importance in the early-stages of GM1-mediated Aβ-oligomerisation on membrane surface.
Project description:Interleukin-2-inducible T-cell kinase (ITK) regulates the response to T-cell receptor signaling and is a drug target for inflammatory and immunological diseases. Molecules that bind preferentially to the active form of ITK have low selectivity between kinases, whereas those that bind preferentially to the inactive form have high selectivity for ITK. Therefore, computational methods to predict the conformational selectivity of compounds are required to design highly selective ITK inhibitors. In this study, we performed absolute binding free-energy perturbation (ABFEP) simulations for 11 compounds on both active and inactive forms of ITK, and the calculated binding free energies were compared with experimental data. The conformational selectivity of 10 of the 11 compounds was correctly predicted using ABFEP. To investigate the mechanism underlying the stabilization of the active and inactive structures by the compounds, we performed extensive, conventional molecular dynamics simulations, which revealed that the compound-induced stabilization of the P-loop and linkage of conformational changes in L489, V419, F501, and M410 upon compound binding were critical factors. A guideline for designing inactive-form binders is proposed based on these key structural factors. The ABFEP and the created guidelines are expected to facilitate the discovery of highly selective ITK inhibitors.
Project description:Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.
Project description:The conformational ensemble of intrinsically disordered proteins, such as α-synuclein, are responsible for their function and malfunction. Misfolding of α-synuclein can lead to neurodegenerative diseases, and the ability to study their conformations and those of other intrinsically disordered proteins under varying physiological conditions can be crucial to understanding and preventing pathologies. In contrast to well-folded peptides, a consensus feature of IDPs is their low hydropathy and high charge, which makes their conformations sensitive to pH perturbation. We examine a prominent member of this subset of IDPs, α-synuclein, using a divide-and-conquer scheme that provides enhanced sampling of IDP structural ensembles. We constructed conformational ensembles of α-synuclein under neutral (pH ~ 7) and low (pH ~ 3) pH conditions and compared our results with available information obtained from smFRET, SAXS, and NMR studies. Specifically, α-synuclein has been found to in a more compact state at low pH conditions and the structural changes observed are consistent with those from experiments. We also characterize the conformational and dynamic differences between these ensembles and discussed the implication on promoting pathogenic fibril formation. We find that under low pH conditions, neutralization of negatively charged residues leads to compaction of the C-terminal portion of α-synuclein while internal reorganization allows α-synuclein to maintain its overall end-to-end distance. We also observe different levels of intra-protein interaction between three regions of α-synuclein at varying pH and a shift towards more hydrophilic interactions with decreasing pH.