L1 Loop Research Articles

Isolation and Molecular Characterization of Fowl Adenoviruses Associated with Inclusion Body Hepatitis-hydropericardium Syndrome in Broiler Chickens in Egypt

Avian adenoviruses are an extremely diversified group of pathogens that recently triggering a variety of problems for poultry production. In particular, Inclusion Body Hepatitis-Hydropericardium Syndrome (IBH-HPS), which has been observed in broiler birds from 3 to 6 weeks of age and is associated with anemia, hemorrhagic disorders, hydropericardium, and high mortality. The disease has been reported worldwide, and recently it was reported in many Egyptian farms, causing severe economic losses. Therefore, the current study aimed to isolate, and genetically type the most common Adenovirus serotypes associated with this syndrome in Egyptian farms. A total of 50 broiler chicken farms (3-6 weeks old) located in different Egyptian governorates were examined. Macroscopically, the diseased flock revealed hydropericardium, enlarged friable livers with ecchymotic hemorrhages, and varying mortality rates (1 to 7.5%). Histopathologically, severe diffuse necrotizing enteritis, hepatitis, pericarditis, and diffuse lymphoid depletion of the spleen were the most prominent lesions. Liver tissues and cloacal swabs were collected from all examined flocks for FAdVs detection by conventional polymerase chain reaction (PCR) targeting the L1 loop in the hexon gene. The PCR products were sequenced for typing of the detected viruses. It was found that 10 out of 50 flocks examined were PCR positive for FAdVs (20%). Phylogenetic analysis of the sequenced genes revealed that the obtained viruses clustered with reference strains belonging to FAdV type D and E serotype 2, 11, and 8a respectively. The isolation of both FAdV type D and FAdV type E were carried out on a primary cell culture chicken embryo liver cell (CEL) and the presence of these viruses was confirmed by PCR after the appearance of cytopathic effect (CPE). From this study, it could be concluded that both FAdVs types D and E are the most common adenoviruses circulating in poultry farms suffering from hydropericardium and inclusion body hepatitis.

Conformational Insights into the Control of CNF1 Toxin Activity by Peptidyl-Prolyl Isomerization: A Molecular Dynamics Perspective.

The cytotoxic necrotizing factor 1 (CNF1) toxin from uropathogenic Escherichia coli constitutively activates Rho GTPases by catalyzing the deamidation of a critical glutamine residue located in the switch II (SWII). In crystallographic structures of the CNF1 catalytic domain (CNF1CD), surface-exposed P768 and P968 peptidyl-prolyl imide bonds (X-Pro) adopt an unusual cis conformation. Here, we show that mutation of each proline residue into glycine abrogates CNF1CD in vitro deamidase activity, while mutant forms of CNF1 remain functional on RhoA in cells. Using molecular dynamics simulations coupled to protein-peptide docking, we highlight the long-distance impact of peptidyl-prolyl cis-trans isomerization on the network of interactions between the loops bordering the entrance of the catalytic cleft. The energetically favorable isomerization of P768 compared with P968, induces an enlargement of loop L1 that fosters the invasion of CNF1CD catalytic cleft by a peptide encompassing SWII of RhoA. The connection of the P968 cis isomer to the catalytic cysteine C866 via a ladder of stacking interactions is alleviated along the cis-trans isomerization. Finally, the cis-trans conversion of P768 favors a switch of the thiol side chain of C866 from a resting to an active orientation. The long-distance impact of peptidyl-prolyl cis-trans isomerizations is expected to have implications for target modification.

The SARS-CoV-2 spike protein is vulnerable to moderate electric fields

Most of the ongoing projects aimed at the development of specific therapies and vaccines against COVID-19 use the SARS-CoV-2 spike (S) protein as the main target. The binding of the spike protein with the ACE2 receptor (ACE2) of the host cell constitutes the first and key step for virus entry. During this process, the receptor binding domain (RBD) of the S protein plays an essential role, since it contains the receptor binding motif (RBM), responsible for the docking to the receptor. So far, mostly biochemical methods are being tested in order to prevent binding of the virus to ACE2. Here we show, with the help of atomistic simulations, that external electric fields of easily achievable and moderate strengths can dramatically destabilise the S protein, inducing long-lasting structural damage. One striking field-induced conformational change occurs at the level of the recognition loop L3 of the RBD where two parallel beta sheets, believed to be responsible for a high affinity to ACE2, undergo a change into an unstructured coil, which exhibits almost no binding possibilities to the ACE2 receptor. We also show that these severe structural changes upon electric-field application also occur in the mutant RBDs corresponding to the variants of concern (VOC) B.1.1.7 (UK), B.1.351 (South Africa) and P.1 (Brazil). Remarkably, while the structural flexibility of S allows the virus to improve its probability of entering the cell, it is also the origin of the surprising vulnerability of S upon application of electric fields of strengths at least two orders of magnitude smaller than those required for damaging most proteins. Our findings suggest the existence of a clean physical method to weaken the SARS-CoV-2 virus without further biochemical processing. Moreover, the effect could be used for infection prevention purposes and also to develop technologies for in-vitro structural manipulation of S. Since the method is largely unspecific, it can be suitable for application to other mutations in S, to other proteins of SARS-CoV-2 and in general to membrane proteins of other virus types.

Clinico-genomic profiling and clonal dynamic modeling of TP53-aberrant myelodysplastic syndrome and acute myeloid leukemia

TP53-aberrant myelodysplastic syndrome and acute myeloid leukemia have dismal outcomes. Here, we define the clinico-genomic landscape of TP53 disruptions in 40 patients and employ clonal dynamic modeling to map the mutational hierarchy against clinical outcomes. Most TP53 mutations (45.2%) localized to the L3 loop or LSH motif of the DNA-binding domain. TP53 disruptions had high co-occurrence with mutations in epigenetic regulators, spliceosome machinery, and cohesin complex and low co-occurrence with mutations in proliferative signaling genes. Ancestral and descendant TP53 mutations constituted measurable residual disease and fueled relapse. High mutant TP53 gene dosage predicted low durability of remission. The median overall survival (OS) was 280 days. Hypomethylating agent-based therapy served as an effective bridge to transplant, leading to improved median OS compared to patients who did not receive a transplant (14.7 vs. 5.1 months). OS was independent of the genomic location of TP53 disruption, which has implications for rational therapeutic design.

Structure of a Mutant Form of Translation Regulator Hfq with the Extended Loop L4

Structure of a Mutant Form of Translation Regulator Hfq with the Extended Loop L4

Hidden electrostatic energy contributions define dynamic allosteric communications within p53 during molecular recognition

Molecular recognition is fundamental to transcription regulation. As a transcription factor, the tumor suppressor p53 has to recognize either specific DNA sequences or repressor protein partners. However, the molecular mechanism underlying the p53 conformational switch from the DNA-bound to repressor-bound states is not fully characterized. The highly charged nature of these interacting molecules prompted us to explore the nonbonded energy contributions behind molecular recognition of either a DNA or the repressor protein iASPP by p53 DNA binding domain (p53DBD), using molecular dynamics simulation followed by rigorous analyses of energy terms. Our results illuminate the allosteric pathway by which iASPP binding to p53 diminishes binding affinity between p53 and DNA. Even though the p53DBD uses a common framework of residues for recognizing both DNA and iASPP, a comparison of the electrostatics in the two p53DBD complexes revealed significant differences in residue-wise contributions to the electrostatic energy. We found that an electrostatic allosteric communication path exists in the presence of both substrates. It consists of evolutionarily conserved residues, from residue K120 of the binding loop L1 to a distal residue R213 of p53DBD. K120 is near the DNA in the p53DBD-DNA complex, whereas iASPP binding moves it away from its DNA binding position in the p53DBD-iASPP complex. The “energy hubs” (the residues show a higher degree of connectivity with other residues in the electrostatic networks) determined from the electrostatic network analysis established that this conformational change in K120 completely rewires the electrostatic network from K120 to R213, thereby impeding DNA binding. Furthermore, we found shifting populations of hydrogen bonds and salt bridges reduce pairwise electrostatic energies within p53DBD in its DNA-bound state.

The two redox states of the human NEET proteins\u2019 [2Fe\u20132S] clusters

The NEET proteins constitute a unique class of [2Fe–2S] proteins. The metal ions bind to three cysteines and one histidine. The proteins’ clusters exist in two redox states; the oxidized protein (containing two FeIII ions) can transfer the cluster to apo-acceptor protein(s), while the reduced form (containing one ferrous ion) remains bound to the protein frame. Here, we perform in silico and in vitro studies on human NEET proteins in both reduced and oxidized forms. Quantum chemical calculations on all available human NEET proteins structures suggest that reducing the cluster weakens the Fe–NHis and Fe–SCys bonds, similar to what is seen in other Fe–S proteins (e.g., ferredoxin and Rieske protein). We further show that the extra electron in the [2Fe–2S]+ clusters of one of the NEET proteins (mNT) is localized on the His-bound iron ion, consistently with our previous spectroscopic studies. Kinetic measurements demonstrate that the mNT [2Fe–2S]+ is released only by an increase in temperature. Thus, the reduced state of human NEET proteins [2Fe–2S] cluster is kinetically inert. This previously unrecognized kinetic inertness of the reduced state, along with the reactivity of the oxidized state, is unique across all [2Fe–2S] proteins. Finally, using a coevolutionary analysis, along with molecular dynamics simulations, we provide insight on the observed allostery between the loop L2 and the cluster region. Specifically, we show that W75, R76, K78, K79, F82 and G85 in the latter region share similar allosteric characteristics in both redox states.Graphic abstract

Fowl adenovirus strains 1/A and 11/D isolated from birds with reovirus infection.

Inclusion body hepatitis (IBH) is, in some cases, a fatal disease affecting fowl by adenovirus strains which are subdivided into 5 species (A-E). In the current study, we investigated sequences from the Loop L1 region of the hexon gene of sequences of adenovirus field stains 1/A and 11/D isolated from a poultry flock co-infected with IBH and avian reoviruses ARVs. In early 2021, an epidemiologic survey highlighted the coinfection adenoviruses with other viruses (orthoreovirus infection) as being particularly deleterious within the poultry industry. Here, we investigated the Loop L1 HVR1-4 region of the hexon gene with relative synonymous codon usage (RSCU) designation and RSCU inclusive of all the mutations. These are the first results that have been presented on fowl adenovirus species A and D with simultaneous reovirus infection in 38-days old broiler chickens in Poland.

Cryo-EM structure of the nucleosome core particle containing Giardia lamblia histones.

Giardia lamblia is a pathogenic unicellular eukaryotic parasite that causes giardiasis. Its genome encodes the canonical histones H2A, H2B, H3, and H4, which share low amino acid sequence identity with their human orthologues. We determined the structure of the G. lamblia nucleosome core particle (NCP) at 3.6 Å resolution by cryo-electron microscopy. G. lamblia histones form a characteristic NCP, in which the visible 125 base-pair region of the DNA is wrapped in a left-handed supercoil. The acidic patch on the G. lamblia octamer is deeper, due to an insertion extending the H2B α1 helix and L1 loop, and thus cannot bind the LANA acidic patch binding peptide. The DNA and histone regions near the DNA entry-exit sites could not be assigned, suggesting that these regions are asymmetrically flexible in the G. lamblia NCP. Characterization by thermal unfolding in solution revealed that both the H2A–H2B and DNA association with the G. lamblia H3–H4 were weaker than those for human H3–H4. These results demonstrate the uniformity of the histone octamer as the organizing platform for eukaryotic chromatin, but also illustrate the unrecognized capability for large scale sequence variations that enable the adaptability of histone octamer surfaces and confer internal stability.

Molecular basis for the selective recognition and ubiquitination of centromeric histone H3 by yeast E3 ligase Psh1

Centromeres are chromosomal loci marked by histone variant CenH3 (centromeric histone H3) and essential for genomic stability and cell division. The budding yeast E3 ubiquitin ligase Psh1 selectively recognizes the yeast CenH3 (Cse4) for ubiquitination and controls the cellular level of Cse4 for proteolysis, but the underlying mechanism remains largely unknown. Here, we show that Psh1 uses a Cse4-binding domain (CBD, residues 1–211) to interact with Cse4-H4 instead of H3-H4, yielding a dissociation constant (Kd) of 27 nM. Psh1 recognizes Cse4-specific residues in the L1 loop and α2 helix to ensure Cse4 binding and ubiquitination. We map the Psh1-binding region of Cse4-H4 and identify a wide range of Cse4-specific residues required for the Psh1-mediated Cse4 recognition and ubiquitination. Further analyses reveal that histone chaperone Scm3 can impair Cse4 ubiquitination by abrogating Psh1-Cse4 binding. Together, our study reveals a novel Cse4-binding mode distinct from those of known CenH3 chaperones and elucidates the mechanism by which Scm3 competes with Psh1 for Cse4 binding.

Probing the physiological roles of the extracellular loops of chitoporin from Vibrio campbellii

VhChiP, a sugar-specific porin found on the outer membrane of Vibrio campbellii, is responsible for the transport of chitooligosaccharides, allowing the bacterium to thrive in aquatic environments using chitin as a nutrient. We previously showed that VhChiP is composed of three identical subunits, each containing a 16-stranded β-barrel connected by eight extracellular loops and eight short periplasmic turns. This study is focused on the specific roles of three prominent extracellular loops of VhChiP-L2, L3, and L8. The deletion of L2 completely disrupted the L2-L2 interactions, thus destabilizing the protein trimers as well as the integrity of the secondary structure. The deletion of L3 caused a drastic loss in the binding affinity for sugar substrates because of the absence of a cluster of key amino acid residues that form the affinity sites. The removal of L8 induced pronounced gating, which is highly responsive to elevated potentials. Our data provide further information on the important roles of the three prominent loops of VhChiP: loop L2 maintains the trimeric structure and the integrity of secondary structure, loop L3 controls the binding affinity for sugar substrates, and loop L8 retains the stably open state of the channel.

Protein Dynamics Influence the Enzymatic Activity of Phospholipase A/Acyltransferases 3 and 4.

Phospholipase A/acyltransferase 3 (PLAAT3) and PLAAT4 are enzymes involved in the synthesis of bioactive lipids. Despite sequential and structural similarities, the two enzymes differ in activity and specificity. The relation between the activity and dynamics of the N-terminal domains of PLAAT3 and PLAAT4 was studied. PLAAT3 has a much higher melting temperature and exhibits less nanosecond and millisecond dynamics in the active site, in particular in loop L2(B6), as shown by NMR spectroscopy and molecular dynamics calculations. Swapping the L2(B6) loops between the two PLAAT enzymes results in strongly increased phospholipase activity in PLAAT3 but no reduction in PLAAT4 activity, indicating that this loop contributes to the low activity of PLAAT3. The results show that, despite structural similarity, protein dynamics differ substantially between the PLAAT variants, which can help to explain the activity and specificity differences.

Allosteric Site on SHIP2 Identified Through Fluorescent Ligand Screening and Crystallography: A Potential New Target for Intervention

Src homology 2 domain-containing inositol phosphate phosphatase 2 (SHIP2) is one of the 10 human inositol phosphate 5-phosphatases. One of its physiological functions is dephosphorylation of phosphatidylinositol 3,4,5-trisphosphate, PtdIns(3,4,5)P3. It is therefore a therapeutic target for pathophysiologies dependent on PtdIns(3,4,5)P3 and PtdIns(3,4)P2. Therapeutic interventions are limited by the dearth of crystallographic data describing ligand/inhibitor binding. An active site-directed fluorescent probe facilitated screening of compound libraries for SHIP2 ligands. With two additional orthogonal assays, several ligands including galloflavin were identified as low micromolar Ki inhibitors. One ligand, an oxo-linked ethylene-bridged dimer of benzene 1,2,4-trisphosphate, was shown to be an uncompetitive inhibitor that binds to a regulatory site on the catalytic domain. We posit that binding of ligands to this site restrains L4 loop motions that are key to interdomain communications that accompany high catalytic activity with phosphoinositide substrate. This site may, therefore, be a future druggable target for medicinal chemistry.

Trichoderma reesei Rad51 tolerates mismatches in hybrid meiosis with diverse genome sequences

Most eukaryotes possess two RecA-like recombinases (ubiquitous Rad51 and meiosis-specific Dmc1) to promote interhomolog recombination during meiosis. However, some eukaryotes have lost Dmc1. Given that mammalian and yeast Saccharomyces cerevisiae (Sc) Dmc1 have been shown to stabilize recombination intermediates containing mismatches better than Rad51, we used the Pezizomycotina filamentous fungus Trichoderma reesei to address if and how Rad51-only eukaryotes conduct interhomolog recombination in zygotes with high sequence heterogeneity. We applied multidisciplinary approaches (next- and third-generation sequencing technology, genetics, cytology, bioinformatics, biochemistry, and single-molecule biophysics) to show that T. reesei Rad51 (TrRad51) is indispensable for interhomolog recombination during meiosis and, like ScDmc1, TrRad51 possesses better mismatch tolerance than ScRad51 during homologous recombination. Our results also indicate that the ancestral TrRad51 evolved to acquire ScDmc1-like properties by creating multiple structural variations, including via amino acid residues in the L1 and L2 DNA-binding loops.

Structural determinants of the IRF4/DNA homodimeric complex.

Interferon regulatory factor 4 (IRF4) is a key transcription factor (TF) in the regulation of immune cells, including B and T cells. It acts by binding DNA as both a homodimer and, in conjunction with other TFs, as a heterodimer. The choice of homo and heterodimeric/ DNA interactions is a critical aspect in the control of the transcriptional program and cell fate outcome. To characterize the nature of this interaction in the homodimeric complex, we have determined the crystal structure of the IRF4/ISRE homodimeric complex. We show that the complex formation is aided by a substantial DNA deformation with co-operative binding achieved exclusively through protein–DNA contact. This markedly contrasts with the heterodimeric form where DNA bound IRF4 is shown to physically interact with PU.1 TF to engage EICE1. We also show that the hotspot residues (Arg98, Cys99 and Asn102) contact both consensus and non-consensus sequences with the L1 loop exhibiting marked flexibility. Additionally, we identified that IRF4L116R, a mutant associated with chronic lymphocytic leukemia, binds more robustly to DNA thereby providing a rationale for the observed gain of function. Together, we demonstrate key structural differences between IRF4 homo and heterodimeric complexes, thereby providing molecular insights into IRF4-mediated transcriptional regulation.

Lipid binding by the N-terminal motif mediates plasma membrane localization of Bordetella effector protein BteA

The respiratory pathogens Bordetella pertussis and Bordetella bronchiseptica employ a type III secretion system (T3SS) to inject a 69-kDa BteA effector protein into host cells. This effector is known to contain two functional domains, including an N-terminal lipid raft targeting (LRT) domain and a cytotoxic C-terminal domain that induces nonapoptotic and caspase-1–independent host cell death. However, the exact molecular mechanisms underlying the interaction of BteA with plasma membrane (PM) as well as its cytotoxic activity in the course of Bordetella infections remain poorly understood. Using a protein–lipid overlay assay and surface plasmon resonance, we show here that the recombinant LRT domain binds negatively charged membrane phospholipids. Specifically, we determined that the dissociation constants of the LRT domain–binding liposomes containing phosphatidylinositol 4,5-bisphosphate, phosphatidic acid, and phosphatidylserine were ∼450 nM, ∼490 nM, and ∼1.2 μM, respectively. Both phosphatidylserine and phosphatidylinositol 4,5-bisphosphate were required to target the LRT domain and/or full-length BteA to the PM of yeast cells. The membrane association further involved electrostatic and hydrophobic interactions of LRT and depended on a leucine residue in the L1 loop between the first two helices of the four-helix bundle. Importantly, charge-reversal substitutions within the L1 region disrupted PM localization of the BteA effector without hampering its cytotoxic activity during B. bronchiseptica infection of HeLa cells. The LRT-mediated targeting of BteA to the cytosolic leaflet of the PM of host cells is, therefore, dispensable for effector cytotoxicity.

Terahertz physical mechanism of arsenic (III) restoring p53 mutant protein activity

Recovering the protein activity of p53 mutants through small molecule ligand binding (eg. arsenic) is an important strategy for tumor suppressor therapy. However, the mechanistic basis on the changes of collective dynamics and their roles of p53 protein in functional recovery process has not been fully elucidated. Herein, the normal mode calculations based on all-atom elastic network model are employed to characterize the terahertz low frequency motions of core DNA-binding domain (p53C) which is essential for p53 protein activities in transcriptional transactivation. We find that the lowest-frequency collective vibration mode of the p53C mutant is effectively restored by the binding of arsenic (III) ligand. In R249S mutant, the L1 loop is stabilized through restricting the swing-out movement. The results obtained from atomic backbone fluctuations suggest that the arsenic binding can significantly improve the L1 loop and L2 loop fluctuations. The statistical analysis of low frequency vibration mode reflects that the arsenic-bound R249S mutant has an apparent recovery of frequency shift in the terahertz range. The residue-residue motion correlation also suggests that structural components binding to arsenic are dynamically coupled. In the H2 helix with arsenic-binding residues, the motions of C124, C135, M133 and C141, are correlated with the arsenic recovery. These results provide the terahertz biophysical mechanism for the recovery effect of arsenic (III) on the p53 protein activity and new evidence for the coupling of the low-frequency vibration characteristics of protein structures with its function, thus giving a new physical insight into the p53 related cancer therapies.

Markov state models and NMR uncover an overlooked allosteric loop in p53.

The tumor suppressor p53 is the most frequently mutated gene in human cancer, and thus reactivation of mutated p53 is a promising avenue for cancer therapy. Analysis of wildtype p53 and the Y220C cancer mutant long-timescale molecular dynamics simulations with Markov state models and validation by NMR relaxation studies has uncovered the involvement of loop L6 in the slowest motions of the protein. Due to its distant location from the DNA-binding surface, the conformational dynamics of this loop has so far remained largely unexplored. We observe mutation-induced stabilization of alternate L6 conformations, distinct from all experimentally-determined structures, in which the loop is both extended and located further away from the DNA-interacting surface. Additionally, the effect of the L6-adjacent Y220C mutation on the conformational landscape of the functionally-important loop L1 suggests an allosteric role to this dynamic loop and the inactivation mechanism of the mutation. Finally, the simulations reveal a novel Y220C cryptic pocket that can be targeted for p53 rescue efforts. Our approach exemplifies the power of the MSM methodology for uncovering intrinsic dynamic and kinetic differences among distinct protein ensembles, such as for the investigation of mutation effects on protein function.

A recent perspective on fiber and hexon genes proteins analyses of fowl adenovirus toward virus infectivity-A review.

Fowl adenovirus (FAdV) is a double-stranded DNA virus with a non-enveloped structure comprising three major proteins known as hexon, penton, and fiber. Molecular analysis which emphasizes on hexon and fiber proteins is currently the major focus of curiosity for FAdV antigenicity and pathogenicity. Recently, disease outbreaks associated with FAdV infections such as inclusion body hepatitis, hepatitis hydropericardium syndrome, and gizzard erosion, were commonly reported and continue to increase worldwide. Studies on the virulence gene of the virus were intensively conducted to provide a better understanding on the role of these major capsid proteins in the development of a safe and effective vaccine against the disease in the poultry industry. This paper highlights the variations of the fiber and hexon genes, their importance in genotypes and serotypes differentiation, and infectivity between FAdV strains. It appears that the L1 loop of hexon and the knob of fiber genes are the infectivity markers for FAdV infection. The fiber-2 protein plays a major role in FAdV pathogenicity than the hexon protein, while the fiber-1 protein is important for viral replication and assembly, regardless of virulence capability instead of infectivity. The hexon protein plays a major role in virus infectivity and tissue tropism. These findings could further enhance the knowledge of FAdV strains’ classification and evolution, diagnosis, and strategies to prevent and control FAdV infection and outbreaks in chicken farms.

Mechanistic insights into the effect of phosphorylation on Ras conformational dynamics and its interactions with cell signaling proteins

Ras undergoes interconversion between the active GTP-bound state and the inactive GDP-bound state. This GTPase cycle, which controls the activities of Ras, is accelerated by Ras GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (SOS). Oncogenic Ras mutations could affect the GTPase cycle and impair Ras functions. Additionally, Src-induced K-Ras Y32/64 dual phosphorylation has been reported to disrupt GTPase cycle and hinder Ras downstream signaling. However, the underlying mechanisms remain unclear. To address this, we performed molecular dynamics simulations (~30 μs in total) on unphosphorylated and phosphorylated K-Ras4B in GTP- and GDP-bound states, and on their complexes with GTPase cycle regulators (GAP and SOS) and the effector protein Raf. We found that K-Ras4B dual phosphorylation mainly alters the conformation at the nucleotide binding site and creates disorder at the catalytic site, resulting in the enlargement of GDP binding pocket and the retard of Ras-GTP intrinsic hydrolysis. We observed phosphorylation-induced shift in the distribution of Ras-GTP inactive-active sub-states and recognized potential druggable pockets in the phosphorylated Ras-GTP. Moreover, decreased catalytic competence or signal delivery abilities due to reduced binding affinities and/or distorted catalytic conformations of GAP, SOS and Raf were observed. In addition, the allosteric pathway from Ras/Raf interface to the distal Raf L4 loop was compromised by Ras phosphorylation. These results reveal the mechanisms by which phosphorylation influences the intrinsic or GAP/SOS catalyzed transformations between GTP- and GDP-bound states of Ras and its signal transduction to Raf. Our findings project Ras phosphorylation as a target for cancer drug discovery.