SimRNAweb v2.0: a web server for RNA folding simulations and 3D structure modeling, with optional restraints and enhanced analysis of folding trajectories

Research on ribonucleic acid (RNA) structures and functions benefits from easy-to-use tools for computational prediction and analyses of RNA three-dimensional (3D) structure. The SimRNAweb server version 2.0 offers an enhanced, user-friendly platform for RNA 3D structure prediction and analysis of RNA folding trajectories based on the SimRNA method. SimRNA employs a coarse-grained model, Monte Carlo sampling and statistical potentials to explore RNA conformational space, optionally guided by spatial restraints. Recognized for its accuracy in RNA 3D structure prediction in RNA-Puzzles and CASP competitions, SimRNA is particularly useful for incorporating restraints based on experimental data. The new server version introduces performance optimizations and extends user control over simulations and the processing of results. It allows the application of various hard and soft restraints, accommodating alternative structures involving canonical and noncanonical base pairs and unpaired residues, while also integrating data from chemical probing methods. Enhanced features include an improved analysis of folding trajectories, offering advanced clustering options and multiple analyses of the generated trajectories. These updates provide comprehensive tools for detailed RNA structure analysis. SimRNAweb v2.0 significantly broadens the scope of RNA modeling, emphasizing flexibility and user-defined parameter control. The web server is available at https://genesilico.pl/SimRNAweb.

RNA tertiary structure prediction in CASP15 by the GeneSilico group: Folding simulations based on statistical potentials and spatial restraints

Ribonucleic acid (RNA) molecules serve as master regulators of cells by encoding their biological function in the ribonucleotide sequence, particularly their ability to interact with other molecules. To understand how RNA molecules perform their biological tasks and to design new sequences with specific functions, it is of great benefit to be able to computationally predict how RNA folds and interacts in the cellular environment. Our workflow for computational modeling of the 3D structures of RNA and its interactions with other molecules uses a set of methods developed in our laboratory, including MeSSPredRNA for predicting canonical and non-canonical base pairs, PARNASSUS for detecting remote homology based on comparisons of sequences and secondary structures, ModeRNA for comparative modeling, the SimRNA family of programs for modeling RNA 3D structure and its complexes with other molecules, and QRNAS for model refinement. In this study, we present the results of testing this workflow in predicting RNA 3D structures in the CASP15 experiment. The overall high score of the computational models predicted by our group demonstrates the robustness of our workflow and its individual components in terms of predicting RNA 3D structures of acceptable quality that are close to the target structures. However, the variance in prediction quality is still quite high, and the results are still too far from the level of protein 3D structure predictions. This exercise led us to consider several improvements, especially to better predict and enforce stacking interactions and non-canonical base pairs.

Role of a cysteine residue in substrate entry and catalysis in MtHIBADH: Analysis by chemical modifications and site-directed mutagenesis

Despite sharing conserved substrate-binding residues, members of 3-hydroxyisobutyrate dehydrogenase (HIBADH) superfamily show remarkable differences in substrate preference. Cysteine residues were identified within a radius of 6 Å surrounding both the active site and the substrate entry site of HIBADH enzyme from Mycobacterium tuberculosis (MtHIBADH). Chemical modification with thiol-modifying reagents, pCMB and DTNB, abrogated the dehydrogenase activity of the enzyme. The loss in activity followed pseudo-first-order kinetics as a function of the concentration of pCMB. S-HIBA (substrate) binding provided partial protection, while NAD (cofactor) binding provided ~70% protection from thiol-modifying reagent. Site-directed mutagenesis of cysteine residues present in the MtHIBADH enzyme identified the indispensable role of Cys-210 residue, located at C-terminal domain, for its dehydrogenase activity. Cys-210 mutation to serine reduced the dehydrogenase activity by ~2-fold while mutation to alanine strikingly reduced the activity by ~140-fold. C210A mutation did not perturb the state of oligomerization of the enzyme but perturbed the secondary structure content. Structural analysis revealed the involvement of Cys-210 residue in inter-chain interaction with Gln-178, which acts as hydrogen bond donor and coordinates with Cys-210 and Gly-208 of the adjacent subunit. The data demonstrate a critical role of Cys-210 residue in maintaining the conformation and rigidity of loop composed of substrate-interacting residues involved in the entry of S-HIBA substrate in MtHIBADH.