Computational Analysis of Protein Tunnels and Channels

Computational investigation of cis-1,4-polyisoprene binding to the latex-clearing protein LcpK30

Latex clearing proteins (Lcps) catalyze the oxidative cleavage of the C = C bonds in cis-1,4-polyisoprene (natural rubber), producing oligomeric compounds that can be repurposed to other materials. The active catalytic site of Lcps is buried inside the protein structure, thus raising the question of how the large hydrophobic rubber chains can access the catalytic center. To improve our understanding of hydrophobic polymeric substrate binding to Lcps and subsequent catalysis, we investigated the interaction of a substrate model containing ten carbon-carbon double bonds with the structurally characterized LcpK30, using multiple computational tools. Prediction of the putative tunnels and cavities in the LcpK30 structure, using CAVER-Pymol plugin 3.0.3, fpocket and Molecular Dynamic (MD) simulations provided valuable insights on how substrate enters from the surface to the buried active site. Two dominant tunnels were discovered that provided feasible routes for substrate binding, and the presence of two hydrophobic pockets was predicted near the heme cofactor. The larger of these pockets is likely to accommodate the substrate and to determine the size distribution of the oligomers. Protein-ligand docking was carried out using GOLD software to predict the conformations and interactions of the substrate within the protein active site. Deeper insight into the protein-substrate interactions, including close-contacts, binding energies and potential cleavage sites in the cis-1,4-polyisoprene, were obtained from MD simulations. Our findings provide further justification that the protein-substrate complexation in LcpK30 is mainly driven by the hydrophobic interactions accompanied by mutual conformational changes of both molecules. Two potential binding modes were identified, with the substrate in either extended or folded conformations. Whilst binding in the extended conformation was most favorable, the folded conformation suggested a preference for cleavage of a central double bond, leading to a preference for oligomers with 5 to 6 C = C bonds. The results provide insight into further enzyme engineering studies to improve catalytic activity and diversify the substrate and product scope of Lcps.

4-(Azolyl)-Benzamidines as a Novel Chemotype for ASIC1a Inhibitors

Acid-sensing ion channels (ASICs) play a key role in the perception and response to extracellular acidification changes. These proton-gated cation channels are critical for neuronal functions, like learning and memory, fear, mechanosensation and internal adjustments like synaptic plasticity. Moreover, they play a key role in neuronal degeneration, ischemic neuronal injury, seizure termination, pain-sensing, etc. Functional ASICs are homo or heterotrimers formed with (ASIC1-ASIC3) homologous subunits. ASIC1a, a major ASIC isoform in the central nervous system (CNS), possesses an acidic pocket in the extracellular region, which is a key regulator of channel gating. Growing data suggest that ASIC1a channels are a potential therapeutic target for treating a variety of neurological disorders, including stroke, epilepsy and pain. Many studies were aimed at identifying allosteric modulators of ASIC channels. However, the regulation of ASICs remains poorly understood. Using all available crystal structures, which correspond to different functional states of ASIC1, and a molecular dynamics simulation (MD) protocol, we analyzed the process of channel inactivation. Then we applied a molecular docking procedure to predict the protein conformation suitable for the amiloride binding. To confirm the effect of its sole active blocker against the ASIC1 state transition route we studied the complex with another MD simulation run. Further experiments evaluated various compounds in the Enamine library that emerge with a detectable ASIC inhibitory activity. We performed a detailed analysis of the structural basis of ASIC1a inhibition by amiloride, using a combination of in silico approaches to visualize its interaction with the ion pore in the open state. An artificial activation (otherwise, expansion of the central pore) causes a complex modification of the channel structure, namely its transmembrane domain. The output protein conformations were used as a set of docking models, suitable for a high-throughput virtual screening of the Enamine chemical library. The outcome of the virtual screening was confirmed by electrophysiological assays with the best results shown for three hit compounds.

Mutations in STARD8 (DLC3) May Cause 46,XY Gonadal Dysgenesis

INTRODUCTION

46,XY gonadal dysgenesis is a condition that is characterised by undeveloped testes in individuals with a male karyotype. Mutations in many genes that underlie this condition have been identified; however, there are still a considerable number of patients with an unknown genetic background. Recently, a mutation in the STARD8 X-linked gene in two sisters with 46,XY gonadal dysgenesis has been reported. It was localised within the START domain, whose homologue in Drosophila is responsible for maintaining testes integrity during their development.

METHODS

We analysed the potential pathogenicity of another STARD8 mutation, p.R887C, that was identified in a patient with 46,XY asymmetric gonadal dysgenesis. For this purpose, molecular dynamics simulations were performed.

RESULTS

These simulations revealed the full rearrangement of the helix containing the p.R887C substitution upstream from the START domain, which may cause STARD8 protein dysfunction and contribute to 46,XY gonadal dysgenesis. A comparison of the phenotypes of the three described 46,XY gonadal dysgenesis patients that harbour STARD8 mutations indicated that alterations of this gene can result in a partial or complete gonadal dysgenesis phenotype.

CONCLUSION

Based on these and previous results, it is reasonable to include STARD8 in gene panels for 46,XY gonadal dysgenesis.

Analysis and Visualization of Protein Channels, Tunnels, and Pores with MOLEonline and ChannelsDB 2.0

Channels, tunnels, and pores serve as pathways for the transport of molecules and ions through protein structures, thus participating to their functions. MOLEonline ( https://mole.upol.cz ) is an interactive web-based tool with enhanced capabilities for detecting and characterizing channels, tunnels, and pores within protein structures. MOLEonline has two distinct calculation modes for analysis of channel and tunnels or transmembrane pores. This application gives researchers rich analytical insights into channel detection, structural characterization, and physicochemical properties. ChannelsDB 2.0 ( https://channelsdb2.biodata.ceitec.cz/ ) is a comprehensive database that offers information on the location, geometry, and physicochemical characteristics of tunnels and pores within macromolecular structures deposited in Protein Data Bank and AlphaFill databases. These tunnels are sourced from manual deposition from literature and automatic detection using software tools MOLE and CAVER. MOLEonline and ChannelsDB visualization is powered by the LiteMol Viewer and Mol* viewer, ensuring a user-friendly workspace. This chapter provides an overview of user applications and usage.

Structures of channelrhodopsin paralogs in peptidiscs explain their contrasting K+ and Na+ selectivities

Kalium channelrhodopsin 1 from Hyphochytrium catenoides (HcKCR1) is a light-gated channel used for optogenetic silencing of mammalian neurons. It selects K+ over Na+ in the absence of the canonical tetrameric K+ selectivity filter found universally in voltage- and ligand-gated channels. The genome of H. catenoides also encodes a highly homologous cation channelrhodopsin (HcCCR), a Na+ channel with >100-fold larger Na+ to K+ permeability ratio. Here, we use cryo-electron microscopy to determine atomic structures of these two channels embedded in peptidiscs to elucidate structural foundations of their dramatically different cation selectivity. Together with structure-guided mutagenesis, we show that K+ versus Na+ selectivity is determined at two distinct sites on the putative ion conduction pathway: in a patch of critical residues in the intracellular segment (Leu69/Phe69, Ile73/Ser73 and Asp116) and within a cluster of aromatic residues in the extracellular segment (primarily, Trp102 and Tyr222). The two filters are on the opposite sides of the photoactive site involved in channel gating.

sMolBoxes: Dataflow Model for Molecular Dynamics Exploration

We present sMolBoxes, a dataflow representation for the exploration and analysis of long molecular dynamics (MD) simulations. When MD simulations reach millions of snapshots, a frame-by-frame observation is not feasible anymore. Thus, biochemists rely to a large extent only on quantitative analysis of geometric and physico-chemical properties. However, the usage of abstract methods to study inherently spatial data hinders the exploration and poses a considerable workload. sMolBoxes link quantitative analysis of a user-defined set of properties with interactive 3D visualizations. They enable visual explanations of molecular behaviors, which lead to an efficient discovery of biochemically significant parts of the MD simulation. sMolBoxes follow a node-based model for flexible definition, combination, and immediate evaluation of properties to be investigated. Progressive analytics enable fluid switching between multiple properties, which facilitates hypothesis generation. Each sMolBox provides quick insight to an observed property or function, available in more detail in the bigBox View. The case studies illustrate that even with relatively few sMolBoxes, it is possible to express complex analytical tasks, and their use in exploratory analysis is perceived as more efficient than traditional scripting-based methods.

Dynamic Long-Range Interactions Influence Substrate Binding and Catalysis by Human Histidine Triad Nucleotide-Binding Proteins (HINTs), Key Regulators of Multiple Cellular Processes and Activators of Antiviral ProTides

Human histidine triad nucleotide-binding (hHINT) proteins catalyze nucleotide phosphoramidase and acyl-phosphatase reactions that are essential for the activation of antiviral proTides, such as Sofosbuvir and Remdesivir. hHINT1 and hHINT2 are highly homologous but exhibit disparate roles as regulators of opioid tolerance (hHINT1) and mitochondrial activity (hHINT2). NMR studies of hHINT1 reveal a pair of dynamic surface residues (Q62, E100), which gate a conserved water channel leading to the active site 13 Å away. hHINT2 crystal structures identify analogous residues (R99, D137) and water channel. hHINT1 Q62 variants significantly alter the steady-state k_cat and K_m for turnover of the fluorescent substrate (TpAd), while stopped-flow kinetics indicate that K_D also changes. hHINT2, like hHINT1, exhibits a burst phase of adenylation, monitored by fluorescent tryptamine release, prior to rate-limiting hydrolysis and nucleotide release. hHINT2 exhibits a much smaller burst-phase amplitude than hHINT1, which is further diminished in hHINT2 R99Q. Kinetic simulations suggest that amplitude variations can be accounted for by a variable fluorescent yield of the E·S complex from changes in the environment of bound TpAd. Isothermal titration calorimetry measurements of inhibitor binding show that these hHINT variants also alter the thermodynamic binding profile. We propose that these altered surface residues engender long-range dynamic changes that affect the orientation of bound ligands, altering the thermodynamic and kinetic characteristics of hHINT active site function. Thus, studies of the cellular roles and proTide activation potential by hHINTs should consider the importance of long-range interactions and possible protein binding surfaces far from the active site.

Crystal Structure and Biochemical Analysis of a Cytochrome P450 CYP101D5 from Sphingomonas echinoides

Cytochrome P450 enzymes (CYPs) are heme-containing enzymes that catalyze hydroxylation with a variety of biological molecules. Despite their diverse activity and substrates, the structures of CYPs are limited to a tertiary structure that is similar across all the enzymes. It has been presumed that CYPs overcome substrate selectivity with highly flexible loops and divergent sequences around the substrate entrance region. Here, we report the newly identified CYP101D5 from Sphingomonas echinoides. CYP101D5 catalyzes the hydroxylation of β-ionone and flavonoids, including naringenin and apigenin, and causes the dehydrogenation of α-ionone. A structural investigation and comparison with other CYP101 families indicated that spatial constraints at the substrate-recognition site originate from the B/C loop. Furthermore, charge distribution at the substrate binding site may be important for substrate selectivity and the preference for CYP101D5.

Virtual Quasi-2D Intermediates as Building Blocks for Plausible Structural Models of Amyloid Fibrils from Proteins with Complex Topologies: A Case Study of Insulin

Conformational transitions of globular proteins into amyloid fibrils are complex multistage processes exceedingly challenging to simulate using molecular dynamics (MD). Slow monomer diffusion rates and rugged free energy landscapes disfavor swift self-assembly of orderly amyloid architectures within timescales accessible to all-atom MD. Here, we conduct a multiscale MD study of the amyloidogenic self-assembly of insulin: a small protein with a complex topology defined by two polypeptide chains interlinked by three disulfide bonds. To avoid kinetic traps, unconventional preplanarized insulin conformations are used as amyloid building blocks. These starting conformers generated through uniaxial compression of the native monomer in various spatial directions represent 6 distinct (out of 16 conceivable) two-dimensional (2D) topological classes varying in N-/C-terminal segments of insulin's A- and B-chains being placed inside or outside of the central loop constituted by the middle sections of both chains and Cys7A-Cys7B/Cys19B-Cys20A disulfide bonds. Simulations of the fibrillar self-assembly are initiated through a biased in-register alignment of two, three, or four layers of flat conformers belonging to a single topological class. The various starting topologies are conserved throughout the self-assembly process resulting in polymorphic amyloid fibrils varying in structural features such as helical twist, presence of cavities, and overall stability. Some of the protofilament structures obtained in this work are highly compatible with the earlier biophysical studies on insulin amyloid and high-resolution studies on insulin-derived amyloidogenic peptide models postulating the presence of steric zippers. Our approach provides in silico means to study amyloidogenic tendencies and viable amyloid architectures of larger disulfide-constrained proteins with complex topologies.

A role for the periplasmic adaptor protein AcrA in vetting substrate access to the RND efflux transporter AcrB

Tripartite resistance-nodulation-division (RND) efflux pumps, such as AcrAB-TolC of Salmonella Typhimurium, contribute to antibiotic resistance and comprise an inner membrane RND-transporter, an outer membrane factor, and a periplasmic adaptor protein (PAP). The role of the PAP in the assembly and active transport process remains poorly understood. Here, we identify the functionally critical residues involved in PAP-RND-transporter binding between AcrA and AcrB and show that the corresponding RND-binding residues in the closely related PAP AcrE, are also important for its interaction with AcrB. We also report a residue in the membrane-proximal domain of AcrA, that when mutated, differentially affects the transport of substrates utilising different AcrB efflux channels, namely channels 1 and 2. This supports a potential role for the PAP in sensing the substrate-occupied state of the proximal binding pocket of the transporter and substrate vetting. Understanding the PAP’s role in the assembly and function of tripartite RND pumps can guide novel ways to inhibit their function to combat antibiotic resistance.

The PshX subunit of the photochemical reaction center from Heliobacterium modesticaldum acts as a low-energy antenna

The anoxygenic phototrophic bacterium Heliobacterium modesticaldum contains a photochemical reaction center protein complex (called the HbRC) consisting of a homodimer of the PshA polypeptide and two copies of a newly discovered polypeptide called PshX, which is a single transmembrane helix that binds two bacteriochlorophyll g molecules. To assess the function of PshX, we produced a ∆pshX strain of Hbt. modesticaldum by leveraging the endogenous Hbt. modesticaldum Type I-A CRISPR-Cas system to aid in mutant selection. We optimized this system by separating the homologous recombination and CRISPR-based selection steps into two plasmid transformations, allowing for markerless gene replacement. Fluorescence and low-temperature absorbance of the purified HbRC from the wild-type and ∆pshX strains showed that the bacteriochlorophylls bound by PshX have the lowest site energies in the entire HbRC. This indicates that PshX acts as a low-energy antenna subunit, participating in entropy-assisted uphill energy transfer toward the P₈₀₀ special bacteriochlorophyll g pair. We further discuss the role that PshX may play in stability of the HbRC, its conservation in other heliobacterial species, and the evolutionary pressure to produce and maintain single-TMH subunits in similar locations in other reaction centers.

SARS-CoV-2 spike evolutionary behaviors; simulation of N501Y mutation outcomes in terms of immunogenicity and structural characteristic

Since the emergence of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2), a large number of mutations in its genome have been reported. Some of the mutations occur in noncoding regions without affecting the pathobiology of the virus, while mutations in coding regions are significant. One of the regions where a mutation can occur, affecting the function of the virus is at the receptor‐binding domain (RBD) of the spike protein. RBD interacts with angiotensin‐converting enzyme 2 (ACE2) and facilitates the entry of the virus into the host cells. There is a lot of focus on RBD mutations, especially the displacement of N501Y which is observed in the UK/Kent, South Africa, and Brazilian lineages of SARS‐CoV‐2. Our group utilizes computational biology approaches such as immunoinformatics, protein–protein interaction analysis, molecular dynamics, free energy computation, and tertiary structure analysis to disclose the consequences of N501Y mutation at the molecular level. Surprisingly, we discovered that this mutation reduces the immunogenicity of the spike protein; also, displacement of Asn with Tyr reduces protein compactness and significantly increases the stability of the spike protein and its affinity to ACE2. Moreover, following the N501Y mutation secondary structure and folding of the spike protein changed dramatically.

Structural Analysis of the Ancestral Haloalkane Dehalogenase AncLinB-DmbA

Haloalkane dehalogenases (EC 3.8.1.5) play an important role in hydrolytic degradation of halogenated compounds, resulting in a halide ion, a proton, and an alcohol. They are used in biocatalysis, bioremediation, and biosensing of environmental pollutants and also for molecular tagging in cell biology. The method of ancestral sequence reconstruction leads to prediction of sequences of ancestral enzymes allowing their experimental characterization. Based on the sequences of modern haloalkane dehalogenases from the subfamily II, the most common ancestor of thoroughly characterized enzymes LinB from Sphingobium japonicum UT26 and DmbA from Mycobacterium bovis 5033/66 was in silico predicted, recombinantly produced and structurally characterized. The ancestral enzyme AncLinB-DmbA was crystallized using the sitting-drop vapor-diffusion method, yielding rod-like crystals that diffracted X-rays to 1.5 Å resolution. Structural comparison of AncLinB-DmbA with their closely related descendants LinB and DmbA revealed some differences in overall structure and tunnel architecture. Newly prepared AncLinB-DmbA has the highest active site cavity volume and the biggest entrance radius on the main tunnel in comparison to descendant enzymes. Ancestral sequence reconstruction is a powerful technique to study molecular evolution and design robust proteins for enzyme technologies.

Development of an optogenetic gene sensitive to daylight and its implications in vision restoration

Optogenetic gene-mediated therapy for restoring vision is thought to be a useful treatment for blind patients. However, light sensitivity achieved using this gene therapy is inferior to that of daylight vision. To increase light sensitivity, we designed three mutants using a bioinformatics approach. Nucleotide sequences encoding two sites in the extracellular loops (ex1, ex3) of mVChR1 close to simulated ion-conducting pathways were replaced by homologous amino acid-encoding sequences of ChR1 or ChR2. The light sensitivity of ex3mV1 was higher than that of mVChR1 at 405–617 nm. Visual responses were restored in Royal College of Surgeons rats with genetically degenerating photoreceptor cells transfected with ex3mV1Co, wherein transmembrane of sixth (TM6) in ex3mV1 was additionally replaced with the corresponding domain of CoChR; these rats responded to light in the order of μW/mm². Thus, ex3mV1Co might be useful for the restoration of advanced visual function.

Insight into molecular characteristics of SARS-CoV-2 spike protein following D614G point mutation, a molecular dynamics study

Undoubtedly, the SARS-CoV-2 has become a major concern for all societies due to its catastrophic effects on public health. In addition, mutations and changes in the structure of the virus make it difficult to design effective treatment. Moreover, the amino acid sequence of a protein is a major factor in the formation of the second and tertiary structure in a protein. Amino acid replacement can have noticeable effects on the folding of a protein, especially if an asymmetric change (substitution of polar residue with non-polar, charged with an uncharged, positive charge with a negative charge, or large residue with small residue) occurs. D614G as a spike mutant of SARS-CoV-2 previously identified as an associated risk factor with a high mortality rate of this virus. Using structural bioinformatics, our group determined that D614G mutation could cause extensive changes in SARS-CoV-2 behavior including the secondary structure, receptor binding pattern, 3D conformation, and stability of it.

Communicated by Ramaswamy H. Sarma

Engineering enzyme access tunnels

Enzymes are efficient and specific catalysts for many essential reactions in biotechnological and pharmaceutical industries. Many times, the natural enzymes do not display the catalytic efficiency, stability or specificity required for these industrial processes. The current enzyme engineering methods offer solutions to this problem, but they mainly target the buried active site where the chemical reaction takes place. Despite being many times ignored, the tunnels and channels connecting the environment with the active site are equally important for the catalytic properties of enzymes. Changes in the enzymatic tunnels and channels affect enzyme activity, specificity, promiscuity, enantioselectivity and stability. This review provides an overview of the emerging field of enzyme access tunnel engineering with case studies describing design of all the aforementioned properties. The software tools for the analysis of geometry and function of the enzymatic tunnels and channels and for the rational design of tunnel modifications will also be discussed. The combination of new software tools and enzyme engineering strategies will provide enzymes with access tunnels and channels specifically tailored for individual industrial processes.

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.