An evolutionarily conserved N-terminal leucine is essential for MX1 GTPase antiviral activity against different families of RNA viruses.

Myxovirus resistance protein 1 (MX1) and MX2 are homologous, dynamin-like large GTPases, induced upon interferon exposure. Human MX1 (HsMX1) is known to inhibit many viruses, including influenza A virus, by likely acting at various steps of their life cycles. Despite decades of studies, the mechanism(s) of action with which MX1 proteins manage to inhibit target viruses is not fully understood. MX1 proteins are mechano-enzymes and share a similar organization to dynamin, with a GTPase domain and a carboxy-terminal stalk domain, connected by a bundle signaling element. These three elements are known to be essential for antiviral activity. HsMX1 has two unstructured regions, the L4 loop, also essential for antiviral activity, and a short amino (N)-terminal region, which greatly varies between MX1 proteins of different species. The role of this N-terminal domain in antiviral activity is not known. Herein, using mutagenesis, imaging, and biochemical approaches, we demonstrate that the N-terminal domain of HsMX1 is essential for antiviral activity against influenza A virus, Vesicular Stomatitis Virus, and La Crosse virus. Furthermore, we pinpoint a highly conserved leucine within this region, which is absolutely crucial for human, mouse, and bat MX1 protein antiviral activity. Importantly, mutation of this leucine does not compromise GTPase activity or oligomerization capabilities but does modify MX1 protein subcellular localization. The discovery of this essential and highly conserved residue defines this region as key for antiviral activity and may reveal insights as to the mechanism(s) of action of MX1 proteins.

Singular Interface Dynamics of the SARS-CoV-2 Delta Variant Explained with Contact Perturbation Analysis.

Emerging SARS-CoV-2 variants raise concerns about our ability to withstand the Covid-19 pandemic, and therefore, understanding mechanistic differences of those variants is crucial. In this study, we investigate disparities between the SARS-CoV-2 wild type and five variants that emerged in late 2020, focusing on the structure and dynamics of the spike protein interface with the human angiotensin-converting enzyme 2 (ACE2) receptor, by using crystallographic structures and extended analysis of microsecond molecular dynamics simulations. Dihedral angle principal component analysis (PCA) showed the strong similarities in the spike receptor binding domain (RBD) dynamics of the Alpha, Beta, Gamma, and Delta variants, in contrast with those of WT and Epsilon. Dynamical perturbation networks and contact PCA identified the peculiar interface dynamics of the Delta variant, which cannot be directly imputable to its specific L452R and T478K mutations since those residues are not in direct contact with the human ACE2 receptor. Our outcome shows that in the Delta variant the L452R and T478K mutations act synergistically on neighboring residues to provoke drastic changes in the spike/ACE2 interface; thus a singular mechanism of action eventually explains why it dominated over preceding variants.

Dihydropyrimidinase protects from DNA replication stress caused by cytotoxic metabolites.

Imbalance in the level of the pyrimidine degradation products dihydrouracil and dihydrothymine is associated with cellular transformation and cancer progression. Dihydropyrimidines are degraded by dihydropyrimidinase (DHP), a zinc metalloenzyme that is upregulated in solid tumors but not in the corresponding normal tissues. How dihydropyrimidine metabolites affect cellular phenotypes remains elusive. Here we show that the accumulation of dihydropyrimidines induces the formation of DNA-protein crosslinks (DPCs) and causes DNA replication and transcriptional stress. We used Xenopus egg extracts to recapitulate DNA replication invitro. We found that dihydropyrimidines interfere directly with the replication of both plasmid and chromosomal DNA. Furthermore, we show that the plant flavonoid dihydromyricetin inhibits human DHP activity. Cellular exposure to dihydromyricetin triggered DPCs-dependent DNA replication stress in cancer cells. This study defines dihydropyrimidines as potentially cytotoxic metabolites that may offer an opportunity for therapeutic-targeting of DHP activity in solid tumors.

4-Substituted-1,2,3-triazolo nucleotide analogues as CD73 inhibitors, their synthesis, in vitro screening, kinetic and in silico studies.

Three series of nucleotide analogues were synthesized and evaluated as potential CD73 inhibitors. Nucleobase replacement consisted in connecting the appropriate aromatic or purine residues through a triazole moiety that is generated from 1,3-dipolar cycloaddition. The first series is related to 4-substituted-1,2,3-triazolo-ß-hydroxyphosphonate ribonucleosides. Additional analogues were also obtained, in which the phosphonate group was replaced by a bisphosphonate pattern (P-C-P-C, series 2) or the ribose moiety was removed leading to acyclic derivatives (series 3). The ß-hydroxyphosphonylphosphonate ribonucleosides (series 2) were found to be potent inhibitors of CD73 using both purified recombinant protein and cell-based assays. Two compounds (2a and 2b) that contained a bis(trifluoromethyl)phenyl or a naphthyl substituents proved to be the most potent inhibitors, with IC50 values of 4.8 ± 0.8 µM and 0.86 ± 0.2 µM, compared to the standard AOPCP (IC50 value of 3.8 ± 0.9 µM), and were able to reverse the adenosine-mediated immune suppression on human T cells. This series of compounds illustrates a new type of CD73 inhibitors.

Identification of allosteric inhibitors of the ecto-5'-nucleotidase (CD73) targeting the dimer interface.

The ecto-5'-nucleotidase CD73 plays an important role in the production of immune-suppressive adenosine in tumor micro-environment, and has become a validated drug target in oncology. Indeed, the anticancer immune response involves extracellular ATP to block cell proliferation through T-cell activation. However, in the tumor micro-environment, two extracellular membrane-bound enzymes (CD39 and CD73) are overexpressed and hydrolyze efficiently ATP into AMP then further into immune-suppressive adenosine. To circumvent the impact of CD73-generated adenosine, we applied an original bioinformatics approach to identify new allosteric inhibitors targeting the dimerization interface of CD73, which should impair the large dynamic motions required for its enzymatic function. Several hit compounds issued from virtual screening campaigns showed a potent inhibition of recombinant CD73 with inhibition constants in the low micromolar range and exhibited a non-competitive inhibition mode. The structure-activity relationships studies indicated that several amino acid residues (D366, H456, K471, Y484 and E543 for polar interactions and G453-454, I455, H456, L475, V542 and G544 for hydrophobic contacts) located at the dimerization interface are involved in the tight binding of hit compounds and likely contributed for their inhibitory activity. Overall, the gathered information will guide the upcoming lead optimization phase that may lead to potent and selective CD73 inhibitors, able to restore the anticancer immune response.

An evolutionary conserved zinc finger protein is involved in Toxoplasma gondii mRNA nuclear export.

Apicomplexan parasites are responsible for some of the most deadly parasitic diseases affecting humans and livestock. There is an urgent need for new medicines that will target apicomplexan-specific pathways. We characterized a Toxoplasma gondii C2H2 zinc finger protein, named TgZNF2, which is conserved among eukaryotes. We constructed an inducible KO strain (iKO-TgZNF2) for this gene where the tgznf2 gene expression is repressed in the presence of a tetracycline analog (ATc). We showed that the iKO-TgZNF2 parasites are unable to proliferate after depletion of the TgZNF2 protein. Complementation with a full length copy of the gene restores the phenotype Moreover, the homolog of this protein in the related apicomplexan Plasmodium falciparum was shown to efficiently rescue the phenotype, suggesting that this pathway is likely conserved among apicomplexan parasites. We demonstrated that the iKO-mutant lacking TgZNF2 are arrested during the cell cycle during the G1 phase. We identified potential protein partners of this protein among which are spliceosomal complex and mRNA nuclear export components. We confirmed that TgZNF2 is able to bind in vivo to transcripts but splicing is not perturbed in the ATc-treated parasites. Instead, we demonstrated that TgZNF2 depletion leads to the sequestration of polyA+ mRNAs in the nucleus while ribosomal RNAs are not affected. We discovered a conserved protein with specific apicomplexan functional properties that is essential for the survival of T. gondii. TgZNF2 may be crucial to ensure the correct polyA+ mRNA nuclear export, a function that is conserved in P. falciparum.

Aminoglycoside binding and catalysis specificity of aminoglycoside 2″-phosphotransferase IVa: A thermodynamic, structural and kinetic study.

BACKGROUND: Aminoglycoside O-phosphotransferases make up a large class of bacterial enzymes that is widely distributed among pathogens and confer a high resistance to several clinically used aminoglycoside antibiotics. Aminoglycoside 2″-phosphotransferase IVa, APH(2″)-IVa, is an important member of this class, but there is little information on the thermodynamics of aminoglycoside binding and on the nature of its rate-limiting step. METHODS: We used isothermal titration calorimetry, electrostatic potential calculations, molecular dynamics simulations and X-ray crystallography to study the interactions between the enzyme and different aminoglycosides. We determined the rate-limiting step of the reaction by the means of transient kinetic measurements. RESULTS: For the first time, Kd values were determined directly for APH(2″)-IVa and different aminoglycosides. The affinity of the enzyme seems to anti-correlate with the molecular weight of the ligand, suggesting a limited degree of freedom in the binding site. The main interactions are electrostatic bonds between the positively charged amino groups of aminoglycosides and Glu or Asp residues of APH. In spite of the significantly different ratio Kd/Km, there is no large difference in the transient kinetics obtained with the different aminoglycosides. We show that a product release step is rate-limiting for the overall reaction. CONCLUSIONS: APH(2″)-IVa has a higher affinity for aminoglycosides carrying an amino group in 2' and 6', but tighter bindings do not correlate with higher catalytic efficiencies. As with APH(3')-IIIa, an intermediate containing product is preponderant during the steady state. GENERAL SIGNIFICANCE: This intermediate may constitute a good target for future drug design.

Kinetic characterization and molecular docking of novel allosteric inhibitors of aminoglycoside phosphotransferases.

BACKGROUND: Bacterial antibiotic resistance often leads to treatment failure which may have serious consequences, especially in critically sick patients. Resistance to aminoglycosides is mainly due to the expression of antibiotic-modifying enzymes. One important mechanism of aminoglycoside modification is the ATP/GTP-dependent O-phosphorylation catalyzed by aminoglycoside phosphotransferases, APHs. The aim of this study is to identify specific inhibitors of APHs that could restore bacterial susceptibility to aminoglycosides. METHODS: We focused on the search for allosteric inhibitors that bind to small cavities of the protein and block the enzyme function by perturbing its dynamics. RESULTS: From normal mode analysis, a cavity of variable volume belonging to a large groove which splits the protein into two parts was chosen as target. By molecular docking, we screened a large library of commercially available compounds. Seventeen of the highest ranked compounds were tested by in vitro kinetic experiments in order to evaluate their ability to inhibit APHs. Site-directed mutagenesis was carried out with the aim of confirming the inhibition mechanism determined kinetically and the interactions with the protein predicted by in silico studies. These interactions were also confirmed by the use of structurally-related molecules. CONCLUSIONS: Two compounds showed interesting inhibition properties, and one was able to block two different classes of APH. GENERAL SIGNIFICANCE: This study gives new insights into the inhibition of APHs by such allosteric inhibitors, and provides the basis for the future development of combined therapies, antibiotic plus APH inhibitor, which may reverse the resistance to aminoglycosides in a clinical context.

Neurotoxicity of a Biopesticide Analog on Zebrafish Larvae at Nanomolar Concentrations.

Despite the ever-increasing role of pesticides in modern agriculture, their deleterious effects are still underexplored. Here we examine the effect of A6, a pesticide derived from the naturally-occurring α-terthienyl, and structurally related to the endocrine disrupting pesticides anilinopyrimidines, on living zebrafish larvae. We show that both A6 and an anilinopyrimidine, cyprodinyl, decrease larval survival and affect central neurons at micromolar concentrations. Focusing on a superficial and easily observable sensory system, the lateral line system, we found that defects in axonal and sensory cell regeneration can be observed at much lower doses, in the nanomolar range. We also show that A6 accumulates preferentially in lateral line neurons and hair cells. We examined whether A6 affects the expression of putative target genes, and found that genes involved in apoptosis/cell proliferation are down-regulated, as well as genes reflecting estrogen receptor activation, consistent with previous reports that anilinopyrimidines act as endocrine disruptors. On the other hand, canonical targets of endocrine signaling are not affected, suggesting that the neurotoxic effect of A6 may be due to the binding of this compound to a recently identified, neuron-specific estrogen receptor.

Beta-hydroxyphosphonate ribonucleoside analogues derived from 4-substituted-1,2,3-triazoles as IMP/GMP mimics: synthesis and biological evaluation.

A series of seventeen ß-hydroxyphosphonate ribonucleoside analogues containing 4-substituted-1,2,3-triazoles was synthesized and fully characterized. Such compounds were designed as potential inhibitors of the cytosolic 5'-nucleotidase II (cN-II), an enzyme involved in the regulation of purine nucleotide pools. NMR and molecular modelling studies showed that a few derivatives adopted similar structural features to IMP or GMP. Five derivatives were identified as modest inhibitors with 53 to 64% of cN-II inhibition at 1 mM.

Restricted mobility of side chains on concave surfaces of solenoid proteins may impart heightened potential for intermolecular interactions.

Significant progress has been made in the determination of the protein structures with their number today passing over a hundred thousand structures. The next challenge is the understanding and prediction of protein-protein and protein-ligand interactions. In this work we address this problem by analyzing curved solenoid proteins. Many of these proteins are considered as "hub molecules" for their high potential to interact with many different molecules and to be a scaffold for multisubunit protein machineries. Our analysis of these structures through molecular dynamics simulations reveals that the mobility of the side-chains on the concave surfaces of the solenoids is lower than on the convex ones. This result provides an explanation to the observed preferential binding of the ligands, including small and flexible ligands, to the concave surface of the curved solenoid proteins. The relationship between the landscapes and dynamic properties of the protein surfaces can be further generalized to the other types of protein structures and eventually used in the computer algorithms, allowing prediction of protein-ligand interactions by analysis of protein surfaces.

Coarse-grained simulations of the HIV-1 matrix protein anchoring: revisiting its assembly on membrane domains.

In the accepted model for human immunodeficiency virus preassembly in infected host cells, the anchoring to the intracellular leaflet of the membrane of the matrix domain (MA) that lies at the N-terminus of the viral Gag protein precursor appears to be one of the crucial steps for particle assembly. In this study, we simulated the membrane anchoring of human immunodeficiency virus-1 myristoylated MA protein using a coarse-grained representation of both the protein and the membrane. Our calculations first suggest that the myristoyl group could spontaneously release from its initial hydrophobic pocket before MA protein interacts with the lipid membrane. All-atom simulations confirmed this possibility with a related energy cost estimated to be ~5 kcal.mol(-1). The phosphatidylinositol (4,5) bisphosphate (PI(4,5)P2) head binds preferentially to the MA highly basic region as described in available NMR data, but interestingly without flipping of its 2' acyl chain into the MA protein. Moreover, MA was able to confine PI(4,5)P2 lipids all around its molecular surface after having found a stable orientation at the membrane surface. Our results suggest that this orientation is dependent on Myr anchoring and that this confinement induces a lateral segregation of PI(4,5)P2 in domains. This is consistent with a PI(4,5)P2 enrichment of the virus envelope as compared to the host cell membrane.

The IbeA protein from adherent invasive Escherichia coli is a flavoprotein sharing structural homology with FAD-dependent oxidoreductases.

Invasion of brain endothelium protein A (IbeA) is a virulence factor specific to pathogenic Escherichia coli. Originally identified in the K1 strain causing neonatal meningitis, it was more recently found in avian pathogenic Escherichia coli (APEC) and adherent invasive Escherichia coli (AIEC). In these bacteria, IbeA facilitates host cell invasion and intracellular survival, in particular, under harsh conditions like oxidative stress. Furthermore, IbeA from AIEC contributes to intramacrophage survival and replication, thus enhancing the inflammatory response within the intestine. Therefore, this factor is a promising drug target for anti-AIEC strategies in the context of Crohn's disease. Despite such an important role, the biological function of IbeA remains largely unknown. In particular, its exact nature and cellular localization, i.e., membrane-bound invasin versus cytosolic factor, are still of debate. Here, we developed an efficient protocol for recombinant expression of IbeA under native conditions and demonstrated that IbeA from AIEC is a soluble, homodimeric flavoprotein. Using mass spectrometry and tryptophan fluorescence measurements, we further showed that IbeA preferentially binds flavin adenine dinucleotide (FAD), with an affinity in the one-hundred nanomolar range and optimal binding under reducing conditions. 3D-modeling with AlphaFold revealed that IbeA shares strong structural homology with FAD-dependent oxidoreductases. Finally, we used ligand docking, mutational analyses, and molecular dynamics simulations to identify the FAD binding pocket within IbeA and characterize possible conformational changes occurring upon ligand binding. Overall, we suggest that the role of IbeA in the survival of AIEC within host cells, notably macrophages, is linked to modulation of redox processes.

Structural basis of the mycobacterial stress-response RNA polymerase auto-inhibition via oligomerization.

Self-assembly of macromolecules into higher-order symmetric structures is fundamental for the regulation of biological processes. Higher-order symmetric structure self-assembly by the gene expression machinery, such as bacterial DNA-dependent RNA polymerase (RNAP), has never been reported before. Here, we show that the stress-response σB factor from the human pathogen, Mycobacterium tuberculosis, induces the RNAP holoenzyme oligomerization into a supramolecular complex composed of eight RNAP units. Cryo-electron microscopy revealed a pseudo-symmetric structure of the RNAP octamer in which RNAP protomers are captured in an auto-inhibited state and display an open-clamp conformation. The structure shows that σB is sequestered by the RNAP flap and clamp domains. The transcriptional activator RbpA prevented octamer formation by promoting the initiation-competent RNAP conformation. Our results reveal that a non-conserved region of σ is an allosteric controller of transcription initiation and demonstrate how basal transcription factors can regulate gene expression by modulating the RNAP holoenzyme assembly and hibernation.

Second-Generation CD73 Inhibitors Based on a 4,6-Biaryl-2-thiopyridine Scaffold.

Various series of 4,6-biaryl-2-thiopyridine derivatives were synthesized and evaluated as potential ecto-5'-nucleotidase (CD73) inhibitors. Two synthetic routes were explored and the coupling of 4,6-disubstituted 3-cyano-2-chloro-pyridines with selected thiols allowed us to explore the structural diversity. Somehow divergent results were obtained in biological assays on CD73 inhibition using either the purified recombinant protein or cell-based assays, highlighting the difficulty to target protein-protein interface on proteins existing as soluble and membrane-bound forms. Among the 18 new derivatives obtained, three derivatives incorporating morpholino substituents on the 4,6-biaryl-2-thiopyridine core were shown to be able to reverse the adenosine-mediated immune suppression on human T cells. The higher blockade efficiency was observed for 2-((3-cyano-4,6-bis(4-morpholinophenyl)pyridin-2-yl)thio)-N-(isoxazol-3-yl)acetamide (with total reversion at 100 µM) and methyl 2-((3-cyano-4,6-bis(4-morpholinophenyl)pyridin-2-yl)thio)acetate (with partial reversion at 10 µM). Thus, this series of compounds illustrates a new chemotype of CD73 allosteric inhibitors.

Structural insights into the inhibition of cytosolic 5'-nucleotidase II (cN-II) by ribonucleoside 5'-monophosphate analogues.

Cytosolic 5'-nucleotidase II (cN-II) regulates the intracellular nucleotide pools within the cell by catalyzing the dephosphorylation of 6-hydroxypurine nucleoside 5'-monophosphates. Beside this physiological function, high level of cN-II expression is correlated with abnormal patient outcome when treated with cytotoxic nucleoside analogues. To identify its specific role in the resistance phenomenon observed during cancer therapy, we screened a particular class of chemical compounds, namely ribonucleoside phosphonates to predict them as potential cN-II inhibitors. These compounds incorporate a chemically and enzymatically stable phosphorus-carbon linkage instead of a regular phosphoester bond. Amongst them, six compounds were predicted as better ligands than the natural substrate of cN-II, inosine 5'-monophosphate (IMP). The study of purine and pyrimidine containing analogues and the introduction of chemical modifications within the phosphonate chain has allowed us to define general rules governing the theoretical affinity of such ligands. The binding strength of these compounds was scrutinized in silico and explained by an impressive number of van der Waals contacts, highlighting the decisive role of three cN-II residues that are Phe 157, His 209 and Tyr 210. Docking predictions were confirmed by experimental measurements of the nucleotidase activity in the presence of the three best available phosphonate analogues. These compounds were shown to induce a total inhibition of the cN-II activity at 2 mM. Altogether, this study emphasizes the importance of the non-hydrolysable phosphonate bond in the design of new competitive cN-II inhibitors and the crucial hydrophobic stacking promoted by three protein residues.

Identification of a non-canonical G3BP-binding sequence in a Mayaro virus nsP3 hypervariable domain.

Ras-GTPase-activating SH3 domain-binding-proteins 1 (G3BP1) and 2 (G3BP2) are multifunctional RNA-binding proteins involved in stress granule nucleation, previously identified as essential cofactors of Old World alphaviruses. They are recruited to viral replication complexes formed by the Chikungunya virus (CHIKV), Semliki Forest virus (SFV), and Sindbis virus (SINV) via an interaction with a duplicated FGxF motif conserved in the hypervariable domain (HVD) of virus-encoded nsP3. According to mutagenesis studies, this FGxF duplication is strictly required for G3BP binding and optimal viral growth. Contrasting with this scenario, nsP3 encoded by Mayaro virus (MAYV), an arthritogenic virus grouped with Old World alphaviruses, contains a single canonical FGxF sequence. In light of this unusual feature, we questioned MAYV nsP3/G3BPs relationships. We report that G3BP1 and G3BP2 are both required for MAYV growth in human cells and bind nsP3 protein. In infected cells, they are recruited to nsP3-containing cytosolic foci and active replication complexes. Unexpectedly, deletion of the single FGxF sequence in MAYV nsP3 did not abolish these phenotypes. Using mutagenesis and in silico modeling, we identify an upstream FGAP amino acid sequence as an additional MAYV nsP3/G3BP interaction motif required for optimal viral infectivity. These results, therefore, highlight a non-conventional G3BP binding sequence in MAYV nsP3.

Cytotoxic and Antitumoral Activity of N-(9H-purin-6-yl) Benzamide Derivatives and Related Water-soluble Prodrugs.

BACKGROUND: The development of small molecules as cancer treatments is still of both interest and importance. OBJECTIVE: Having synthesized and identified the initial cytotoxic activity of a series of chemically related N-(9H-purin-6-yl) benzamide derivatives, we continued their evaluation on cancer cell models. We also synthesized water-soluble prodrugs of the main compound and performed in vivo experiments. METHOD: We used organic chemistry to obtain compounds of interest and prodrugs. The biological evaluation included MTT assays, synergy experiments, proliferation assays by CFSE, cell cycle distribution and in vivo antitumoral activity. RESULTS: Our results show activities on cancer cell lines ranging from 3-39 µM for the best compounds, with both induction of apoptosis and decrease in cell proliferation. Two compounds evaluated in vivo showed weak antitumoral activity. In addition, the lead compound and its prodrug had a synergistic activity with the nucleoside analogue fludarabine in vitro and in vivo. CONCLUSION: Our work allowed us to gain better knowledge on the activity of N-(9H-purin-6-yl) benzamide derivatives and showed new examples of water-soluble prodrugs. More research is warranted to decipher the molecular mechanisms of the molecules.

Ligand chirality effects on the dynamics of human 3-phosphoglycerate kinase: comparison between D- and L-nucleotides.

l-nucleoside analogues are now largely used as antiviral drugs for the treatment of viral infections like HBV, HCV and HIV. However, in order to be fully active, they need to be phosphorylated by cellular or viral kinases. Human 3-phosphogycerate kinase (hPGK) was shown to catalyze the last step of activation of l-enantiomers and thus constitutes an attractive target for theoretical predictions of its phosphorylation efficiency. Molecular dynamics simulations were carried out with four different nucleotides (d-/l-ADP and d-/l-CDP) in complex with hPGK and 1,3-bisphospho-d-glycerate (bPG). The binding affinities of CDPs (both enantiomers) for hPGK were found very weak while d- and l-ADP were better substrates. Interestingly, the binding affinity of the bPG substrate was found to be lower in presence of d-ADP than l-ADP which indicates a potential antagonistic effect on one substrate to the other. A detailed analysis of the simulations unravels important dynamic conditions for efficient phosphorylation. Indeed, as previously described for the natural substrate, the hinge bending motion of the domains upon substrates binding should be more correlated and directional. Interestingly, the unforeseen finding was the larger dynamics freedoms observed for the substrates that was favored by the protein atoms flexibility around the nucleobase binding site.

Electrostatic repulsion between HIV-1 capsid proteins modulates hexamer plasticity and in vitro assembly.

Capsid protein (CA) is the major component of the human immunodeficiency virus type 1 (HIV-1) core. Three major phosphorylation sites have been identified at positions S(109), S(149) and S(178) in the amino-acid sequence of CA. Here, we investigated the possible consequences of phosphorylation at these sites on the CA hexamer organization and plasticity using in silico approaches. The biological relevance of molecular modeling was then evaluated by analyzing the in vitro assembly properties of bacterially expressed CA bearing S(109)D, S(149)D, or S(178)D substitutions that mimic constitutive phosphorylation at these sites. We found that a constitutive negative charge at position 109 or 149 impaired the capacity of mature CA to assemble in vitro. In vivo, HIV-1 mutants bearing the corresponding mutation showed dramatic alterations of core morphology. At the level of CA hexamer, S(149) phosphorylation generates inter-monomer repulsions, while phosphorylation at position 109 resulted in cleavage of important bonds required for preserving the stability of the edifice. Addition of a negative charge at position 178 allowed efficient assembly of CA into core-like structures in vitro and in vivo and significantly increased CA hexamer stability when modeled in silico. All mutant viruses studied lacked infectivity since they were unable to produce proviral DNA. Altogether our data indicate that negative charges, that mimic phosphorylation, modulate assembling capacity of CA and affect structural properties of CA hexamers and of HIV-1 cores. In the context of the assembled core, phosphorylation at these sites may be considered as an event interfering with core organization and HIV-1 replicative cycle.