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.

Spontaneous Dimerization and Distinct Packing Modes of Transmembrane Domains in Receptor Tyrosine Kinases.

The insulin receptor (IR) and the insulin-like growth factor-1 receptor (IGF1R) are homodimeric transmembrane glycoproteins that transduce signals across the membrane on binding of extracellular peptide ligands. The structures of IR/IGF1R fragments in apo and liganded states have revealed that the extracellular subunits of these receptors adopt Λ-shaped configurations to which are connected the intracellular tyrosine kinase (TK) domains. The binding of peptide ligands induces structural transitions in the extracellular subunits leading to potential dimerization of transmembrane domains (TMDs) and autophosphorylation in TKs. However, the activation mechanisms of IR/IGF1R, especially the role of TMDs in coordinating signal-inducing structural transitions, remain poorly understood, in part due to the lack of structures of full-length receptors in apo or liganded states. While atomistic simulations of IR/IGF1R TMDs showed that these domains can dimerize in single component membranes, spontaneous unbiased dimerization in a plasma membrane having physiologically representative lipid composition has not been observed. We address this limitation by employing coarse-grained (CG) molecular dynamics simulations to probe the dimerization propensity of IR/IGF1R TMDs. We observed that TMDs in both receptors spontaneously dimerized independent of their initial orientations in their dissociated states, signifying their natural propensity for dimerization. In the dimeric state, IR TMDs predominantly adopted X-shaped configurations with asymmetric helical packing and significant tilt relative to the membrane normal, while IGF1R TMDs adopted symmetric V-shaped or parallel configurations with either no tilt or a small tilt relative to the membrane normal. Our results suggest that IR/IGF1R TMDs spontaneously dimerize and adopt distinct dimerized configurations.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop.

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

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.

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.

Molecular architecture and conservation of an immature human endogenous retrovirus.

The human endogenous retrovirus K (HERV-K) is the most recently acquired endogenous retrovirus in the human genome and is activated and expressed in many cancers and amyotrophic lateral sclerosis. We present the immature HERV-K capsid structure at 3.2 Å resolution determined from native virus-like particles using cryo-electron tomography and subtomogram averaging. The structure shows a hexamer unit oligomerized through a 6-helix bundle, which is stabilized by a small molecule analogous to IP6 in immature HIV-1 capsid. The HERV-K immature lattice is assembled via highly conserved dimer and trimer interfaces, as detailed through all-atom molecular dynamics simulations and supported by mutational studies. A large conformational change mediated by the linker between the N-terminal and the C-terminal domains of CA occurs during HERV-K maturation. Comparison between HERV-K and other retroviral immature capsid structures reveals a highly conserved mechanism for the assembly and maturation of retroviruses across genera and evolutionary time.

Molecular architecture and conservation of an immature human endogenous retrovirus.

A significant part of the human genome consists of endogenous retroviruses sequences. Human endogenous retrovirus K (HERV-K) is the most recently acquired endogenous retrovirus, is activated and expressed in many cancers and amyotrophic lateral sclerosis and possibly contributes to the aging process. To understand the molecular architecture of endogenous retroviruses, we determined the structure of immature HERV-K from native virus-like particles (VLPs) using cryo-electron tomography and subtomogram averaging (cryoET STA). The HERV-K VLPs show a greater distance between the viral membrane and immature capsid lattice, correlating with the presence of additional peptides, SP1 and p15, between the capsid (CA) and matrix (MA) proteins compared to the other retroviruses. The resulting cryoET STA map of the immature HERV-K capsid at 3.2 Å resolution shows a hexamer unit oligomerized through a 6-helix bundle which is further stabilized by a small molecule in the same way as the IP6 in immature HIV-1 capsid. The HERV-K immature CA hexamer assembles into the immature lattice via highly conserved dimmer and trimer interfaces, whose interactions were further detailed through all-atom molecular dynamics simulations and supported by mutational studies. A large conformational change mediated by the flexible linker between the N-terminal and the C-terminal domains of CA occurs between the immature and the mature HERV-K capsid protein, analogous to HIV-1. Comparison between HERV-K and other retroviral immature capsid structures reveals a highly conserved mechanism for the assembly and maturation of retroviruses across genera and evolutionary time.

The capsid lattice engages a bipartite NUP153 motif to mediate nuclear entry of HIV-1 cores.

Increasing evidence has suggested that the HIV-1 capsid enters the nucleus in a largely assembled, intact form. However, not much is known about how the cone-shaped capsid interacts with the nucleoporins (NUPs) in the nuclear pore for crossing the nuclear pore complex. Here, we elucidate how NUP153 binds HIV-1 capsid by engaging the assembled capsid protein (CA) lattice. A bipartite motif containing both canonical and noncanonical interaction modules was identified at the C-terminal tail region of NUP153. The canonical cargo-targeting phenylalanine-glycine (FG) motif engaged the CA hexamer. By contrast, a previously unidentified triple-arginine (RRR) motif in NUP153 targeted HIV-1 capsid at the CA tri-hexamer interface in the capsid. HIV-1 infection studies indicated that both FG- and RRR-motifs were important for the nuclear import of HIV-1 cores. Moreover, the presence of NUP153 stabilized tubular CA assemblies in vitro. Our results provide molecular-level mechanistic evidence that NUP153 contributes to the entry of the intact capsid into the nucleus.

The Drug-Induced Interface That Drives HIV-1 Integrase Hypermultimerization and Loss of Function.

Allosteric HIV-1 integrase (IN) inhibitors (ALLINIs) are an emerging class of small molecules that disrupt viral maturation by inducing the aberrant multimerization of IN. Here, we present cocrystal structures of HIV-1 IN with two potent ALLINIs, namely, BI-D and the drug candidate Pirmitegravir. The structures reveal atomistic details of the ALLINI-induced interface between the HIV-1 IN catalytic core and carboxyl-terminal domains (CCD and CTD). Projecting from their principal binding pocket on the IN CCD dimer, the compounds act as molecular glue by engaging a triad of invariant HIV-1 IN CTD residues, namely, Tyr226, Trp235, and Lys266, to nucleate the CTD-CCD interaction. The drug-induced interface involves the CTD SH3-like fold and extends to the beginning of the IN carboxyl-terminal tail region. We show that mutations of HIV-1 IN CTD residues that participate in the interface with the CCD greatly reduce the IN-aggregation properties of Pirmitegravir. Our results explain the mechanism of the ALLINI-induced condensation of HIV-1 IN and provide a reliable template for the rational development of this series of antiretrovirals through the optimization of their key contacts with the viral target. IMPORTANCE Despite the remarkable success of combination antiretroviral therapy, HIV-1 remains among the major causes of human suffering and loss of life in poor and developing nations. To prevail in this drawn-out battle with the pandemic, it is essential to continue developing advanced antiviral agents to fight drug resistant HIV-1 variants. Allosteric integrase inhibitors (ALLINIs) are an emerging class of HIV-1 antagonists that are orthogonal to the current antiretroviral drugs. These small molecules act as highly specific molecular glue, which triggers the aggregation of HIV-1 integrase. In this work, we present high-resolution crystal structures that reveal the crucial interactions made by two potent ALLINIs, namely, BI-D and Pirmitegravir, with HIV-1 integrase. Our results explain the mechanism of drug action and will inform the development of this promising class of small molecules for future use in antiretroviral regimens.

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.

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.

Reaction Coordinate and Thermodynamics of Base Flipping in RNA.

Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.

Ligand Recognition in Viral RNA Necessitates Rare Conformational Transitions.

Ribonucleic acids (RNAs) are conformationally flexible molecules that fold into three-dimensional structures and play an important role in different cellular processes as well as in the development of many diseases. RNA has therefore become an important target for developing novel therapeutic approaches. The biophysical processes underlying RNA function are often associated with rare structural transitions that play a key role in ligand recognition. In this work, we probe these rarely occurring transitions using nonequilibrium simulations by characterizing the dissociation of a ligand molecule from an HIV-1 viral RNA element. Specifically, we observed base-flipping rare events that are coupled with ligand binding/unbinding and also provided mechanistic details underlying these transitions.