Stem Rot of Pearl Millet Prevalence, Symptomatology, Disease Cycle, Disease Rating Scale and Pathogen Characterization in Pearl Millet-Klebsiella Pathosystem

The oldest and most extensively cultivated form of millet, known as pearl millet (Pennisetum glaucum (L.) R. Br. Syn. Pennisetum americanum (L.) Leeke), is raised over 312.00 lakh hectares in Asian and African countries. India is regarded as the significant hotspot for pearl millet diversity. In the Indian state of Haryana, where pearl millet is grown, a new and catastrophic bacterial disease known as stem rot of pearl millet spurred by the bacterium Klebsiella aerogenes (formerly Enterobacter) was first observed during fall 2018. The disease appears in form of small to long streaks on leaves, lesions on stem, and slimy rot appearance of stem. The associated bacterium showed close resemblance to Klebsiella aerogenes that was confirmed by a molecular evaluation based on 16S rDNA and gyrA gene nucleotide sequences. The isolates were also identified to be Klebsiella aerogenes based on biochemical assays, where Klebsiella isolates differed in D-trehalose and succinate alkalisation tests. During fall 2021-2023, the disease has spread all the pearl millet-growing districts of the state, extending up to 70% disease incidence in the affected fields. The disease is causing considering grain as well as fodder losses. The proposed scale, consisting of six levels (0-5), is developed where scores 0, 1, 2, 3, 4, and 5 have been categorized as highly resistant, resistant, moderately resistant, moderately susceptible, susceptible, and highly susceptible disease reaction, respectively. The disease cycle, survival of pathogen, and possible losses have also been studied to understand other features of the disease.

Molecular dynamics simulation of RAC1 protein and its de novo variants related to developmental disorders

Neurodevelopmental disorders (NDDs) are conceptualized as childhood disability, but it has increasingly been recognized as lifelong neurological conditions that could notably impact adult functioning and quality of life. About 1%-3% of the general population suffers from NDDs including ADHD, ASD, IDD, communication disorders, motor disorders, etc. Studies suggest that Rho GTPases are key in neuronal development, highlighting the importance of altered GTPase signaling in NDDs. RAC1, a member of the Rho GTPase family, plays a critical role in neurogenesis, migration, synapse formation, axon growth, and regulation of actin cytoskeleton dynamics. We performed 6µs all-atom molecular dynamics simulation of native RAC1 (PDB: 3TH5) and three-point mutations (C18Y, N39S, and Y64D) related to developmental disorders to understand the impact of mutations on protein stability and functional dynamics. Our analysis, which included root mean square deviation (RMSD), root mean square fluctuation (RMSF), solvent accessible surface area (SASA), radius of gyration (Rg), free energy landscape (FEL), and principal component analysis (PCA), revealed that the N39S and Y64D mutations induced significant structural changes in RAC1. These alterations primarily occurred in the functional region adjacent to switch II, a region crucial for complex conformational rearrangements during the GDP and GTP exchange cycle.Communicated by Ramaswamy H. Sarma.

Markov State Models Reconcile Conformational Plasticity of GTPase with Its Substrate Binding Event

Ras GTPase is an enzyme that catalyzes the hydrolysis of guanosine triphosphate (GTP) and plays an important role in controlling crucial cellular signaling pathways. However, this enzyme has always been believed to be undruggable due to its strong binding affinity with its native substrate GTP. To understand the potential origin of high GTPase/GTP recognition, here we reconstruct the complete process of GTP binding to Ras GTPase via building Markov state models (MSMs) using a 0.1 ms long all-atom molecular dynamics (MD) simulation. The kinetic network model, derived from the MSM, identifies multiple pathways of GTP en route to its binding pocket. While the substrate stalls onto a set of non-native metastable GTPase/GTP encounter complexes, the MSM accurately discovers the native pose of GTP at its designated catalytic site in crystallographic precision. However, the series of events exhibit signatures of conformational plasticity in which the protein remains trapped in multiple non-native conformations even when GTP has already located itself in its native binding site. The investigation demonstrates mechanistic relays pertaining to simultaneous fluctuations of switch 1 and switch 2 residues which remain most instrumental in maneuvering the GTP-binding process. Scanning of the crystallographic database reveals close resemblance between observed non-native GTP binding poses and precedent crystal structures of substrate-bound GTPase, suggesting potential roles of these binding-competent intermediates in allosteric regulation of the recognition process.

Resolving Protein Conformational Plasticity and Substrate Binding via Machine Learning

A long-standing target in elucidating the biomolecular recognition process is the identification of binding-competent conformations of the receptor protein. However, protein conformational plasticity and the stochastic nature of the recognition processes often preclude the assignment of a specific protein conformation to an individual ligand-bound pose. Here, we demonstrate that a computational framework coined as RF-TICA-MD, which integrates an ensemble decision-tree-based Random Forest (RF) machine learning (ML) technique with an unsupervised dimension reduction approach time-structured independent component analysis (TICA), provides an efficient and unambiguous solution toward resolving protein conformational plasticity and the substrate binding process. In particular, we consider multimicrosecond-long molecular dynamics (MD) simulation trajectories of a ligand recognition process in solvent-inaccessible cavities of archetypal proteins T4 lysozyme and cytochrome P450cam. We show that in a scenario in which clear correspondence between protein conformation and binding-competent macrostates could not be obtained via an unsupervised dimension reduction approach, an a priori decision-tree-based supervised classification of the simulated recognition trajectories via RF would help characterize key amino acid residue pairs of the protein that are deemed sensitive for ligand binding. A subsequent unsupervised dimensional reduction of the selected residue pairs via TICA would then delineate a conformational landscape of protein which is able to demarcate ligand-bound poses from unbound ones. The proposed RF-TICA-MD approach is shown to be data agnostic and found to be robust when using other ML-based classification methods such as XGBoost. As a promising spinoff of the protocol, the framework is found to be capable of identifying distal protein locations which would be allosterically important for ligand binding and would characterize their roles in recognition pathways. A Python implementation of a proposed ML workflow is available in GitHub https://github.com/navjeet0211/rf-tica-md.

First report of Klebsiella Leaf Streak on Sorghum Caused by Klebsiella variicola in Haryana, India

Sorghum (Sorghum bicolor [L.] Moench) is one of the top ten cereal crops in the world and is grown for fodder and seed purposes. During the fall of 2019 to 2022, a disease causing small to long streaks on leaves was observed in sorghum fields of Hisar (29° 9' 6.6996'' N, 75° 43' 16.0428'' E), Rohtak (28° 53' 43.8540'' N, 76° 36' 23.8068'' E) and Mohindergarh (28° 16' 6.0492'' N, 76° 9' 3.3552'' E) regions of Haryana between July and October. The reddish brown streaks were observed in the interveinal spaces of upper and lower leaves. The disease incidence reached 20-30% of plants in affected fields. The diseased leaf tissues were disinfected with 70% alcohol and placed in a tube with sterile water. After 30 minutes, 100 µl of the suspension was inoculated onto nutrient agar medium, incubated at 28 ± 2°C for three days, and a pure culture was obtained by restreaking on nutrient agar (Janse, 2005). The rod-shaped gram-negative bacterium with round, cream to white colonies was positive for methyl red, citrate utilization, urease activity, and glucose, lactose, sorbitol, rhamnose and sucrose fermentation tests. The genomic DNA of the bacterial suspension was extracted and 16S rDNA was amplified using universal 27F/1492R primers (Marchesi et al. 1998), resulting in tentative identification as Klebsiella sp. It was further confirmed with PCR amplification of Klebsiella specific primers (F:5'-CGCGTACTATACGCCATGAACGTA-3'; R:5'-ACCGTTGATCACTTCGGTCAGG-3') for gyrA gene (Brisse and Verhoef 2001). Discrete PCR amplicons of 1,500 (16S rDNA) and 300 bp (gyrA) were observed in a 1% (w/v) agarose gel. Forward and reverse DNA sequencing of both amplicons of the Hisar isolate (VMKV101) was carried out using a BDT v3.1 Cycle sequencing kit and consensus sequences were generated by using the program SeqMan Ultra (DNASTAR Lasergene). Sequences of the PCR products were deposited in GenBank with accession numbers MZ569433 (16S rDNA) and OP390080 (gyrA). The 16S rDNA sequence was 97.32% similar to K. variicola strain 13450 (CP026013; 1,450/1,490 bp) and the gyrA sequence had 99.66% similarity to K. variicola strain FH-1 (CP054254; 297/298 bp). A 16S RNA-based phylogenetic tree done by MEGA11 (Tamura et al. 2021) using the Maximum Likelihood method showed that strain VMKV101 clustered with K. variicola type strain F2R9. The complete bacterial genome (GCA025629215), sequenced by the Ion GeneStudio S5 system using Ion 530 chips (Thermo Fisher Scientific), was 99.03% identical by average nucleotide identity (ANI) to the type genome (CP045783) of Klebsiella variicola, with 87.8% genome coverage. For pathogenicity testing, a bacterial suspension (10 ml, 1×107 colony forming units/ml) was injected into the whole whorl after mechanical injury on 15-20 days old seedlings of the susceptible genotype HC-171, then plants were incubated at 35 ± 2°C, >80% relative humidity. Control plants were injected with sterile distilled water. Initial symptoms were observed on leaves of inoculated plants after 5 to 7 days as narrow, small longitudinal reddish brown streaks. As the disease progressed, the streaks on the leaf blade increased in number and size maintaining the reddish brown color. These streaks had slightly wavy margins and were surrounded by bright yellow halos. After 15 to 20 days, the streaks were 0.5 to 2.0 mm wide and 1.0 to 5.0 cm long, occasionally up to 10.0 cm long on both side of the leaves. Over time, neighboring streaks coalesced to form large necrotic areas. All inoculated plants exhibited identical symptoms. No symptoms were observed on control leaves. The reisolated bacterium from diseased sorghum leaves showed exactly the same morphological, biochemical and 16S RNA and gyrA molecular characteristics. To our knowledge, this is the first report of K. variicola causing a leaf streak disease on sorghum. Klebsiella species primarily cause diseases in humans and animals, but K. variicola has been found to incite banana soft rot (Fan et al. 2015) and K. aerogenes to cause stem rot in pearl millet (Malik et al. 2022). Differences of prevalence, spread and control between K. variicola and two other bacteria (Xanthomonas vasicola pv. holcicola causing Bacterial leaf streak; Paraburkholderia andropogonis causing Bacterial leaf stripe) causing leaf streak diseases on sorghum need to be determined. The identification of Klebsiella leaf streak disease lays the groundwork for future investigations into epidemiology and management of K. variicola on sorghum.

Structural and functional insights into the candidate genes associated with different developmental stages of flag leaf in bread wheat (Triticum aestivum L.)

Grain yield is one of the most important aims for combating the needs of the growing world population. The role of development and nutrient transfer in flag leaf for higher yields at the grain level is well known. It is a great challenge to properly exploit this knowledge because all the processes, starting from the emergence of the flag leaf to the grain filling stages of wheat (Triticum aestivum L.), are very complex biochemical and physiological processes to address. This study was conducted with the primary goal of functionally and structurally annotating the candidate genes associated with different developmental stages of flag leaf in a comprehensive manner using a plethora of in silico tools. Flag leaf-associated genes were analyzed for their structural and functional impacts using a set of bioinformatics tools and algorithms. The results revealed the association of 17 candidate genes with different stages of flag leaf development in wheat crop. Of these 17 candidate genes, the expression analysis results revealed the upregulation of genes such as TaSRT1-5D, TaPNH1-7B, and TaNfl1-2B and the downregulation of genes such as TaNAP1-7B, TaNOL-4D, and TaOsl2-2B can be utilized for the generation of high-yielding wheat varieties. Through MD simulation and other in silico analyses, all these proteins were found to be stable. Based on the outcome of bioinformatics and molecular analysis, the identified candidate genes were found to play principal roles in the flag leaf development process and can be utilized for higher-yield wheat production.

Reconciling conformational heterogeneity and substrate recognition in cytochrome P450

Cytochrome P450, the ubiquitous metalloenzyme involved in detoxification of foreign components, has remained one of the most popular systems for substrate-recognition process. However, despite being known for its high substrate specificity, the mechanistic basis of substrate-binding by archetypal system cytochrome P450cam has remained at odds with the contrasting reports of multiple diverse crystallographic structures of its substrate-free form. Here, we address this issue by elucidating the probability of mutual dynamical transition to the other crystallographic pose of cytochrome P450cam and vice versa via unbiased all-atom computer simulation. A robust Markov state model, constructed using adaptively sampled 84-μs-long molecular dynamics simulation trajectories, maps the broad and heterogenous P450cam conformational landscape into five key substates. In particular, the Markov state model identifies an intermediate-assisted dynamic equilibrium between a pair of conformations of P450cam, in which the substrate-recognition sites remain "closed" and "open," respectively. However, the estimate of a significantly higher stationary population of closed conformation, coupled with faster rate of open → closed transition than its reverse process, dictates that the net conformational equilibrium would be swayed in favor of "closed" conformation. Together, the investigation quantitatively infers that although a potential substrate of cytochrome P450cam would, in principle, explore a diverse array of conformations of substrate-free protein, it would mostly encounter a "closed" or solvent-occluded conformation and hence would follow an induced-fit-based recognition process. Overall, the work reconciles multiple precedent crystallographic, spectroscopic investigations and establishes how a statistical elucidation of conformational heterogeneity in protein would provide crucial insights in the mechanism of potential substrate-recognition process.

An Appraisal of Computer Simulation Approaches in Elucidating Biomolecular Recognition Pathways

Computer simulation approaches in biomolecular recognition processes have come a long way. In this Perspective, we highlight a series of recent success stories in which computer simulations have played a remarkable role in elucidating the atomic resolution mechanism of kinetic processes of protein-ligand binding in a quantitative fashion. In particular, we show that a robust combination of unbiased simulation, harnessed by a high-fidelity computing environment, and Markov state modeling approaches has been instrumental in revealing novel protein-ligand recognition pathways in multiple systems. We also elucidate the role of recent developments in enhanced sampling approaches in providing the much-needed impetus in accelerating simulation of the ligand recognition process. We identify multiple key issues, including force fields and the sampling bottleneck, which are currently preventing the field from achieving quantitative reconstruction of experimental measurements. Finally, we suggest a possible way forward via adoption of multiscale approaches and coarse-grained simulations as next steps toward efficient elucidation of ligand binding kinetics.

Role of allosteric switches and adaptor domains in long-distance cross-talk and transient tunnel formation

Transient tunnels that assemble and disassemble to facilitate passage of unstable intermediates in enzymes containing multiple reaction centers are controlled by allosteric cues. Using the 140-kDa purine biosynthetic enzyme PurL as a model system and a combination of biochemical and x-ray crystallographic studies, we show that long-distance communication between ~25-Å distal active sites is initiated by an allosteric switch, residing in a conserved catalytic loop, adjacent to the synthetase active site. Further, combinatory experiments seeded from molecular dynamics simulations help to delineate transient states that bring out the central role of nonfunctional adaptor domains. We show that carefully orchestrated conformational changes, facilitated by interplay of dynamic interactions at the allosteric switch and adaptor-domain interface, control reactivity and concomitant formation of the ammonia tunnel. This study asserts that substrate channeling is modulated by allosteric hotspots that alter protein energy landscape, thereby allowing the protein to adopt transient conformations paramount to function.

Mapping the Substrate Recognition Pathway in Cytochrome P450

Cytochrome P450s are ubiquitous metalloenzymes involved in the metabolism and detoxification of foreign components via catalysis of the hydroxylation reactions of a vast array of organic substrates. However, the mechanism underlying the pharmaceutically critical process of substrate access to the catalytic center of cytochrome P450 is a long-standing puzzle, further complicated by the crystallographic evidence of a closed catalytic center in both substrate-free and substrate-bound cytochrome P450. Here, we address a crucial question whether the conformational heterogeneity prevalent in cytochrome P450 translates to heterogeneous pathways for substrate access to the catalytic center of these metalloenzymes. By atomistically capturing the full process of spontaneous substrate association from bulk solvent to the occluded catalytic center of an archetypal system P450cam in multi-microsecond-long continuous unbiased molecular dynamics simulations, we here demonstrate that the substrate recognition in P450cam always occurs through a single well-defined dominant pathway. The simulated final bound pose resulting from these unguided simulations is in striking resemblance with the crystallographic bound pose. Each individual binding trajectory reveals that the substrate, initially placed at random locations in bulk solvent, spontaneously lands on a single key channel on the protein-surface of P450cam and resides there for an uncharacteristically long period, before correctly identifying the occluded target-binding cavity. Surprisingly, the passage of substrate to the closed catalytic center is not accompanied by any large-scale opening in protein. Rather, the unbiased simulated trajectories (∼57 μs) and underlying Markov state model, in combination with free-energy analysis, unequivocally show that the substrate recognition process in P450cam needs a substrate-induced side-chain displacement coupled with a complex array of dynamical interconversions of multiple metastable substrate conformations. Further, the work reconciles multiple precedent experimental and theoretical observations on P450cam and establishes a comprehensive view of substrate-recognition in cytochrome P450 that only occurs via substrate-induced structural rearrangements.

Assessment and optimization of collective variables for protein conformational landscape: GB1 β-hairpin as a case study

Collective variables (CVs), when chosen judiciously, can play an important role in recognizing rate-limiting processes and rare events in any biomolecular systems. However, high dimensionality and inherent complexities associated with such biochemical systems render the identification of an optimal CV a challenging task, which in turn precludes the elucidation of an underlying conformational landscape in sufficient details. In this context, a relevant model system is presented by a 16-residue β-hairpin of GB1 protein. Despite being the target of numerous theoretical and computational studies for understanding the protein folding, the set of CVs optimally characterizing the conformational landscape of the β-hairpin of GB1 protein has remained elusive, resulting in a lack of consensus on its folding mechanism. Here we address this by proposing a pair of optimal CVs which can resolve the underlying free energy landscape of the GB1 hairpin quite efficiently. Expressed as a linear combination of a number of traditional CVs, the optimal CV for this system is derived by employing the recently introduced time-structured independent component analysis approach on a large number of independent unbiased simulations. By projecting the replica-exchange simulated trajectories along these pair of optimized CVs, the resulting free energy landscape of this system is able to resolve four distinct well-separated metastable states encompassing the extensive ensembles of folded, unfolded, and molten globule states. Importantly, the optimized CVs were found to be capable of automatically recovering a novel partial helical state of this protein, without needing to explicitly invoke helicity as a constituent CV. Furthermore, a quantitative sensitivity analysis of each constituent in the optimized CV provided key insights on the relative contributions of the constituent CVs in the overall free energy landscapes. Finally, the kinetic pathways connecting these metastable states, constructed using a Markov state model, provide an optimum description of the underlying folding mechanism of the peptide. Taken together, this work offers a quantitatively robust approach toward comprehensive mapping of the underlying folding landscape of a quintessential model system along its optimized CV.

Atomistic Insights into Structural Differences between E3 and E4 Isoforms of Apolipoprotein E

Among various isoforms of Apolipoprotein E (ApoE), the E4 isoform (ApoE4) is considered to be the strongest risk factor for Alzheimer's disease, whereas the E3 isoform (ApoE3) is neutral to the disease. Interestingly, the sequence of ApoE4 differs from its wild-type ApoE3 by a single amino acid C112R in the 299-amino-acid-long sequence. Hence, the puzzle remains: how a single-amino-acid difference between the ApoE3 and ApoE4 sequences can give rise to structural dissimilarities between the two isoforms, which can potentially lead to functional differences with significant pathological consequences. The major obstacle in addressing this question has been the lack of a 3D atomistic structure of ApoE4 to date. In this work, we resolve the issue by computationally modeling a plausible atomistic 3D structure of ApoE4. Our microsecond-long atomistic simulations elucidate key structural differences between monomeric ApoE3 and ApoE4, which renders ApoE4 thermodynamically less stable, less structured, and topologically less rigid compared to ApoE3. Consistent with an experimental report of the molten globule state of ApoE4, simulations identify multiple partially folded intermediates for ApoE4, which are implicated in the stronger aggregation propensity of ApoE4.

Molecular Mechanism of Nucleotide-Dependent Allosteric Regulation in AMP-Activated Protein Kinase

The AMP-activated protein kinase (AMPK), a central enzyme in the regulation of energy homeostasis, is an important drug target for type 2 diabetes, obesity, and cancer. Binding of adenosine nucleotides to the regulatory γ-subunit tightly regulates the activity of this enzyme. Though recent crystal structures of AMPK have provided important insights into the allosteric activation of AMPK, molecular details of the regulatory mechanism of AMPK activation is still elusive. Here, we have performed extensive all-atom molecular dynamics (MD) simulations and shown that the kinase domain (KD) and γ-subunit come closer resulting in a more compact heterotrimeric AMPK complex in AMP-bound state compared to the ATP-bound state. The binding of ATP at site 3 of regulatory γ-subunit allosterically inhibits AMPK by destabilizing different regulatory regions of α-subunit: the autoinhibitory domain, the linker region, and the activation loop of the kinase core. The catalytically important residues experience a change in mechanical stress, and major rearrangements in community structure derived from residue-residue interaction energy-based network are observed in KD and α-linker region upon binding of different nucleotides. Our results also highlight the role of conserved charged residues forming an ionic network near the site 3 of γ-subunit in allosteric communications.

Activation of corticotropin-releasing factor 1 receptor: insights from molecular dynamics simulations

G-protein-coupled receptors (GPCRs) constitute the largest family of membrane-bound proteins involved in translation of extracellular signals into intracellular responses. They regulate diverse physiological and pathophysiological processes, and hence, they are prime drug targets for therapeutic intervention. In spite of the recent advancements in membrane protein crystallography, limited information is available on the molecular signatures of activation of GPCRs. Although few studies have been reported for class A GPCRs, the activation mechanism of class B GPCRs remains unexplored. Corticotropin-releasing factor 1 receptor (CRF1R), a class B GPCR, is associated with various disease conditions including stress, anxiety, and irritable bowel syndrome. Here, we report the activation of CRF1R using accelerated molecular dynamics simulations of the apo receptor. The breakage of His155(2.50)-Glu209(3.50) and Glu209(3.50)-Thr316(6.42) interactions is found to be crucial in transition of the receptor to its active conformation. Compared to the inactive crystal structure, major structural rearrangements occurred in the intracellular region of the transmembrane (TM) domain upon activation: TM3 twisted away from TM2, and an opening of the G-protein binding site occurred as a result of the outward movements of TM5 and TM6 from the helical bundle. Further, an inward tilt of TM7 toward the helical core is observed at the extracellular side, in agreement with recent findings (Coin et al. Cell 2013, 155, 1258-1269), where it is proposed that this movement helps in establishing favorable interactions with peptide agonist. Moreover, different allosteric pathways in the inactive and active states are identified using the correlations in torsion angle space. The inactive state is found to be less dynamic as compared to the putative active state of the receptor. Results from the current study could present a model for class B GPCRs activation and aid in the design of CRF1R modulators against brain and metabolic disorders.

An ebselen like catalyst with enhanced GPx activity via a selenol intermediate

Benzamide ring-substituted, quinine-derived ebselen analogue is synthesized which exists in selenol form upon addition of PhSH. It catalyses oxidation of PhSH with H₂O₂ faster (10³-fold) than ebselen.