Conserved C-Terminal Tail Is Responsible for Membrane Localization and Function of Pseudomonas aeruginosa Hemerythrin.

Many bacteria have hemerythrin (Hr) proteins that bind O2, including Pseudomonas aeruginosa, in which microoxia-induced Hr (Mhr) provide fitness advantages under microoxic conditions. Mhr has a 23 amino-acid extension at its C-terminus relative to a well-characterized Hr from Methylococcus capsulatus, and similar extensions are also found in Hrs from other bacteria. The last 11 amino acids of this extended, C-terminal tail are highly conserved in gammaproteobacteria and predicted to form a helix with positively charged and hydrophobic faces. In cellular fractionation assays, wild-type (WT) Mhr was found in both membrane and cytosolic fractions, while a MhrW143* variant lacking the last 11 residues was largely in the cytosol and did not complement Mhr function in competition assays. MhrL112Y, a variant that has a much longer-lived O2-bound form, was fully functional and had a similar localization pattern to that of WT Mhr. Both MhrW143* and MhrL112Y had secondary structures, stabilities, and O2-binding kinetics similar to those of WT Mhr. Fluorescence studies revealed that the C-terminal tail, and particularly the fragment corresponding to its last 11 residues, was sufficient and necessary for association with lipid vesicles. Molecular dynamics simulations and subsequent cellular analysis of Mhr variants have demonstrated that conserved, positively charged residues in the tail are important for Mhr interactions with negatively charged membranes and the contribution of this protein to competitive fitness. Together, these data suggest that peripheral interactions of Mhr with membranes are guided by the C-terminal tail and are independent of O2-binding.

Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases.

Understanding kinase-inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A-1) of the insulin-like growth factor-I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine-arginine-aspartic acid (HRD) and aspartic acid-phenylalanine-glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C-helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.

Interactions and Transport of a Bioconjugated Peptide Targeting the Mitomembrane.

The Szeto-Schiller (SS) peptides are a subclass of cell-penetrating peptides that can specifically target mitochondria and mediate conditions caused by mitochondrial dysfunction. In this work, we constructed an iron-chelating SS peptide and studied its interaction with a mitochondrial-mimicking membrane using atomistic molecular dynamics (MD) simulations. We report that the peptide/membrane interaction is thermodynamically favorable, and the localization of the peptide to the membrane is driven by electrostatic interactions between the cationic residues and the anionic phospholipid headgroups. The insertion of the peptide into the membrane is driven by hydrophobic interactions between the aromatic side chains in the peptide and the lipid acyl tails. We also probed the translocation of the peptide across the membrane by applying nonequilibrium steered MD simulations and resolved the translocation pathway, free energy profile, and metastable states. We explored four distinct orientations of the peptide along the translocation pathway and found that one orientation was energetically more favorable than the other orientations. We tested a significantly slower pulling velocity on the most thermodynamically favorable system and compared metastable states during peptide translocation. We found that the peptide can optimize hydrophobic interactions with the membrane by having aromatic side chains interacting with the lipid acyl tails instead of forming π-π interactions with each other. The mechanistic insights emerging from our work will potentially facilitate improved peptide design with enhanced activity.

Quantitative Assessment of Energetic Contributions of Residues in a SARS-CoV-2 Viral Enzyme/Nanobody Interface.

The highly conserved protease enzyme from SARS-CoV-2 (MPro) is crucial for viral replication and is an attractive target for the design of novel inhibitory compounds. MPro is known to be conformationally flexible and has been stabilized in an extended conformation in a complex with a novel nanobody (NB2B4), which inhibits the dimerization of the enzyme via binding to an allosteric site. However, the energetic contributions of the nanobody residues stabilizing the MPro/nanobody interface remain unresolved. We probed these residues using all-atom MD simulations in combination with alchemical free energy calculations by studying the physical residue-residue interactions and discovered the role of hydrophobic and electrostatic interactions in stabilizing the complex. Specifically, we found via mutational analysis that three interfacial nanobody residues (Y59, R106, and L109) contributed significantly, two residues (L107 and P110) contributed moderately, and two residues (H112 and T113) contributed minimally to the overall binding affinity of the nanobody. We also discovered that the nanobody affinity could be enhanced via a charge-reversal mutation (D62R) that alters the local interfacial electrostatic environment of this residue in the complex. These findings are potentially useful in designing novel synthetic nanobodies as allosteric inhibitors of MPro.

Structural and computational studies of HIV-1 RNA.

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

Structural Models for a Series of Allosteric Inhibitors of IGF1R Kinase.

The allosteric inhibition of insulin-like growth factor receptor 1 kinase (IGF1RK) is a potential strategy to overcome selectivity barriers for targeting receptor tyrosine kinases. We constructed structural models of a series of 12 indole-butyl-amine derivatives that have been reported as allosteric inhibitors of IGF1RK. We further studied the dynamics and interactions of each inhibitor in the allosteric pocket via all-atom explicit-solvent molecular dynamics (MD) simulations. We discovered that a bulky carbonyl substitution at the R1 indole ring is structurally unfavorable for inhibitor binding in the IGF1RK allosteric pocket. Moreover, we found that the most potent derivative (termed C11) acquires a distinct conformation: forming an allosteric pocket channel with better shape complementarity and interactions with the receptor. In addition to a hydrogen-bonding interaction with V1063, the cyano derivative C11 forms a stable hydrogen bond with M1156, which is responsible for its unique binding conformation in the allosteric pocket. Our findings show that the positioning of chemical substituents with different pharmacophore features at the R1 indole ring influences molecular interactions and binding conformations of indole-butyl-amine derivatives and, hence, dramatically affects their potencies. Our results provide a structural framework for the design of allosteric inhibitors with improved affinities and specificities against IGF1RK.

Structural models of viral insulin-like peptides and their analogs.

The insulin receptor (IR), the insulin-like growth factor-1 receptor (IGF1R), and the insulin/IGF1 hybrid receptors (hybR) are homologous transmembrane receptors. The peptide ligands, insulin and IGF1, exhibit significant structural homology and can bind to each receptor via site-1 and site-2 residues with distinct affinities. The variants of the Iridoviridae virus family show capability in expressing single-chain insulin/IGF1 like proteins, termed viral insulin-like peptides (VILPs), which can stimulate receptors from the insulin family. The sequences of VILPs lacking the central C-domain (dcVILPs) are known, but their structures in unbound and receptor-bound states have not been resolved to date. We report all-atom structural models of three dcVILPs (dcGIV, dcSGIV, and dcLCDV1) and their complexes with the receptors (µIR, µIGF1R, and µhybR), and probed the peptide/receptor interactions in each system using all-atom molecular dynamics (MD) simulations. Based on the nonbonded interaction energies computed between each residue of peptides (insulin and dcVILPs) and the receptors, we provide details on residues establishing significant interactions. The observed site-1 insulin/µIR interactions are consistent with previous experimental studies, and a residue-level comparison of interactions of peptides (insulin and dcVILPs) with the receptors revealed that, due to sequence differences, dcVILPs also establish some interactions distinct from those between insulin and IR. We also designed insulin analogs and report enhanced interactions between some analogs and the receptors.

Mechanism of Ligand Discrimination by the NMT1 Riboswitch.

Riboswitches are conserved functional domains in mRNA that almost exclusively exist in bacteria. They regulate the biosynthesis and transport of amino acids and essential metabolites such as coenzymes, nucleobases, and their derivatives by specifically binding small molecules. Due to their ability to precisely discriminate between different cognate molecules as well as their common existence in bacteria, riboswitches have become potential antibacterial drug targets that could deliver urgently needed antibiotics with novel mechanisms of action. In this work, we report the recognition mechanisms of four oxidization products (XAN, AZA, UAC, and HPA) generated during purine degradation by an RNA motif termed the NMT1 riboswitch. Specifically, we investigated the physical interactions between the riboswitch and the oxidized metabolites by computing the changes in the free energy on mutating key nucleobases in the ligand binding pocket of the riboswitch. We discovered that the electrostatic interactions are central to ligand discrimination by this riboswitch. The relative binding free energies of the mutations further indicated that some of the mutations can also strengthen the binding affinities of the ligands (AZA, UAC, and HPA). These mechanistic details are also potentially relevant in the design of novel compounds targeting riboswitches.

Role of Dynamics and Mutations in Interactions of a Zinc Finger Antiviral Protein with CG-rich Viral RNA.

Zinc finger antiviral protein (ZAP) is a host antiviral factor that selectively inhibits the replication of a variety of viruses. ZAP recognizes the CG-enriched RNA sequences and activates the viral RNA degradation machinery. In this work, we investigated the dynamics of a ZAP/RNA complex and computed the energetics of mutations in ZAP that affect its binding to the viral RNA. The crystal structure of a mouse-ZAP/RNA complex showed that RNA interacts with the zinc finger 2 (ZF2) and ZF3 domains. However, we found that due to the dynamic behavior of the single-stranded RNA, the terminal nucleotides C1 and G2 of RNA change their positions from the ZF3 to the ZF1 domain. Moreover, the electrostatic interactions between the zinc ions and the viral RNA provide further stability to the ZAP/RNA complex. We also provide structural and thermodynamic evidence for seven residue pairs (C1-Arg74, C1-Arg179, G2-Arg74, U3-Lys76, C4-Lys76, G5-Arg95, and U6-Glu204) that show favorable ZAP/RNA interactions, although these interactions were not observed in the ZAP/RNA crystal structure. Consistent with the observations from the mouse-ZAP/RNA crystal structure, we found that four residue pairs (C4-Lys89, C4-Leu90, C4-Tyr108, and G5-Lys107) maintained stable interactions in MD simulations. Based on experimental mutagenesis studies and our residue-level interaction analysis, we chose seven residues (Arg74, Lys76, Lys89, Arg95, Lys107, Tyr108, and Arg179) for individual alanine mutations. In addition, we studied mutations in those residues that are only observed in the crystal structures as interacting with RNA (Tyr98, Glu148, and Arg170). Out of these 10 mutations, we found that the Ala mutation in each of the five residues Arg74, Lys76, Lys89, Lys107, and Glu148 significantly reduced the binding affinity of ZAP to RNA.

Structures and interactions of insulin-like peptides from cone snail venom.

The venomous insulin-like peptides released by certain cone snails stimulate hypoglycemic shock to immobilize fish and catch the prey. Compared to human insulin (hIns), the cone snail insulins (Con-Ins) are typically monomeric and shorter in sequence, yet they exhibit moderate hIns-like biological activity. We have modeled six variants of Con-Ins (G3, K1, K2, T1A, T1B, and T2) and carried out explicit-solvent molecular dynamics (MD) simulations of eight types of insulins, two with known structures (hIns and Con-Ins-G1) and six Con-Ins with modeled structures, to characterize key residues of each insulin that interact with the truncated human insulin receptor (µIR). We show that each insulin/µIR complex is stable during explicit-solvent MD simulations and hIns interactions indicate the highest affinity for the "site 1" of IR. The residue contact maps reveal that each insulin preferably interacts with the αCT peptide than the L1 domain of IR. Through analysis of the average nonbonded interaction energy contribution of every residue of each insulin for the µIR, we probe the residues establishing favorable interactions with the receptor. We compared the interaction energy of each residue of every Con-Ins to the µIR and observed that γ-carboxylated glutamate (Gla), His, Thr, Tyr, Tyr/His, and Asn in Con-Ins are favorable substitutions for GluA4, AsnA21, ValB12, LeuB15, GlyB20, and ArgB22 in hIns, respectively. The identified insulin analogs, although lacking the last eight residues of the B-chain of hIns, bind strongly to µIR. Our findings are potentially useful in designing potent fast-acting therapeutic insulin.

Molecular interactions and inhibition of the SARS-CoV-2 main protease by a thiadiazolidinone derivative.

We report molecular interactions and inhibition of the main protease (MPro ) of SARS-CoV-2, a key enzyme involved in the viral life cycle. By using a thiadiazolidinone (TDZD) derivative as a chemical probe, we explore the conformational dynamics of MPro via docking protocols and molecular dynamics simulations in all-atom detail. We reveal the local and global dynamics of MPro in the presence of this inhibitor and confirm the inhibition of the enzyme with an IC50 value of 1.39 ± 0.22 µM, which is comparable to other known inhibitors of this enzyme.

Role of Mutations in Differential Recognition of Viral RNA Molecules by Peptides.

The conserved noncoding RNA elements in viral genomes interact with proteins to regulate various events during viral replication. We report studies on the recognition mechanisms of two helical peptides, namely, a native (Rev) peptide and a lab-evolved (RSG1.2) peptide, by a highly conserved viral RNA element from the human immunodeficiency virus 1 genome. Specifically, we investigated the physical interactions between the viral RNA molecule and helical peptides by computing free energy changes on mutating key amino acid residues involved in recognition of an internal loop in the viral RNA molecule.

Allosteric Pathways Originating at Cysteine Residues in Regulators of G-Protein Signaling Proteins.

Regulators of G-protein signaling (RGS) proteins play a central role in modulating signaling via G-protein coupled receptors (GPCRs). Specifically, RGS proteins bind to activated Gα subunits in G-proteins, accelerate the GTP hydrolysis, and thereby rapidly dampen GPCR signaling. Therefore, covalent molecules targeting conserved cysteine residues among RGS proteins have emerged as potential candidates to inhibit the RGS/Gα protein-protein interaction and enhance GPCR signaling. Although these inhibitors bind to conserved cysteine residues among RGS proteins, we have previously suggested [J. Am. Chem. Soc. 2018;140:3454-3460] that their potencies and specificities are related to differential protein dynamics among RGS proteins. Using data from all-atom molecular dynamics simulations, we reveal these differences in dynamics of RGS proteins by partitioning the protein structural space into a network of communities that allow allosteric signals to propagate along unique pathways originating at inhibitor binding sites and terminating at the RGS/Gα protein-protein interface.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides.

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

Conformational dynamics and energetics of viral RNA recognition by lab-evolved proteins.

The conserved and structured elements in viral RNA genomes interact with proteins to regulate various events in the viral life cycle and have become key targets for developing novel therapeutic approaches. We probe physical interactions between lab-evolved proteins and a viral RNA element from the HIV-1 genome. Specifically, we study the role of an arginine-rich loop in recognition of designed proteins by the viral RNA element. We report free energy calculations to quantitatively estimate the protein/RNA binding energetics, focusing on the mutations of arginine residues involved in recognition of the major groove of RNA by proteins.

Role of conformational heterogeneity in ligand recognition by viral RNA molecules.

Ribonucleic acid (RNA) molecules are known to undergo conformational changes in response to various environmental stimuli including temperature, pH, and ligands. In particular, viral RNA molecules are a key example of conformationally adapting molecules that have evolved to switch between many functional conformations. The transactivation response element (TAR) RNA from the type-1 human immunodeficiency virus (HIV-1) is a viral RNA molecule that is being increasingly explored as a potential therapeutic target due to its role in the viral replication process. In this work, we have studied the dynamics in TAR RNA in apo and liganded states by performing explicit-solvent molecular dynamics (MD) simulations initiated with 27 distinct structures. We determined that the TAR RNA structure is significantly stabilized on ligand binding with especially decreased fluctuations in its two helices. This rigidity is further coupled with the decreased flipping of bulge nucleotides, which were observed to flip more frequently in the absence of ligands. We found that initially-distinct structures of TAR RNA converged to similar conformations on removing ligands. We also report that conformational dynamics in unliganded TAR structures leads to the formation of binding pockets capable of accommodating ligands of various sizes.

Self-assembly behavior of experimentally realizable lobed patchy particles.

We report simulation studies on the self-assembly behavior of five different types of lobed patchy particles of different shapes (snowman, dumbbell, trigonal planar, square planar, and tetrahedral). Inspired by an experimental method of synthesizing patchy particles (Wang et al., Nature, 2012, 491, 51-55), we control the lobe size indirectly by gradually varying the seed diameter and study its effect on self-assembled structures at different temperatures. Snowman shaped particles self-assemble only at a lower temperature and form two-dimensional sheets, elongated micelles, and spherical micelles, depending on the seed diameter. Each of the four other lobed particles self-assemble into four distinct morphologies (random aggregates, spherical aggregates, liquid droplets, and crystalline structures) for a given lobe size and temperature. We observed temperature-dependent transitions between two morphologies depending on the type of the lobed particle. The self-assembled structures formed by these four types of particles are porous. We show that their porosities can be tuned by controlling the lobe size and temperature.

Self-assembly of lobed particles into amorphous and crystalline porous structures.

We report simulation studies on the self-assembly of hard-lobed particles (patchy particles where patches appear as lobes around a seed) of different shapes and show that various types of self-assembled morphologies can be achieved by tuning inter-lobe interactions. On self-assembly, the linear building blocks having two lobes around the seed formed rings, the trigonal planar building blocks formed cylindrical hollow tubes and two-dimensional sheets, and the square planar building blocks formed spherical clathrates. The tetrahedral, trigonal bipyramidal, and the octahedral-shaped particles formed compact porous crystalline structures which are constituted by either hexagonal close packed or face centered cubic lattices. The pore size distributions revealed that linear, trigonal planar, and square planar building blocks create highly porous self-assembled structures. Our results suggest that these self-assembled morphologies will potentially find applications in tissue engineering, host-guest chemistry, adsorption, and catalysis.

Role of Entropy in Colloidal Self-Assembly.

Entropy plays a key role in the self-assembly of colloidal particles. Specifically, in the case of hard particles, which do not interact or overlap with each other during the process of self-assembly, the free energy is minimized due to an increase in the entropy of the system. Understanding the contribution of entropy and engineering it is increasingly becoming central to modern colloidal self-assembly research, because the entropy serves as a guide to design a wide variety of self-assembled structures for many technological and biomedical applications. In this work, we highlight the importance of entropy in different theoretical and experimental self-assembly studies. We discuss the role of shape entropy and depletion interactions in colloidal self-assembly. We also highlight the effect of entropy in the formation of open and closed crystalline structures, as well as describe recent advances in engineering entropy to achieve targeted self-assembled structures.

An Interhelical Salt Bridge Controls Flexibility and Inhibitor Potency for Regulators of G-protein Signaling Proteins 4, 8, and 19.

Regulators of G-protein signaling (RGS) proteins modulate receptor signaling by binding to activated G-protein α-subunits, accelerating GTP hydrolysis. Selective inhibition of RGS proteins increases G-protein activity and may provide unique tissue specificity. Thiadiazolidinones (TDZDs) are covalent inhibitors that act on cysteine residues to inhibit RGS4, RGS8, and RGS19. There is a correlation between protein flexibility and potency of inhibition by the TDZD 4-[(4- fluorophenyl)methyl]-2-(4-methylphenyl)-1,2,4-thiadiazolidine-3,5-dione (CCG-50014). In the context of a single conserved cysteine residue on the α 4 helix, RGS19 is the most flexible and most potently inhibited by CCG-50014, followed by RGS4 and RGS8. In this work, we identify residues responsible for differences in both flexibility and potency of inhibition among RGS isoforms. RGS19 lacks a charged residue on the α 4 helix that is present in RGS4 and RGS8. Introducing a negative charge at this position (L118D) increased the thermal stability of RGS19 and decreased the potency of inhibition of CCG-50014 by 8-fold. Mutations eliminating salt bridge formation in RGS8 and RGS4 decreased thermal stability in RGS8 and increased potency of inhibition of both RGS4 and RGS8 by 4- and 2-fold, respectively. Molecular dynamics simulations with an added salt bridge in RGS19 (L118D) showed reduced RGS19 flexibility. Hydrogen-deuterium exchange studies showed striking differences in flexibility in the α 4 helix of RGS4, 8, and 19 with salt bridge-modifying mutations. These results show that the α 4 salt bridge-forming residue controls flexibility in several RGS isoforms and supports a causal relationship between RGS flexibility and the potency of TDZD inhibitors. SIGNIFICANCE STATEMENT: Inhibitor potency is often viewed in relation to the static structure of a target protein binding pocket. Using both experimental and computation studies we assess determinants of dynamics and inhibitor potency for three different RGS proteins. A single salt bridge-forming residue determines differences in flexibility between RGS isoforms; mutations either increase or decrease protein motion with correlated alterations in inhibitor potency. This strongly suggests a causal relationship between RGS protein flexibility and covalent inhibitor potency.