Molecular simulation of PcCel45A protein expressed from Aspergillus nidulans to understand its structure, dynamics, and thermostability

Discovery of andrographolide hit analog as a potent cyclooxygenase-2 inhibitor through consensus MD-simulation, electrostatic potential energy simulation and ligand efficiency metrics

Cyclooxygenase-2 (COX-2) is the key enzyme responsible for the conversion of arachidonic acid to prostaglandins that display pro-inflammatory properties and thus, it is a potential target protein to develop anti-inflammatory drugs. In this study, chemical and bio-informatics approaches have been employed to find a novel potent andrographolide (AGP) analog as a COX-2 inhibitor having better pharmacological properties than aspirin and rofecoxib (controls). The full amino acid sequenced human Alpha fold (AF) COX-2 protein (604AA) was selected and validated for its accuracy against the reported COX-2 protein structures (PDB ID: 5F19, 5KIR, 5F1A, 5IKQ and 1V0X) followed by multiple sequence alignment analysis to establish the sequence conservation. The systematic virtual screening of 237 AGP analogs against AF-COX-2 protein yielded 22 lead compounds based on the binding energy score (< - 8.0 kcal/mol). These were further screened out to 7 analogs by molecular docking analysis and investigated further for ADMET prediction, ligand efficiency metrics calculations, quantum mechanical analysis, MD simulation, electrostatic potential energy (EPE) docking simulation, and MM/GBSA. In-depth analysis revealed that AGP analog A3 (3-[2-[(1R,4aR,5R,6R,8aR)-6-hydroxy-5,6,8a-trimethyl-2-methylidene-3,4,4a,5,7,8-hexahydro-1H-naphthalen-1-yl]ethylidene]-4-hydroxyoxolan-2-one) forms the most stable complex with the AF-COX-2 showing the least RMSD value (0.37 ± 0.03 nm), a good number of hydrogen bonds (protein-ligand H-bond = 11, and protein H-bond = 525), minimum EPE score (- 53.81 kcal/mol), and lowest MM-GBSA before and after simulation (- 55.37 and - 56.25 kcal/mol, respectively) value compared to other analogs and controls. Thus, we suggest that the identified A3 AGP analog could be developed as a promising plant-based anti-inflammatory drug by inhibiting COX-2.

Fungal cellulases: protein engineering and post-translational modifications

Enzymatic degradation of lignocelluloses into fermentable sugars to produce biofuels and other biomaterials is critical for environmentally sustainable development and energy resource supply. However, there are problems in enzymatic cellulose hydrolysis, such as the complex cellulase composition, low degradation efficiency, high production cost, and post-translational modifications (PTMs), all of which are closely related to specific characteristics of cellulases that remain unclear. These problems hinder the practical application of cellulases. Due to the rapid development of computer technology in recent years, computer-aided protein engineering is being widely used, which also brings new opportunities for the development of cellulases. Especially in recent years, a large number of studies have reported on the application of computer-aided protein engineering in the development of cellulases; however, these articles have not been systematically reviewed. This article focused on the aspect of protein engineering and PTMs of fungal cellulases. In this manuscript, the latest literatures and the distribution of potential sites of cellulases for engineering have been systematically summarized, which provide reference for further improvement of cellulase properties. KEY POINTS: •Rational design based on virtual mutagenesis can improve cellulase properties. •Modifying protein side chains and glycans helps obtain superior cellulases. •N-terminal glutamine-pyroglutamate conversion stabilizes fungal cellulases.

Residues 315 and 369 in HN Protein Contribute to the Thermostability of Newcastle Disease Virus

Thermostable Newcastle disease virus (NDV) vaccines have been widely used in areas where a “cold-chain” is not reliable. However, the molecular mechanism of NDV thermostability remains poorly understood. In this work, we constructed chimeric viruses by exchanging viral fusion (F) and/or hemagglutinin-neuraminidase (HN) genes between the heat-resistant strain HR09 and thermolabile strain La Sota utilizing a reverse genetic system. The results showed that only chimeras with HN derived from the thermostable virus exhibited a thermostable phenotype at 56°C. The hemagglutinin (HA) and neuraminidase (NA) activities of chimeras with HN derived from the HR09 strain were more thermostable than those containing HN from the La Sota strain. Then, we used molecular dynamics simulation at different temperatures (310 K and 330 K) to measure the HN protein of the La Sota strain. The conformation of an amino acid region (residues 315–375) was observed to fluctuate. Sequence alignment of the HN protein revealed that residues 315, 329, and 369 in the La Sota strain and thermostable strains differed. Whether the three amino acid substitutions affected viral thermostability was investigated. Three mutant viruses based on the thermolabile strain were generated by substituting one, two or three amino acids at positions 315, 369, and 329 in the HN protein. In comparison with the parental virus, the mutant viruses containing mutations S315P and I369V possessed higher thermostablity and HA titers, NA and fusion activities. Taken together, these data indicate that the HN gene of NDV is a major determinant of thermostability, and residues 315 and 369 have important effects on viral thermostability.