Why Does the E1219V Mutation Expand T-Rich PAM Recognition in Cas9 from Streptococcus pyogenes?

Popular RNA-guided DNA endonuclease Cas9 from Streptococcus pyogenes (SpCas9) recognizes the canonical 5'-NGG-3' protospacer adjacent motif (PAM) and triggers double-stranded DNA cleavage activity. Mutations in SpCas9 were demonstrated to expand the PAM readability and hold promise for therapeutic and genome editing applications. However, the energetics of the PAM recognition and its relation to the atomic structure remain unknown. Using the X-ray structure (precatalytic SpCas9:sgRNA:dsDNA) as a template, we calculated the change in the PAM binding affinity in response to SpCas9 mutations using computer simulations. The E1219V mutation in SpCas9 fine-tunes the water accessibility in the PAM binding pocket and promotes new interactions in the SpCas9:noncanonical T-rich PAM, thus weakening the PAM stringency. The nucleotide-specific interaction of two arginine residues (i.e., R1333 and R1335 of SpCas9) ensured stringent 5'-NGG-3' PAM recognition. R1335A substitution (SpCas9R1335A) completely disrupts the direct interaction between SpCas9 and PAM sequences (canonical or noncanonical), accounting for the loss of editing activity. Interestingly, the double mutant (SpCas9R1335A,E1219V) boosts DNA binding affinity by favoring protein:PAM electrostatic contact in a desolvated pocket. The underlying thermodynamics explain the varied DNA cleavage activity of SpCas9 variants. A direct link between the energetics, structures, and activity is highlighted, which can aid in the rational design of improved SpCas9-based genome editing tools.

Effect of Spacer Length Modification of the Cationic Side Chain on the Energetics of Antimicrobial Peptide Binding to Membrane-Mimetic Bilayers.

Understanding the mechanism of action of the antimicrobial peptide (AMP) in terms of its structure and energetics is the key to designing new potent and selective AMPs. Recently, we reported a membranolytic 14-residue-long lysine-rich cationic antimicrobial peptide (LL-14: NH3+-LKWLKKLLKWLKKL-CONH2) against Pseudomonas aeruginosa, Klebsiella pneumoniae, and Staphylococcus aureus, which is limited by cytotoxicity and expected to undergo facile protease degradation. Aliphatic side-chain-length modification of the cationic amino-acid residues (Lys and Arg) is a popular strategy for designing protease-resistant AMPs. However, the effect of the peptide side-chain length modifications on the membrane binding affinity and its relation to the atomic structure remain an unsolved problem. We report computer simulations that quantitatively calculated the difference in peptide binding affinity to membrane-mimetic-bilayer models (bacterial: 1-palmitoyl-2-oleoyl-phosphatidylethanolamine (POPE)/1-palmitoyl-2-oleoyl-phosphatidylglycerol (POPG) bilayer and mammalian: 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer) upon decreasing or increasing the spacer length of the cationic lysine residues of LL-14 (as well as their arginine analogues). We show that the peptide/bilayer interaction energetics varies drastically in response to spacer length modification. The strength of peptide discrimination depends strongly on the nature of the bilayer (bacterial or mammalian mimetic model). An increase in the lysine spacer length by one carbon (i.e., homolysine analogue of LL-14) is weakly/strongly disfavored by the bacterial/mammalian-membrane-mimetic bilayer. Recently, we have demonstrated an excellent correlation between the antimicrobial activity of the membranolytic cationic peptides and their binding affinity to membrane-mimetic-bilayer models. Thus, the homolysine analogue of LL-14 is a promising noncytotoxic AMP with conserved activity. On the other hand, homoarginine analogue (arginine spacer length increment by a single carbon) was preferred by both the bacteria and the mammalian mimetic bilayers and displayed the strongest affinity for the former among the peptides studied in this work. Thus, the promising most potent homoarginine analogue is likely to be cytotoxic. Shortening the Lys/Arg side chain to a three-carbon spacer (Dab/Agb) improves the binding affinity to bacterial and mammalian-membrane-mimetic bilayers. Arginine and arginine-derivative peptides exhibited stronger binding affinity to the bilayers relative to the lysine analogue. The results provide a plausible explanation to the previous experimental observations, viz., superior antimicrobial activity of the arginine peptides relative to Lys peptides and the improvement of antimicrobial activity upon substitution of Lys with Dab in the cationic peptides. The simulations revealed that the small change in the peptide hydrophobicity by Lys/Arg spacer length modification could drastically alter the energetics of peptide/bilayer binding by fine-tuning the electrostatic interactions. The energetics underlying the peptide selectivity by simple membrane-mimetic bilayer models may be beneficial for designing new selective and protease-resistant AMPs.

Bridging Thermodynamics, Antimicrobial Activity, and pH Sensitivity of Cationic Membranolytic Heptapeptides-A Computational and Experimental Study.

Understanding the thermodynamics of peptide:membrane binding and the factors that alter the stability is the key to designing potent and selective small antimicrobial peptides. Here, we report the thermodynamics, antimicrobial activity, and mechanism of a de novo designed seven residue long cationic antimicrobial peptide (P4: NH3+-LKWLKKL-CONH2, Charge +4) and its analogs (P5: Lysine's → Arginine's; P6: Lysine's → Uncharged-Histidine's; P7: Tryptophan → Leucine) by combining computation and experiments. Computer simulations predicted the order of decreasing peptide binding affinity to the membrane-mimetic systems (micelle/bilayer) as P5 > P4 > P7 ≫ P6. Antimicrobial assays of these peptides against P. aeruginosa and E. coli at physiological pH 7.4 confirmed P5 as the most potent peptide (followed by P4), whereas P6 showed inferior activity. P7 was found to be inactive against E. coli. Substitution of the uncharged-histidine (P6) by the charged-histidine (P6*) significantly favored micelle/bilayer binding. Thus, P6 was predicted to be an effective antimicrobial peptide only at low pH. Noticeable improvement in the antimicrobial activity of the histidine-peptide (P6) against E. coli (an acid-resistant bacteria) upon lowering the pH was demonstrated and validated the computational claim. The peptides displayed a membranolytic mode of action. The link between the structure and calculated energetics (ΔΔG) has been established, and the correlation between the calculated energetics and the antimicrobial activity has been highlighted. The histidine-peptide (P6) is reported to be active against acid-resistant bacteria, thus, a promising membranolytic pH-sensitive AMP.

Structure-based thermodynamics of ion selectivity (Mg2+versus Ca2+ and K+versus Na+) in the active site of the eukaryotic lariat group II intron from algae Pylaiella littoralis.

Group II introns are metalloenzymes that can catalyze self-splicing. Recently, the crystal structures of the eukaryotic group IIB lariat intron from the brown algae Pylaiella littoralis have been reported for two intermediate states (pre-hydrolytic (2s) and post-hydrolytic) along the self-splicing pathway. Three characteristic metal-ion binding sites (M1 and M2 sites for catalytic Mg2+ ions, and K1 site for K+) in the catalytic pocket of the lariat intron have been identified and proposed to be crucial for self-splicing. Using the X-ray structures as a template, we quantitatively estimated the energetics of divalent (Mg2+versus Ca2+) and monovalent (K+versus Na+) ion selectivity and established a direct link between the energetics and structures of this lariat intron (bound to cognate and near-cognate metal ions). Molecular dynamics (MD) free energy simulations showed that the lariat intron was strongly selective between divalent metal ions. The strength of divalent metal-ion selectivity was noticeably high in the post-hydrolytic state (ΔΔG ≈ 20 kcal mol-1) relative to its pre-hydrolytic (2s) state (ΔΔG ≈ 13 kcal mol-1). Quantum chemical calculations ensured that the sign of the estimated divalent metal-ion selectivity was correct. The M1-binding pocket was less solvent-exposed in the case of the post-hydrolytic state relative to the pre-hydrolytic (2s) state, which boosted the metal-ion selectivity of the former. Surprisingly, in contrast to the bacterial linear group II intron, the lariat intron was found to be non-selective between monovalent ions (K+versus Na+). The interaction network in the first coordination shell of Ca2+ in the M1-binding pocket was different relative to Mg2+. Mg2+ substitution by Ca2+ resulted in the substitution of a single M1-RNA interaction by the M1-water interaction. In the pre-hydrolytic (2s) state, Ca2+ substitution completely disrupted the M1⋯5'-exon interaction; thus, the nature of the divalent metal ion is critical for catalysis. The interaction network in the M2 site was independent of the nature of the divalent metal ions (Mg2+ or Ca2+). The monovalent ion was loosely bound in the wet binding pocket (K1 site) of the lariat intron; thus, the substitution of K+ by Na+ could not significantly alter the free energy of the complex. The metal ion selectivity was dependent on the solvent accessibility of the metal-ion-binding-pocket, dry pocket enhanced the selectivity.

Effect of a monovalent salt on the energetics of an antimicrobial-peptide: micelle dissociation.

Antimicrobial peptides (AMPs) are promising antimicrobial and therapeutic agents. Recently, we synthesized a cationic 14 residue AMP (LL-14: LKWLKKLLKWLKKL), which showed high broad-spectrum antimicrobial activity. However, the antimicrobial activity of LL-14 was compromised in the presence of NaCl. Salt sensitivity of antimicrobial potency is one of the fundamental limitations of AMP therapeutics. Thus, understanding the thermodynamics of AMP binding to simple membrane-mimetic systems and the effect of NaCl that contributes to their stability is crucial for designing promising AMPs against microbial infection. In this work, we reported computational analysis of LL-14 binding to SDS micelles (the simplest bacterial membrane mimic) at various NaCl concentrations (0.0%, 0.5%, 1.0% w/v). The thermodynamics of LL-14 dissociation from the SDS micelles was estimated by employing steered molecular dynamics (SMD) simulation followed by umbrella sampling. The results indicated that the increase in NaCl concentration systematically disfavoured the LL-14:SDS binding, primarily by stabilizing the dissociative state (i.e., free LL-14 and free micelles in water). We proposed a kinetic scheme in which the salt-induced selective stabilization of the dissociative state increased the activation barrier for the peptide:micelle binding event resulting in reduced affinity. Center-of-mass pulling indicated that the interactions involving the N-terminal of the LL-14 (residues 1-6) and SDS micelle were crucial for the stability of the LL-14:SDS complex, and LL-14 underwent a conformational change (helix → unstructured) before dissociating from the SDS micelle. The observed structural features from the peptide:micelle dissociation pathway corroborate our previous simulations as well as circular dichroism (CD), and fluorescence experiments.

Effect of Differential Geminal Substitution of γ Amino Acid Residues at the (i + 2) Position of αγ Turn Segments on the Conformation of Template ß-Hairpin Peptides.

The effect of insertion of three geminally dimethyl substituted γ amino acid residues [γ2,2 (4-amino-2,2-dimethylbutanoic acid), γ3,3 (4-amino-3,3-dimethylbutanoic acid), and γ4,4 (4-amino-4,4-dimethylbutanoic acid)] at the (i + 2) position of a two-residue αγ C12 turn segment in a model octapeptide sequence Leu-Phe-Val-Aib-Xxx-Leu-Phe-Val (where Xxx = γ amino acid residues) has been investigated in this study. Solution conformational studies (NMR, CD, and IR) and ab initio calculations indicated that γ3,3 and γ4,4 residues were well accommodated in the ß-hairpin nucleating αγ C12 turns, which gave rise to well-registered hairpins, in contrast to γ2,2, which was unable to form a tight C12 ß-hairpin nucleating turn and promote a well-registered ß-hairpin. Geminal disubstitution at the Cα carbon in γ2,2 led to unfavorable steric contacts, disabling its accommodation in the αγ C12 hairpin nucleating turn unlike the γ3,3 and γ4,4 residues. Geminal substitutions at different carbons along the backbone constrained backbone torsion angles for the three γ amino acid residues differently, generating diverse conformational preferences in them. Folded hairpins were energetically more stable (∼8 to 9 kcal/mol) than the unfolded peptides. Conformational preference of the peptides was independent of the N-terminal protecting group. Such fundamental understanding will instrumentalize the future directed design of foldamers.

Free energy calculations of RNA interactions.

This review discusses the use of molecular dynamics free energy calculations for characterizing RNA interactions, with particular emphasis on molecular recognition events involved in mRNA translation on the ribosome. The general methodology for efficient free energy calculations is outlined and our specific implementation for binding free energy changes due to base mutations in mRNA and tRNA is described. We show that there are a number of key problems related to the accuracy of protein synthesis that can be addressed with this type of computational approach and several such examples are discussed in detail. These include the decoding of mRNA during peptide chain elongation, initiation and termination of translation, as well as the energetic effects of base tautomerization and tRNA modifications. It is shown that free energy calculations can be made sufficiently reliable to allow quantitative conclusions to be drawn regarding the energetics of cognate versus non-cognate interactions and its structural origins.

Dicyclohexylurea derivatives of amino acids as dye absorbent organogels and anion sensors.

Dicyclohexyl urea (DCU) derivatives of amino acids Fmoc-Phe-DCU (M1), Fmoc-Phg-DCU (M2) and Fmoc-Gaba-DCU (M3) have been shown to form phase selective, thermoreversible and mechanically robust gels in a large range of organic solvents. This is the first report of low molecular weight gelators (LMWG) from DCU derivatives of amino acids. The self-assembly mechanism of the organogels has been probed using concentration dependent 1H NMR, DMSO titration 1H NMR, fluorescence, FTIR, PXRD and FESEM techniques. Self-assembly leading to gelation process is mainly driven by hydrophobicity and π-π stacking interactions in between Fmoc groups. Interestingly, the gels can absorb several kinds of organic dyes efficiently and can be reused for dye absorption for multiple cycles. Additionally, M1-M3 act as sensors for anions like fluoride, acetate and hydroxide, for which they have specific fluorescence response. Gel formation by M1-M3 is completely arrested in the presence of fluoride. The possible binding mode of fluoride has been delineated using DFT studies. Calculations suggest, involvement of urea NH in a six membered intramolecular hydrogen bond, rendering it unavailable for fluoride binding. Backbone -NH of the amino acids of M1-M3 is responsible for fluoride binding. The reported small, economically viable, synthetically facile molecules not only enrich the repertoire of LMWG molecules, but can have multifaceted applications.

Mutations in Parkinson's Disease Associated Protein DJ-1 Alter the Energetics of DJ-1 Dimerization.

Patients suffering from familial Parkinson's disease are linked to mutated DJ-1 protein. Wild-type DJ-1 occurs as a homodimer, which appears to be crucial for its function. It has been established that mutation (L166P) in DJ-1 protein could destabilize the DJ-1 homodimer. Hence, dimerization aspect of DJ-1 is fundamentally important for understanding its link to the disease. X-ray structures of wild-type DJ-1 dimer have given an atomic insight into the interaction network at the dimer interface. However, the energetics of dimerization in the wild-type and its mutant protein is unknown. Using the X-ray structure of wild-type DJ-1 as the template, we report ∼1.55 µs of molecular dynamics simulations to quantitatively estimate the relative free energy of DJ-1 dimerization in the disease linked variant (L166P, A104T, and M26I) with respect to its wild-type analogue. The results suggest that dimerization is disfavored for L166P and A104T mutations, severely for the former. Notably, the M26I mutation does not alter the energetics of DJ-1 dimerization. The dynamics of the DJ-1 dimer is significantly altered in response to the L166P and A104T mutations, resulting in the significant loss of interactions at the dimer interface. L166P mutant showed the structural difference and increased flexibility in α6, α7, α8 regions with respect to the WT. A structural difference in the α6 region was noticeable between WT and A104T mutant of DJ-1. The interaction network in the dimer interface is identical for the wild-type protein and the M26I mutant. No significant change in secondary structural content was observed for DJ-1 mutants (L166P, A104T, M26I) with respect to its WT analogue.

Structure-Based Energetics of Stop Codon Recognition by Eukaryotic Release Factor.

In translation termination, the eukaryotic release factor (eRF1) recognizes mRNA stop codons (UAA, UAG, or UGA) in a ribosomal A site and triggers release of the nascent polypeptide chain from P-site tRNA. eRF1 is highly selective for U in the first position and a combination of purines (except two consecutive guanines, i.e., GG) in the second and third positions. Eukaryotes decode all three stop codons with a single release factor eRF1, instead of two (RF1 and RF2), in bacteria. Furthermore, unlike bacterial RF1/RF2, eRF1 stabilizes the compact U-turn mRNA configuration in the ribosomal A site by accommodating four nucleotides instead of three. Despite the available cryo-EM structures (resolution ∼3.5-3.8 Å), the energetic principle for eRF1 selectivity toward a stop codon remains a fundamentally unsolved problem. Using cryo-EM structures of eukaryotic translation termination complexes as templates, we carried out molecular dynamics free energy simulations of cognate and near-cognate complexes to quantitatively address the energetics of stop codon recognition by eRF1. Our results suggest that eRF1 has a higher discriminatory power against sense codons, compared to that reported earlier for RF1/RF2. The compact mRNA formed specific intra-mRNA interactions, which itself contributed to stop codon specificity. Furthermore, the specificity is enhanced by the loss of protein-mRNA interactions and, most importantly, by desolvation of the incorrect codons in the near-cognate complexes. Our work provides a clue to how eRF1 discriminates between cognate and near-cognate codons during protein synthesis.

Electrostatic free energies in translational GTPases: Classic allostery and the rest.

GTPases typically switch between an inactive, OFF conformation and an active, ON conformation when a GDP ligand is replaced by GTP. Their ON/OFF populations and activity thus depend on the stabilities of four protein complexes, two apo-protein forms, and GTP/GDP in solution. A complete characterization is usually not possible experimentally and poses major challenges for simulations. We review the most important methodological challenges and we review thermodynamic data for two GTPases involved in translation of the genetic code: archaeal Initiation Factors 2 and 5B (aIF2, aIF5B). One main challenge is the multiplicity of states and conformations, including those of GTP/GDP in solution. Another is force field accuracy, especially for interactions of GTP/GDP with co-bound divalent Mg(2+) ions. The calculation of electrostatic free energies also poses specific challenges, and requires careful protocols. For aIF2, experiments and earlier simulations showed that it is a "classic" GTPase, with distinct ON/OFF conformations that prefer to bind GTP and GDP, respectively. For aIF5B, we recently proposed a non-classic mechanism, where the ON/OFF states differ only in the protonation state of Glu81 in the nucleotide binding pocket. This model is characterized here using free energy simulations. The methodological analysis should help future studies, while the aIF2, aIF5B examples illustrate the diversity of ATPase/GTPase mechanisms. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Why base tautomerization does not cause errors in mRNA decoding on the ribosome.

The structure of the genetic code implies strict Watson-Crick base pairing in the first two codon positions, while the third position is known to be degenerate, thus allowing wobble base pairing. Recent crystal structures of near-cognate tRNAs accommodated into the ribosomal A-site, however, show canonical geometry even with first and second position mismatches. This immediately raises the question of whether these structures correspond to tautomerization of the base pairs. Further, if unusual tautomers are indeed trapped why do they not cause errors in decoding? Here, we use molecular dynamics free energy calculations of ribosomal complexes with cognate and near-cognate tRNAs to analyze the structures and energetics of G-U mismatches in the first two codon positions. We find that the enol tautomer of G is almost isoenergetic with the corresponding ketone in the first position, while it is actually more stable in the second position. Tautomerization of U, on the other hand is highly penalized. The presence of the unusual enol form of G thus explains the crystallographic observations. However, the calculations also show that this tautomer does not cause high codon reading error frequencies, as the resulting tRNA binding free energies are significantly higher than for the cognate complex.

Structure-based energetics of mRNA decoding on the ribosome.

The origin of high fidelity in bacterial protein synthesis on the ribosome remains a fundamental unsolved problem despite available three-dimensional structures of different stages of the translation process. However, these structures open up the possibility of directly computing the energetics of tRNA selection that is required for an authentic understanding of fidelity in decoding. Here, we report extensive computer simulations that allow us to quantitatively calculate tRNA discrimination and uncover the energetics underlying accuracy in code translation. We show that the tRNA-mRNA interaction energetics varies drastically along the path from initial selection to peptide bond formation. While the selection process is obviously controlled by kinetics, the underlying thermodynamics explains the origin of the high degree of accuracy. The existence of both low- and high-selectivity states provides an efficient mechanism for initial selection and proofreading that does not require codon-dependent long-range structural signaling within the ribosome. It is instead the distinctly unequal population of the high-selectivity states for cognate and noncognate substrates that is the key discriminatory factor. The simulations reveal the essential roles played both by the 30S subunit conformational switch and by the common tRNA modification at position 37 in amplifying the accuracy.

Energetic tuning by tRNA modifications ensures correct decoding of isoleucine and methionine on the ribosome.

Chemical modifications of tRNAs are critical for accurate translation of the genetic code on the ribosome. The discrimination between isoleucine (AUA) and methionine (AUG) codons depends on such modifications of the wobble position in isoleucine tRNA anticodon loops, in all kingdoms of life. Bacteria and archaea employ functionally similar lysine- and agmatine-conjugated cytidine derivatives to ensure decoding fidelity, but the thermodynamics underlying codon discrimination remains unknown. Here, we report structure-based computer simulations that quantitatively reveal the energetics of this decoding strategy in archaea. The results further show that the agmatidine modification confers tRNA specificity primarily by desolvation of the incorrect codon in the non-cognate complex. Tautomerism is found to play no significant role in this decoding system as the usual amino form of the modified tRNA is by far the most stable.

E-Protein Protonation Titration-Induced Single-Particle Chemical Force Spectroscopy for Microscopic Understanding and pI Estimation of Infectious DENV.

The ionization state of amino acids on the outer surface of a virus regulates its physicochemical properties toward the sorbent surface. Serologically different strains of the dengue virus (DENV) show different extents of infectivity depending upon their interactions with a receptor on the host cell. To understand the structural dependence of E-protein protonation over its sequence dependence, we have followed E-protein titration kinetics both experimentally and theoretically for two differentially infected dengue serotypes, namely, DENV-2 and DENV-4. We have performed E-protein protonation titration-induced single-particle chemical force spectroscopy using an atomic force microscope (AFM) to measure the surface chemistry of DENV in physiological aqueous solutions not only to understand the charge distribution dynamics on the virus surface but also to estimate the isoelectric point (pI) accurately for infectious dengue viruses. Cryo-EM structure-based theoretical pI calculations of the DENV-2 surface protein were shown to be consistent with the evaluated pI value from force spectroscopy measurements. We also highlighted here the role of the microenvironment around the titrable residues (in the 3D-folded structure of the protein) in altering the pKa. This is a comprehensive study to understand how the cumulative charge distribution on the outer surface of a specific serotype of DENV regulates a prominent role of infectivity over minute changes at the genetic level.

Development of protease resistant and non-cytotoxic Jelleine analogs with enhanced broad spectrum antimicrobial efficacy.

Short systemic half- life of Antimicrobial Peptides (AMP) is one of the major bottlenecks that limits their successful commercialization as therapeutics. In this work, we have designed analogs of the natural AMP Jelleine, obtained from royal jelly of apis mellifera. Among the designed peptides, J3 and J4 were the most potent with broad spectrum activities against a varied class of ESKAPE pathogens and fungus C. albicans. All the developed peptides were more effective against Gram-negative bacteria in comparison to the Gram-positive pathogens, and were especially effective against P. aeruginosa and C. albicans.J3 and J4 were completely trypsin resistant and serum stable, while retaining the non-cytotoxicity of the parent Jelleine, Jc. The designed peptides were membranolytic in their mode of action. CD and MD simulations in the presence of bilayers, established that J3 and J4 were non-structured even upon membrane binding and suggested that biological properties of the AMPs were innocent of any specific secondary structural requirements. Enhancement of charge to increase the antimicrobial potency, controlling the hydrophobic-hydrophilic balance to maintain non-cytotoxicity and induction of unnatural amino acid residues to impart protease resistance, remains some of the fundamental principles in the design of more effective antimicrobial therapeutics of the future, which may help combat the quickly rising menace of antimicrobial resistance in the microbes.

Mechanism of Protease Resistance of D-Amino Acid Residue Containing Cationic Antimicrobial Heptapeptides.

Antimicrobial peptides (AMPs) have been an alternate promising class of therapeutics in combating global antibiotic resistance threat. However, the short half-life of AMPs, owing to protease degradability, is one of the major bottlenecks in its commercial success. In this study, we have developed all-D-amino acid containing small cationic peptides P4C and P5C, which are completely protease-resistant, noncytotoxic, nonhemolytic, and potent against the ESKAPE pathogens in comparison to their L analogues. MD simulations suggested marginal improvement in the peptide-binding affinity to the membrane-mimetic SDS micelle (∼ 1 kcal/mol) in response to L → D conversion, corroborating the marginal improvement in the antimicrobial activity. However, L → D chirality conversion severely compromised the peptide:protease (trypsin) binding affinity (≥10 kcal/mol). The relative distance between the scissile peptide carbonyl and the catalytic triad of the protease (H57, D102, and S195) was found to be significantly altered in the D-peptide:protease complex (inactive conformation) relative to the active L-peptide:protease complex. Thus, the poor binding affinity between D-peptides and the protease, resulting in the inactive complex formation, explained their experimentally observed proteolytic stability. This mechanistic insight might be extended to the proteolytic stability of the D-peptides in general and stimulate the rational design of protease-resistant AMPs.

Simulating GTP:Mg and GDP:Mg with a simple force field: a structural and thermodynamic analysis.

Di- and tri-phosphate nucleotides are essential cofactors for many proteins, usually in an Mg(2+) -bound form. Proteins like GTPases often detect the difference between NDP and NTP and respond by changing conformations. To study such complexes, simple, fixed charge force fields have been used, which allow long simulations and precise free energy calculations. The preference for NTP or NDP binding depends on many factors, including ligand structure and Mg(2+) coordination and the changes they undergo upon binding. Here, we use a simple force field to examine two Mg(2+) coordination modes for the unbound GDP and GTP: direct, or "Inner Sphere" (IS) coordination by one or more phosphate oxygens and indirect, "Outer Sphere" (OS) coordination involving one or more bridging waters. We compare GTP: and GDP:Mg binding with OS and IS coordination; combining the results with experimental data then indicates that GTP prefers the latter. We also examine different kinds of IS coordination and their sensitivity to a key force field parameter: the optimal Mg:oxygen van der Waals distance Rmin . Increasing Rmin improves the Mg:oxygen distances, the GTP: and GDP:Mg binding affinities, and the fraction of GTP:Mg with ß + γ phosphate coordination, but does not improve or change the GTP/GDP affinity difference, which remains much larger than experiment. It has no effect on the free energy of GDP binding to a GTPase.

Insights into the Mechanism of CRISPR/Cas9-Based Genome Editing from Molecular Dynamics Simulations.

The CRISPR/Cas9 system is a popular genome-editing tool with immense therapeutic potential. It is a simple two-component system (Cas9 protein and RNA) that recognizes the DNA sequence on the basis of RNA:DNA complementarity, and the Cas9 protein catalyzes the double-stranded break in the DNA. In the past decade, near-atomic resolution structures at various stages of the CRISPR/Cas9 DNA editing pathway have been reported along with numerous experimental and computational studies. Such studies have boosted knowledge of the genome-editing mechanism. Despite such advancements, the application of CRISPR/Cas9 in therapeutics is still limited, primarily due to off-target effects. Several studies aim at engineering high-fidelity Cas9 to minimize the off-target effects. Molecular Dynamics (MD) simulations have been an excellent complement to the experimental studies for investigating the mechanism of CRISPR/Cas9 editing in terms of structure, thermodynamics, and kinetics. MD-based studies have uncovered several important molecular aspects of Cas9, such as nucleotide binding, catalytic mechanism, and off-target effects. In this Review, the contribution of MD simulation to understand the CRISPR/Cas9 mechanism has been discussed, preceded by an overview of the history, mechanism, and structural aspects of the CRISPR/Cas9 system. These studies are important for the rational design of highly specific Cas9 and will also be extremely promising for achieving more accurate genome editing in the future.

Ambidexterity and Left-Handedness Induced by Geminally Disubstituted γ Amino Acid Residues in Chiral 310 Helices.

Chirality is an omnipresent feature in nature's architecture starting from simple molecules like amino acids to complex higher-order structures viz. proteins, DNA, and RNA. The L configuration of proteinogenic amino acids gives rise to right-handed helices. Ambidexterity is as rare in organisms as in molecules. There are only a few reports of ambidexterity in single-peptide molecules composed of either mixed L and D or achiral residues. Here, we report, for the first time, the ambidextrous and left-handed helical conformations in the chiral nonapeptides P1-P3 (Boc-LUVUγx,xULUV-OMe where U = Aib, x,x = 2,2/3,3/4,4), containing chiral L α amino acid residues, in addition to the usually observed right-handed helical conformation. The centrally located achiral γ residue, capable of adopting both left and right-handed helical conformations, induces its handedness on the neighboring chiral and achiral residues, leading to the observation of both left and right-handed helices in P2 and P3. The presence of a single water molecule proximal to the γ residue induces the reversal of helix handedness by forming distinct and stable water-mediated hydrogen bonds. This gives rise to ambidextrous helices as major conformers in P1 and P2. The absence of the observation of ambidexterity in P3 might be due to the inability of γ4,4 in the recruitment of a water molecule. Experiments (NMR, X-ray, and CD) and density functional theory (DFT) calculations suggest that the position of geminal disubstitution is crucial for determining the population of the amenable helical conformations (ambidextrous, left and right-handed) in these chiral peptides.