Stüve-Wiedemann syndrome (SWS), a rare autosomal recessive disorder, characterized by diminutive size, curvature of the elongated bones, bent fingers, episodes of heightened body temperature, respiratory distress or periods of breath-holding, and challenges with feeding, especially causes fatality in infants. SWS is an outcome of potential missense mutations in the leukemia inhibitory factor receptor gene reflected as numerous amino acid mutations at protein level. Employing in silico tools and techniques like mutational screening with Pred_MutHTP, I-Mutant2.0, PANTHER.db, PolyPhen, to classify mutations as deleterious/destabilizing, in conjunction with experimental data analysis, P136A and S279P emerged as \'effect\'-causing mutations. Pre-existing knowledge suggests, SWS progression is effectuated conformationally altered and dysfunctional LIFR, unable to bind to LIF and further form the LIF/LIFR/gp130 signalling complex. To gain functional insights into the effect of the said mutations on the wild type protein, an all-atom, explicit, solvent molecular dynamics simulation was performed following docking approaches. Consequently, referring to the RMSD, RMSF, protein dynamic network analysis, energy landscape plots and domain motion analysis, it was revealed that unbound LIFR_WT was more prone to LIF binding as usual whereas the mutants exhibited considerable domain closure to inhibit LIF binding. We conducted binding affinity analysis via MM/GBSA and dissociation constant estimation after LIFR-LIF docking and found the WT_complex to be more stable and compact as a whole when compared to the flexible mutant complexes thus being associated with SWS. Our study offers a route for understanding molecular level implications upon LIFR mutations which opens an avenue for therapeutic interventions.

The platform chemical 3-hydroxypropionic acid is used to synthesize various valuable materials, including bioplastics. Bifunctional malonyl-CoA reductase is a key enzyme in 3-hydroxypropionic acid biosynthesis as it catalyzes the two-step reduction of malonyl-CoA to malonate semialdehyde to 3-hydroxypropionic acid. Here, we report the cryo-EM structure of a full-length malonyl-CoA reductase protein from Chloroflexus aurantiacus (CaMCRFull). The EM model of CaMCRFull reveals a tandem helix architecture comprising an N-terminal (CaMCRND) and a C-terminal (CaMCRCD) domain. The CaMCRFull model also revealed that the enzyme undergoes a dynamic domain movement between CaMCRND and CaMCRCD due to the presence of a flexible linker between these two domains. Increasing the flexibility and extension of the linker resulted in a twofold increase in enzyme activity, indicating that for CaMCR, domain movement is crucial for high enzyme activity. We also describe the structural features of CaMCRND and CaMCRCD. This study reveals the protein structures underlying the molecular mechanism of CaMCRFull and thereby provides valuable information for future enzyme engineering to improve the productivity of 3-hydroxypropionic acid.

The introduction of CRISPR/Cas9 based gene editing has greatly accelerated therapeutic genome editing. However, the off-target DNA cleavage by CRISPR/Cas9 protein hampers its clinical translation, hindering its widespread use as a programmable genome editing tool. Although Cas9 variants with better mismatch discrimination have been developed, they have significantly lower rates of on-target DNA cleavage. Here, we have compared the dynamics of a more specific naturally occurring Cas9 from Francisella novicida (FnCas9) to the most widely used, SpCas9 protein. Long-scale atomistic MD simulation of free and gRNA bound forms of both the Cas9 proteins was performed, and their domain rearrangements and binding affinity with gRNA were compared to decipher the possible reason behind the enhanced specificity of FnCas9 protein. The greater binding affinity with gRNA, high domain electrostatics, and more volatility of FnCas9 than SpCas9 may explain its increased specificity and lower tolerance for mismatches.

CRISPR/Cas system, a newly but extensively investigated genome-editing method, harbors practical solutions for various genetic problems. It relies on short guide RNAs (gRNAs) to recruit the Cas9 protein, a DNA cleaving enzyme, to its genomic target DNAs. The Cas9 enzyme exhibits some unique properties, like the ability to differentiate self vs. non-self - DNA strands using the base-pairing potential of crRNA, i.e., only CRISPR DNA is entirely complementary to the CRISPR repeat sequences at the crRNA whereas the presence of mismatches in the upstream region of the spacer permit CRISPR interference which is inhibited in case of CRISPR-DNA, allosteric regulation in its domains, and domain reorientation on sgRNA binding. Several groups have contributed their efforts in understanding the functioning of the CRISPR/Cas system, but even then, there is a lot more to explore in this area. The structural and sequence-based understanding of the whole CRISPR-associated bacterial ortholog family landscape is still ambiguous. A better understanding of the underlying energetics of the CRISPR/Cas9 system should reveal critical parameters to design better CRISPR/Cas9s.

Substrate-binding proteins (SBPs) mediate ligand translocation and have been classified into seven clusters (A-G). Although the substrate specificities of these clusters are known to some extent, their ligand-binding mechanism(s) remain(s) incompletely understood. In this study, the list of SBPs belonging to different clusters was updated (764 SBPs) compared to the previously reported study (504 SBPs). Furthermore, a new cluster referred to as cluster H was identified. Results reveal that SBPs follow different ligand-binding mechanisms. Intriguingly, the majority of the SBPs follow the \'one domain movement\' rather than the well-known \'Venus Fly-trap\' mechanism. Moreover, SBPs of a few clusters display subdomain conformational movement rather than the complete movement of the N- and C-terminal domains.

Citrate synthase (CS) catalyzes the formation of citrate and coenzyme A from acetyl-CoA and oxaloacetate. CS exists in two forms: type I and type II. We determined the citrate-bound crystal structure of type II CS from the Hymenobacter sp. PAMC 26554 bacterium (HyCS; isolated from Antarctic lichen). Citrate molecules bound to a cleft between the large and small domains of HyCS. Structural comparison of HyCS with other type II CSs revealed that type II CSs have a highly conserved flexible hinge region (residues G264-P265 in HyCS), enabling correct positioning of active site residues. Notably, the catalytic His266 residue of HyCS interacted with Trp262 in the inactive (unliganded open) state of other type II CSs, whereas the His266 residue moved to the active site via a small-domain swing motion, interacting with the bound citrate in the closed conformation of HyCS. However, type I CSs lack this tryptophan residue and face-to-edge interactions. Thus, type II CSs might have a unique domain-motion control mechanism enabling a tight allosteric regulation. An activity assay using a W262A mutant showed a Hill coefficient of 2.4; thus, the interaction between Trp262 and His266 was closely related to the positive cooperative ligand binding of type II CS.

Glutamate dehydrogenase (GDH) releases ammonia in a reversible NAD(P)+-dependent oxidative deamination of glutamate that yields 2-oxoglutarate (2OG). In current perception, GDH contributes to Glu homeostasis and plays a significant role at the junction of carbon and nitrogen assimilation pathways. GDHs are members of a superfamily of ELFV (Glu/Leu/Phe/Val) amino acid dehydrogenases and are subdivided into three subclasses, based on coenzyme specificity: NAD+-specific, NAD+/NADP+ dual-specific, and NADP+-specific. We determined in this work that the mitochondrial AtGDH1 isozyme from A. thaliana is NAD+-specific. Altogether, A. thaliana expresses three GDH isozymes (AtGDH1-3) targeted to mitochondria, of which AtGDH2 has an extra EF-hand motif and is stimulated by calcium. Our enzymatic assays of AtGDH1 established that its sensitivity to calcium is negligible. In vivo the AtGDH1-3 enzymes form homo- and heterohexamers of varied composition. We solved the crystal structure of recombinant AtGDH1 in the apo-form and in complex with NAD+ at 2.59 and 2.03 Å resolution, respectively. We demonstrate also that both in the apo form and in 1:1 complex with NAD+, it forms D 3-symmetric homohexamers. A subunit of AtGDH1 consists of domain I, which is involved in hexamer formation and substrate binding, and of domain II which binds coenzyme. Most of the subunits in our crystal structures, including those in NAD+ complex, are in open conformation, with domain II forming a large (albeit variable) angle with domain I. One of the subunits of the AtGDH1-NAD+ hexamer contains a serendipitous 2OG molecule in the active site, causing a dramatic (∼25°) closure of the domains. We provide convincing evidence that the N-terminal peptide preceding domain I is a mitochondrial targeting signal, with a predicted cleavage site for mitochondrial processing peptidase (MPP) at Leu17-Leu18 that is followed by an unexpected potassium coordination site (Ser27, Ile30). We also identified several MPD [(+/-)-2-methyl-2,4-pentanediol] binding sites with conserved sequence. Although AtGDH1 is insensitive to MPD in our assays, the observation of druggable sites opens a potential for non-competitive herbicide design.

Matrix metaloproteinase-2 (MMP-2) is an extracellular Zn2+ protease specific to type I and IV collagens. Its expression is associated with several inflammatory, degenerative, and malignant diseases. Conformational properties, domain movements, and interactions between MMP-2 and its associated metal ions were characterized using a 1.0 µs molecular dynamics simulation. Dihedral principle component analysis revealed ten families of conformations with the greatest degree of variability occurring in the link region connecting the catalytic and hemopexin domains. Dynamic cross-correlation analysis indicated domain movements corresponding to the opening and closing of the hemopexin domain in relation to the fibronectin and catalytic domains facilitated by the link region. Interaction energies were calculated using the molecular mechanics Poisson Boltzman surface area-interaction entropy (MMPBSA-IE) analysis method and revealed strong binding energies for the catalytic Zn2+ ion 1, Ca2+ ion 1, and Ca2+ ion 3 with significant conformational stability at the binding sites of Zn2+ ion 1 and Ca2+ ion 1. Ca2+ ion 2 diffuses freely away from its crystallographically defined binding site. Zn2+ ion 2 plays a minor role in conformational stability of the catalytic domain while Ca2+ ion 3 is strongly attracted to the highly electronegative sidechains of the Asp residues around the central β-sheet core of the hemopexin domain; however, the interacting residue sidechain carboxyl groups are outside of Ca2+ ion 3\'s coordination sphere.

Staphylococcus aureus is an important cause of resistant healthcare-associated infections. It has been shown that the wall teichoic acid (WTA) may be an important drug target acting on antibiotic-resistant cells. The UDP-N-acetylglucosamine 2-epimerase, MnaA, is one of the first enzymes on the pathway for the biosynthesis of the WTA. Here, detailed molecular dynamics simulations of S. aureus MnaA were used to characterize the conformational changes that occur in the presence of UDP and UDP-GlcNac and also the energetic landscape associated with these changes. Using different simulation techniques, such as ABMD and GAMD, it was possible to assess the energetic profile for the protein with and without ligands in its active site. We found that there is a dynamic energy landscape that has its minimum changed by the presence of the ligands, with a closed structure of the enzyme being more frequently observed for the bound state while the unbound enzyme favors an opened conformation. Further structural analysis indicated that positively charged amino acids associated with UDP and UDP-GlcNac interactions play a major role in the enzyme opening movement. Finally, the energy landscape profiled in this work provides important conclusions for the design of inhibitor candidates targeting S. aureus MnaA.

In eubacteria, cyclic di-GMP (c-di-GMP) signaling is involved in virulence, persistence, motility and generally orchestrates multicellular behavior in bacterial biofilms. Intracellular c-di-GMP levels are maintained by the opposing activities of diguanylate cyclases (DGCs) and cognate phosphodiesterases (PDEs). The c-di-GMP homeostasis in Mycobacterium smegmatis is supported by DcpA, a conserved, bifunctional protein with both DGC and PDE activities. DcpA is a multidomain protein whose GAF-GGDEF-EAL domains are arranged in tandem and are required for these two activities. To gain insight into how interactions among these three domains affect DcpA activity, here we studied its domain dynamics using real-time FRET. We demonstrate that substrate binding in DcpA results in domain movement that prompts a switch from an \"open\" to a \"closed\" conformation and alters its catalytic activity. We found that a single point mutation in the conserved EAL motif (E384A) results in complete loss of the PDE activity of the EAL domain and in a significant decrease in the DGC activity of the GGDEF domain. Structural analyses revealed multiple hydrophobic and aromatic residues around Cys579 that are necessary for proper DcpA folding and maintenance of the active conformation. On the basis of these observations and taking into account additional bioinformatics analysis of EAL domain-containing proteins, we identified a critical putatively conserved motif, GCXXXQGF, that plays an important role in c-di-GMP turnover. We conclude that a substrate-induced conformational switch involving movement of a loop containing a conserved motif in the bifunctional diguanylate cyclase-phosphodiesterase DcpA controls c-di-GMP turnover in M. smegmatis.