Mechanism of voltage gating in the voltage-sensing phosphatase Ci-VSP.

Conformational changes in voltage-sensing domains (VSDs) are driven by the transmembrane electric field acting on the protein charges. Yet, the overall energetics and detailed mechanism of this process are not fully understood. Here, we determined free energy and displacement charge landscapes as well as the major conformations visited during a complete functional gating cycle in the isolated VSD of the phosphatase Ci-VSP (Ci-VSD) comprising four transmembrane helices (segments S1 to S4). Molecular dynamics simulations highlight the extent of S4 movements. In addition to the crystallographically determined activated "Up" and resting "Down" states, the simulations predict two Ci-VSD conformations: a deeper resting state ("down-minus") and an extended activated ("up-plus") state. These additional conformations were experimentally probed via systematic cysteine mutagenesis with metal-ion bridges and the engineering of proton conducting mutants at hyperpolarizing voltages. The present results show that these four states are visited sequentially in a stepwise manner during voltage activation, each step translocating one arginine or the equivalent of ∼1 e0 across the membrane electric field, yielding a transfer of ∼3 e0 charges in total for the complete process.

Passive green space exposure leading to lower aggression: The mediating role of sense of control.

Green spaces, integral to natural environments, have been extensively studied for their positive impact on mental health, yet their influence on social behavior, particularly aggression, is less explored. While prior research has predominantly emphasized the effects of actively engaging with nature, the significant role of passive nature exposure-a more common daily occurrence-has often been overlooked. We conducted two studies to explore the influence of passive green space exposure on aggression and the mediating effect of the sense of control. Study 1 (N = 240) utilized a cross-sectional survey to assess the relationship between passive green space exposure, sense of control, and aggression. Study 2 (N = 260) employed a single-factor between-subjects experimental design to further explore these relationships in a controlled environment. The results from both studies indicated that passive green space exposure is negatively related to aggression, and that this relationship is partially mediated by an increased sense of control. Specifically, passive green space exposure was found to negatively predict aggression by bolstering individuals' sense of control. These findings underscore the potential of enhancing the sense of control through environmental factors like green spaces as an effective strategy to reduce aggression. This study enriches our understanding of the broader impacts of green spaces, extending beyond mental health to include social behaviors. We discussed both the theoretical and practical implications of our findings, highlighting how urban planning and environmental design can incorporate green spaces to foster community well-being and mitigate aggressive behaviors.

Snapshots and ensembles of BTK and cIAP1 protein degrader ternary complexes.

Heterobifunctional chimeric degraders are a class of ligands that recruit target proteins to E3 ubiquitin ligases to drive compound-dependent protein degradation. Advancing from initial chemical tools, protein degraders represent a mechanism of growing interest in drug discovery. Critical to the mechanism of action is the formation of a ternary complex between the target, degrader and E3 ligase to promote ubiquitination and subsequent degradation. However, limited insights into ternary complex structures exist, including a near absence of studies on one of the most widely co-opted E3s, cellular inhibitor of apoptosis 1 (cIAP1). In this work, we use a combination of biochemical, biophysical and structural studies to characterize degrader-mediated ternary complexes of Bruton's tyrosine kinase and cIAP1. Our results reveal new insights from unique ternary complex structures and show that increased ternary complex stability or rigidity need not always correlate with increased degradation efficiency.

Error analysis of a rotating-metasurface polarimeter.

Polarimeters, which measure the polarization states of light directly, are essentially desired in many areas of science and technology. In our previous work, we have constructed a polarimeter based on a rotating-metasurface, and the polarization Stokes parameters of the light were measured with the known Mueller elements of the metasurface. Here, we further perform the error analysis of the metasurface polarimeter. The errors in the measured Stokes parameters have been formulated for the errors in Mueller elements of the metasurface. This analysis can be used to evaluate and minimize the errors of the metasurface polarimeter.

MIUH Inhibits the Hippocampal Neuron Growth in Fetal Rat by Affecting the PTEN Pathway.

Mild intrauterine hypoperfusion (MIUH) can induce placental dysfunction and lead to long-term changes during the process of brain development. A better understanding of the mechanism of MIUH will help in the development of new neuroprotective strategies for the placental chamber. To better understand the mechanism of the effect of MIUH on the neural development of offspring, we constructed a model of MIUH in pregnant rats. The proliferation, apoptosis, and autophagy of hippocampal neurons in fetal rats were studied via flow cytometry, immunofluorescence staining, JC-1 staining, western blotting, and real-time polymerase chain reaction at different time points (6, 24, 48, and 72 h). The results showed that MIUH significantly inhibited the proliferation of hippocampal neurons and promoted their apoptosis and autophagy. Simultaneously, MIUH could promote PTEN expression and affect the PTEN signaling pathway. bpV, an inhibitor of PTEN, could restore the inhibition of hippocampal nerve cell growth caused by MIUH. MIUH may inhibit neuronal proliferation and promote neuronal apoptosis and autophagy by regulating the PTEN signaling pathway.

Pathological mechanism of joint destruction in haemophilic arthropathy.

Haemophilic arthropathy (HA), caused by intra-articular haemorrhage, is one of the most common complications in patients with haemophilia. Factor replacement therapy provides missing coagulation factors to prevent children with haemophilia from joint bleeding and decreases their risk for HA. However, haemophilia patients in developing countries are still suffering from HA due to insufficient replacement therapy. Symptoms such as pain and activity limitations caused by HA seriously affect the functional abilities and quality of life of patients with HA, causing a high disability rate in the haemophilia cohort. The pathological mechanism of HA is complicated because the whole pathological mainly involves hypertrophic synovitis, osteopenia, cartilage and bone destruction, and these pathological changes occur in parallel and interact with each other. Inflammation plays an important role in the whole complex pathological process, and iron, cytokines, growth factors and other factors are involved. This review summarizes the pathological mechanism of HA to provide background for clinical and basic research.

[Diagnosis of a fetus with atelosteogenesis type 2 through combined prenatal ultrasonography and whole exome sequencing].

OBJECTIVE: To explore the genetic basis for fetus with short limbs detected by prenatal ultrasonography. METHODS: Results of clinical imaging of the fetus was collected. Amniotic fluid sample was collected through amniocentesis for the extraction of fetal DNA. Whole exome sequencing was carried out to detect variants related to the clinical phenotypes. Candidate variant was verified by Sanger sequencing. RESULTS: Prenatal ultrasound showed that the fetus had short limbs but no other abnormality. Whole exome sequencing has identified that the fetus carried two heterozygous pathogenic variants c.484G>T and c.1436dupA of the SLC26A2 gene, for which its mother and father were heterozygous carriers, respectively. CONCLUSION: The fetus was diagnosed with atelosteogenesis type 2 by combined prenatal ultrasonography and whole exome sequencing, which may be attributed to the compound heterozygous variants of the SLC26A2 gene. Above findings provided evidence for the diagnosis of the fetus and genetic counseling.

Transition path theory analysis of c-Src kinase activation.

Nonreceptor tyrosine kinases of the Src family are large multidomain allosteric proteins that are crucial to cellular signaling pathways. In a previous study, we generated a Markov state model (MSM) to simulate the activation of c-Src catalytic domain, used as a prototypical tyrosine kinase. The long-time kinetics of transition predicted by the MSM was in agreement with experimental observations. In the present study, we apply the framework of transition path theory (TPT) to the previously constructed MSM to characterize the main features of the activation pathway. The analysis indicates that the activating transition, in which the activation loop first opens up followed by an inward rotation of the αC-helix, takes place via a dense set of intermediate microstates distributed within a fairly broad "transition tube" in a multidimensional conformational subspace connecting the two end-point conformations. Multiple microstates with negligible equilibrium probabilities carry a large transition flux associated with the activating transition, which explains why extensive conformational sampling is necessary to accurately determine the kinetics of activation. Our results suggest that the combination of MSM with TPT provides an effective framework to represent conformational transitions in complex biomolecular systems.

Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases.

Protein kinases are enzymes that catalyze the covalent transfer of the γ-phosphate of an adenosine triphosphate (ATP) molecule onto a tyrosine, serine, threonine, or histidine residue in the substrate and thus send a chemical signal to networks of downstream proteins. They are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Unregulated protein kinase activity is often associated with a wide range of diseases, therefore making protein kinases major therapeutic targets. A prototypical system of central interest to understand the regulation of kinase activity is provided by tyrosine kinase c-Src, which belongs to the family of Src-related non-receptor tyrosine kinases (SFKs). Although the broad picture of autoinhibition via the regulatory domains and via the phosphorylation of the C-terminal tail is well characterized from a structural point of view, a detailed mechanistic understanding at the atomic-level is lacking. Advanced computational methods based on all-atom molecular dynamics (MD) simulations are employed to advance our understanding of tyrosine kinase activation. The computational studies suggest that the isolated kinase domain (KD) is energetically most favorable in the inactive conformation when the activation loop (A-loop) of the KD is not phosphorylated. The KD makes transient visits to a catalytically competent active-like conformation. The process of bimolecular trans-autophosphorylation of the A-loop eventually locks the KD in the active state. Activating point mutations may act by slightly increasing the population of the active-like conformation, enhancing the availability of the A-loop to be phosphorylated. The Src-homology 2 (SH2) and Src-homology 3 (SH3) regulatory domains, depending upon their configuration, either promote the inactive or the active state of the kinase domain. In addition to the roles played by the SH3, SH2, and KD, the Src-homology 4-Unique domain (SH4-U) region also serves as a key moderator of substrate specificity and kinase function. Thus, a fundamental understanding of the conformational propensity of the SH4-U region and how this affects the association to the membrane surface are likely to lead to the discovery of new intermediate states and alternate strategies for inhibition of kinase activity for drug discovery. The existence of a multitude of KD conformations poses a great challenge aimed at the design of specific inhibitors. One promising computational strategy to explore the conformational flexibility of the KD is to construct Markov state models from aggregated MD data.

Explaining why Gleevec is a specific and potent inhibitor of Abl kinase.

Tyrosine kinases present attractive drug targets for specific types of cancers. Gleevec, a well-known therapeutic agent against chronic myelogenous leukemia, is an effective inhibitor of Abl tyrosine kinase. However, Gleevec fails to inhibit closely homologous tyrosine kinases, such as c-Src. Because many structural features of the binding site are conserved, the molecular determinants responsible for binding specificity are not immediately apparent. Some have attributed the difference in binding specificity of Gleevec to subtle variations in ligand-protein interactions (binding affinity control), whereas others have proposed that it is the conformation of the DFG motif, in which ligand binding is only accessible to Abl and not to c-Src (conformational selection control). To address this issue, the absolute binding free energy was computed using all-atom molecular dynamics simulations with explicit solvent. The results of the free energy simulations are in good agreement with experiments, thereby enabling a meaningful decomposition of the binding free energy to elucidate the factors controlling Gleevec's binding specificity. The latter is shown to be controlled by a conformational selection mechanism and also by differences in key van der Waals interactions responsible for the stabilization of Gleevec in the binding pocket of Abl.

Computational study of Gleevec and G6G reveals molecular determinants of kinase inhibitor selectivity.

Gleevec is a potent inhibitor of Abl tyrosine kinase but not of the highly homologous c-Src kinase. Because the ligand binds to an inactive form of the protein in which an Asp-Phe-Gly structural motif along the activation loop adopts a so-called DFG-out conformation, it was suggested that binding specificity was controlled by a "conformational selection" mechanism. In this context, the binding affinity displayed by the kinase inhibitor G6G poses an intriguing challenge. Although it possesses a chemical core very similar to that of Gleevec, G6G is a potent inhibitor of both Abl and c-Src kinases. Both inhibitors bind to the DFG-out conformation of the kinases, which seems to be in contradiction with the conformational selection mechanism. To address this issue and display the hidden thermodynamic contributions affecting the binding selectivity, molecular dynamics free energy simulations with explicit solvent molecules were carried out. Relative to Gleevec, G6G forms highly favorable van der Waals dispersive interactions upon binding to the kinases via its triazine functional group, which is considerably larger than the corresponding pyridine moiety in Gleevec. Upon binding of G6G to c-Src, these interactions offset the unfavorable free energy cost of the DFG-out conformation. When binding to Abl, however, G6G experiences an unfavorable free energy penalty due to steric clashes with the phosphate-binding loop, yielding an overall binding affinity that is similar to that of Gleevec. Such steric clashes are absent when G6G binds to c-Src, due to the extended conformation of the phosphate-binding loop.

Generalized Scalable Multiple Copy Algorithms for Molecular Dynamics Simulations in NAMD.

Computational methodologies that couple the dynamical evolution of a set of replicated copies of a system of interest offer powerful and flexible approaches to characterize complex molecular processes. Such multiple copy algorithms (MCAs) can be used to enhance sampling, compute reversible work and free energies, as well as refine transition pathways. Widely used examples of MCAs include temperature and Hamiltonian-tempering replica-exchange molecular dynamics (T-REMD and H-REMD), alchemical free energy perturbation with lambda replica-exchange (FEP/λ-REMD), umbrella sampling with Hamiltonian replica exchange (US/H-REMD), and string method with swarms-of-trajectories conformational transition pathways. Here, we report a robust and general implementation of MCAs for molecular dynamics (MD) simulations in the highly scalable program NAMD built upon the parallel programming system Charm++. Multiple concurrent NAMD instances are launched with internal partitions of Charm++ and located continuously within a single communication world. Messages between NAMD instances are passed by low-level point-to-point communication functions, which are accessible through NAMD's Tcl scripting interface. The communication-enabled Tcl scripting provides a sustainable application interface for end users to realize generalized MCAs without modifying the source code. Illustrative applications of MCAs with fine-grained inter-copy communication structure, including global lambda exchange in FEP/λ-REMD, window swapping US/H-REMD in multidimensional order parameter space, and string method with swarms-of-trajectories were carried out on IBM Blue Gene/Q to demonstrate the versatility and massive scalability of the present implementation.

Global trends and focuses of GLP-1RA in renal disease: a bibliometric analysis and visualization from 2005 to 2022.

Glucagon-like peptide 1 receptor agonist (GLP-1RA) is a new class of glucose-lowing agents with the kidney benefit effect. This paper aims at finding the current state and hotspots of the research on GLP-1RA in kidney disease by using bibliometric methodologies and visualization maps to analyze publications and provide the direction for future studies on that topic. Literature information was obtained by retrieving the WoSCC database. Then, software like Microsoft Excel, VOSviewer, and CiteSpace was used to analyze and process obtained data. Bibliometric analysis and visualization of nations, authors, organizations, journals, keywords, and references were also done by VOSviewer and CiteSpace. A total of 991 publications written by 4747 authors from 1637 organizations in 75 countries on GLP-1RA in renal disease in Web of Science Core Collection were retrieved. The number of publications and citations kept growing from 2015 to 2022. The USA, Univ Copenhagen, and Rossing Peter are the leading country, organization, and author on this topic, respectively. All literature was published in 346 journals, and DIABETES OBESITY & METABOLISM is the journal with the most contributions. Meanwhile, most references are from DIABETES CARE. "Cardiovascular outcome" is the most frequent keyword in the total publications, and the reference cited most times is "Liraglutide and Cardiovascular Outcomes in Type 2 Diabetes" by Marso SP. The topic of GLP-1RA in renal disease has attracted more and more attention all over the world. Existing studies are mainly about clinical use in patients with diabetes, and studies on the mechanism are lacking.

Relationship between ferroptosis and mitophagy in renal fibrosis: a systematic review.

Renal fibrosis, characterised by glomerulosclerosis and tubulointerstitial fibrosis, is a typical pathological alteration in the progression of chronic kidney disease (CKD) to end-stage renal disease (ESRD). However, the limited and expensive options for treating renal fibrosis place a heavy financial burden on patients and healthcare systems. Therefore, it is significant to find an effective treatment for renal fibrosis. Ferroptosis, a non-traditional form of cell death, has been found to play an important role in acute kidney injury (AKI), tumours, neurodegenerative diseases, and so on. Moreover, a growing body of research suggests that ferroptosis might be a potential target of renal fibrosis. Meanwhile, mitophagy is a type of selective autophagy that can selectively degrade damaged or dysfunctional mitochondria as a form of mitochondrial quality control, reducing the production of reactive oxygen species (ROS), the accumulation of which is the main cause of renal fibrosis. Additionally, as a receptor of mitophagy, NIX can release beclin1 to induce mitophagy, which can also bind to solute carrier family 7 member 11 (SLC7A11) to block the activity of cystine/glutamate antitransporter (system Xc-) and inhibit ferroptosis, thereby suggesting a link between mitophagy and ferroptosis. However, there have been only limited studies on the relationship among mitophagy, ferroptosis and renal fibrosis. In this paper, we review the mechanisms of mitophagy, and describe how ferroptosis and mitophagy are related to renal fibrosis in an effort to identify potential novel targets for the treatment of renal fibrosis.

Enriched environment for offspring improves learning and memory impairments induced by sevoflurane exposure during the second trimester of pregnancy.

Studies in animals indicate that sevoflurane exposure in the second trimester of pregnancy has harmful effects on the learning and memory of offspring. Whether an enriched environment can reverse the damage of sevoflurane exposure in the second trimester of pregnancy on the learning and memory of rat offspring remains unclear. In this study, rats at 14 days of pregnancy were exposed to 3.5% sevoflurane for 2 hours and their offspring were treated with an enriched environment for 20 successive days. We found that the enriched environment for offspring increased nestin and Ki67 levels in hippocampal tissue, increased hippocampal neurogenesis, inhibited glycogen synthase kinase 3ß activity, and increased the expression of cell proliferation-related ß-catenin and apoptosis-related Bcl-2, indicating that an enriched environment reduces sevoflurane-induced damage by increasing the proliferation of stem cells in the hippocampus. These findings suggest that an enriched environment can reverse the effects of sevoflurane inhaled by rats during the second trimester of pregnancy on learning and memory of offspring. This study was approved by the Animal Ethics Committee of Shengjing Hospital of China Medical University (approval No. 2018PS07K) on January 2, 2018.

Immune cell infiltration landscape and immune marker molecular typing in preeclampsia.

Preeclampsia (PE) is an important topic in obstetrics. In this study, we used weighted gene co-expression network analysis (WGCNA) to screen the key modules related to immune cell infiltration and to identify the hub genes for the molecular subtyping of PE. We first downloaded a set of PE transcriptional data (GSE75010; 157 samples: 80 PE and 77 non-PE) from the GEO database. We then analyzed the PE samples and non-PE samples for immune cell infiltration and screened cells with differences in such infiltration. Next, we downloaded the immune-related genes from an immune-related database to screen the expression profile of the immune-related genes. Then, we obtained a candidate gene set by screening the immune-related genes differentially expressed between the two groups. We used WGCNA to construct a weighted co-expression network for these candidate genes, mined co-expression modules, and then calculated the correlation between each module and immune cells with differential infiltration. We screened the modules related to infiltrating immune cells, identified the key modules' hub genes, and determined the key module genes that interacted with each other. Finally, we obtained the hub genes related to the infiltrating immune cells. We classified the preeclampsia patients by unsupervised cluster molecular typing, determined the difference of immune cell infiltration among the different PE subtypes, and calculated the expression of hub genes in these different subtypes. In conclusion, we found 41 hub genes that may be closely related to the molecular typing of PE.

Identification of Druggable Kinase Target Conformations Using Markov Model Metastable States Analysis of apo-Abl.

Kinases are important targets for drug development. However, accounting for the impact of possible structural rearrangements on the binding of kinase inhibitors is complicated by the extensive flexibility of their catalytic domain. The dynamic N-lobe contains four particular mobile structural elements: the Asp-Phe-Gly (DFG) motif, the phosphate (P) positioning loop, the activation (A) loop, and the αC helix. In our previous study [Meng et al. J. Chem. Theory Comput. 2018 14, 2721-2732], we combined various simulation techniques with Markov state modeling (MSM) to explore the free energy landscape of Abl kinase beyond conformations that are known from X-ray crystallography. Here we examine the resulting Markov model in greater detail by analyzing its metastable states. A characterization of the states in terms of their DFG state, P-loop, and αC conformations is presented and compared to existing classification schemes. Several metastable states are found to be structurally close to known crystal structures of different kinases in complex with a variety of inhibitors. These results suggest that the set of conformations accessible to tyrosine kinases may be shared within the entire family and that the conformational dynamics of one kinase in the absence of any ligand can provide meaningful information about possible target conformations for inhibitors of any member of the kinase family.

Effects of mild intrauterine hypoperfusion in the second trimester on memory and learning function in rat offspring.

Mild intrauterine hypoperfusion (MIUH) is a serious pathological event that affects the growth and development of fetuses and offspring. MIUH can lead to growth restriction, low birth weight, neurodevelopmental disorders, and other adverse clinical outcomes. To study the effects of MIUH on learning and memory function in offspring, a model of MIUH was established by placing a coil (length 2.5 mm, diameter 0.24 mm) on the uterine artery and ovarian uterine artery of Sprague-Dawley rats in the second trimester of pregnancy (day 17). Next, 120 mg/kg lithium chloride (the MIUH + Li group) or normal saline (the MIUH group) was injected intraperitoneally into these rats. In addition, 120 mg/kg lithium chloride (the Li group) or normal saline (the SHAM group) was injected intraperitoneally into pregnant rats without coil placement. The Morris water maze was used to detect changes in learning and memory ability in the offspring at 4 weeks after birth. In the MIUH group, the escape latency and journey length before reaching the platform were both increased, and the number of times that the platform was crossed and the activity time in the target quadrant within 90 seconds were both decreased compared with the SHAM group. Immunofluorescence double staining and western blot assays demonstrated that hippocampal nestin and Ki67 (both cell-proliferation-related proteins) expression was significantly downregulated in the MIUH group compared with the SHAM group. Furthermore, western blot assays were conducted to investigate changes in related signaling pathway proteins in the brains of offspring rats, and revealed that glycogen synthase kinase 3ß (GSK3ß) expression was upregulated and ß-catenin expression was downregulated in the MIUH group compared with the SHAM group. In addition, compared with the MIUH group, the expression levels of p-GSK3ß and ß-catenin were upregulated in the MIUH + Li group. These results suggest that MIUH may affect learning and memory function in rat offspring by regulating the GSK3ß signaling pathway. The experimental procedures were approved by Animal Ethics Committee of Shengjing Hospital of China Medical University (approval No. 2018PS07K) in June 2018.

Case Report: Compound Heterozygous Phosphatidylinositol-Glycan Biosynthesis Class N (PIGN) Mutations in a Chinese Fetus With Hypotonia-Seizures Syndrome 1.

Multiple congenital anomalies-hypotonia-seizures syndrome 1 (MCAHS1) caused by phosphatidylinositol-glycan biosynthesis class N (PIGN) mutations is an autosomal recessive disease involving many systems of the body, such as the urogenital, cardiovascular, gastrointestinal, and central nervous systems. Here, compound heterozygous variants NM_012327.6:c.2427-2A > G and c.963G > A in PIGN were identified in a Chinese proband with MCAHS1. The features of the MCAHS1 family proband were evaluated to understand the mechanism of the PIGN mutation leading to the occurrence of MCAHS1. Ultrasound was conducted to examine the fetus, and his clinical manifestations were evaluated. Genetic testing was performed by whole-exome sequencing and the results were verified by Sanger sequencing of the proband and his parents. Reverse transcription-polymerase chain reaction was performed, and the products were subjected to Sanger sequencing. Quantitative PCR (Q-PCR) was conducted to compare gene expression between the patient and wild-type subjects. The compound heterozygous mutation NM_012327.6:c.2427-2A > G and c.963G > A was identified by whole-exome sequencing and was confirmed by Sanger sequencing. The NM_012327.6:c.2427-2A > G mutation led to skipping of exon 26, which resulted in a low expression level of the gene, as measured by Q-PCR. These findings provided a basis for genetic counseling and reproduction guidance in this family. Phenotype-genotype correlations may be defined by an expanded array of mutations.

pH-dependent mechanism of nitric oxide release in nitrophorins 2 and 4.

Nitrophorins are NO carrier proteins that transport and release NO through a pH-dependent conformational change. They bind NO tightly in a low pH environment and release it in a higher pH environment. Experimental evidence shows that the increase in the NO dissociation equilibrium constant, K(d), is due mainly to an increase in the NO release rate. Structural and kinetic data strongly suggest that NPs control NO escape by modulating its migration from the active site to the solvent through a pH-dependent conformational change. NP2 and NP4 are two representative proteins of the family displaying a 39% overall sequence identity, and interestingly, NP2 releases NO slower than NP4. The proposal that NPs' NO release relies mainly on the NO escape rate makes NPs a very peculiar case among typical heme proteins. The connection between the pH-dependent conformational change and ligand release mechanism is not fully understood and the structural basis for the pH induced structural transition and the different NO release patterns in NPs are unresolved, yet interesting issues. In this work, we have used state of the art molecular dynamics simulations to study the NO escape process in NP2 and NP4 in both the low and high pH states. Our results show that both NPs modulate NO release by switching between a "closed" conformation in a low pH environment and an "open" conformation at higher pH. In both proteins, the change is caused by the differential protonation of a common residue Asp30 in NP4 and Asp29 in NP2, and the NO escape route is conserved. Finally, our results show that, in NP2, the conformational change to the "open" conformation is smaller than that for NP4 which results in a higher barrier for NO release.