Excessive production of transforming growth factor β1 (TGF-β1) in activated hepatic stellate cells (aHSCs) promotes liver fibrosis by activating the TGF-β1/Smad signaling pathway. Thus, specifically inhibiting the pro-fibrotic activity of TGF-β1 in aHSCs is an ideal strategy for treating liver fibrosis. Overexpression of platelet-derived growth factor β receptor (PDGFβR) has been demonstrated on the surface of aHSCs relative to normal cells in liver fibrosis. Interferon-gamma peptidomimetic (mIFNγ) and truncated TGF-β receptor type II (tTβRII) inhibit the TGF-β1/Smad signaling pathway by different mechanisms. In this study, we designed a chimeric protein by the conjugation of (1) mIFNγ and tTβRII coupled via plasma protease-cleavable linker sequences (FNPKTP) to (2) PDGFβR-recognizing peptide (BiPPB), namely BiPPB-mIFNγ-tTβRII. This novel protein BiPPB-mIFNγ-tTβRII was effectively prepared using Escherichia coli expression system. The active components BiPPB-mIFNγ and tTβRII were slowly released from BiPPB-mIFNγ-tTβRII by hydrolysis using the plasma protease thrombin in vitro. Moreover, BiPPB-mIFNγ-tTβRII highly targeted to fibrotic liver tissues, markedly ameliorated liver morphology and fibrotic responses in chronic liver fibrosis mice by both inhibiting the phosphorylation of Smad2/3 and inducing the expression of Smad7. Meanwhile, BiPPB-mIFNγ-tTβRII markedly reduced the deposition of collagen fibrils and expression of fibrosis-related proteins in acute liver fibrosis mice. Furthermore, BiPPB-mIFNγ-tTβRII showed a good safety performance in both liver fibrosis mice. Taken together, BiPPB-mIFNγ-tTβRII improved the in vivo anti-liver fibrotic activity due to its high fibrotic liver-targeting potential and the dual inhibition of the TGF-β1/Smad signaling pathway, which may be a potential candidate for targeting therapy on liver fibrosis.

Recombinant fluorescent fusion proteins are fundamental to advancing many aspects of protein science. Such proteins are typically used to enable the visualization of functional proteins in experimental systems, particularly cell biology. An important problem in biotechnology is the production of functional, soluble proteins. Here we report the use of mCherry-fusions of soluble, cysteine-rich, Leptospira-secreted exotoxins in the PF07598 gene family, the so-called virulence modifying (VM) proteins. The mCherry fusion proteins facilitated the visual detection of pink colonies of the VM proteins (LA3490 and LA1402) and following them through lysis and sequential chromatography steps. CD-spectroscopy analysis confirmed the stability and robustness of the mCherry-fusion protein, with a structure comparable to AlphaFold structural predictions. LA0591, a unique member of the PF07598 gene family that lacks N-terminal ricin B-like domains, was produced without mCherry tag that strengthens the recombinant protein production protocol without fusion protein as well. The current study provides the approaches for the synthesis of 50-125 kDa soluble, cysteine-rich, high-quality fast protein liquid chromatography (FPLC)-purified protein, with and without a mCherry tag. The use of mCherry-fusion proteins enables a streamlined, efficient process of protein production and qualitative and quantitative downstream analytical and functional studies. Approaches for troubleshooting and optimization were evaluated to overcome difficulties in recombinant protein expression and purification, demonstrating biotechnology utility in accelerating recombinant protein production.

The lack of specific antiviral therapy and complications associated with the existing peste des petits ruminants (PPR) vaccines accentuates the search of novel antiviral blocking agents in order to curtail the PPR infection at initial level. The synthetic hemagglutinin-neuraminidase (HN) homologous peptides may compete with the natural HN protein of PPR virus for binding to signaling lymphocytic activation molecule (SLAM) receptor, consequently, may disrupt peste des petits ruminants virus (PPRV) at entry level. Therefore, insilico analysis, synthesis, purification and subsequent characterization of HN homologous peptides were conducted in this study. The HN homologous peptides were synthesized by means of solid phase chemistry and were purified by reversed-phase-high performance liquid chromatography. The mass as well as sequence of HN homologous peptides were assessed by mass spectroscopy while its secondary structure was elucidated by circular dichroism spectroscopy. The binding (interaction) efficacy of HN homologous peptides with PPRV antibodies was assessed via indirect enzyme linked immunosorbent assay, visual detection test (red wine to purple), bathochromic shift under UV-Vis spectrophotometry and lateral flow immunochromatographic strip test. The antiviral properties and cytotoxicity of these peptides were also assessed in B95a cell line with changes in cytopathic effect and titer of PPRV (Sungri/96). The presence of green fluorescein isothiocyanate over the B95a cell surface pointed towards the binding of HN homologous peptides with surface SLAM receptor. Moreover, the intact beta sheet configuration in water and lower cytotoxicity [cytotoxic concentration 50 (CC50) > 1000 µg/ml] of these peptides signifies its in vivo use. Among HN homologous peptides, the binding efficacy and antiviral properties of pep A was relatively high in comparison to pep B and Pep ppr peptides. The prerequisite concentration of HN homologous peptides (pep A = 12.5 µg/ml; pep B = 25 µg/ml; pep ppr = 25 µg/ml) to exemplify its antiviral effect was much lower than its CC50 level. Hence, this study signifies the therapeutic potential of synthetic HN homologous peptides.

Chondroitinase ABC I (cABC I) from Proteus vulgaris is an important enzyme in medicinal biotechnology due to its ability to help axon regeneration after spinal cord injury. Its practical application involves solving several problems at the molecular and cellular levels. Structurally, most residues at the C-terminal domain of cABC I are arranged as organized strands, and only a small fraction of residues have helical conformation. The structural and functional features of modified residues on two specific helix fragments have previously been reported. The single mutant M889K has been combined with L679S and L679D mutants to make enzyme variants containing simultaneously modified helix. Here, the pH stability and temperature-based analysis of the transition state structure for the catalysis reaction were investigated. We found that double mutant L679D/M889K is the better choice to use in physiological conditions due to its higher pH stability at physiological pH as well as its different optimum temperature as compared with the (wild-type) WT protein. According to Arrhenius's analysis, the values of the Gibbs free energy of the transition state (∆G#) are not changed upon mutation. However, the relative contribution and absolute values of the enthalpy and entropy change to the total value of ∆G#, varied between the WT and mutants.

Recombinant human keratinocyte growth factor (rhKGF) is a highly aggregation-prone therapeutic protein. The present study aimed to reduce aggregation propensity of rhKGF by engineering the aggregation hotspots. Initially, 21 mutants were designed based on the previously-identified aggregation-prone regions (APRs) and then four of them including mutants No. 4 (L91K, I119K), 7 (V13S, L91K), 14 (L91D, I119D), and 21 (A51E) were selected based on molecular dynamics (MD) simulations for further experimental studies. The recombinantly produced rhKGF and mutants were analyzed regarding secondary structure, thermal stability, aggregation propensity, and biological activity. Far-UV CD spectroscopy showed that the mutants have similar secondary structure with rhKGF. A51E mutant showed enhanced stability and decreased monomer loss under heat stress suggesting its reduced aggregation propensity compared to rhKGF. Mutant No. 14 showed higher stability and less aggregation tendency than mutant No. 4 indicating that only mutations decreasing pI of rhKGF are effective in reducing its aggregation tendency. All of the mutants were at least as potent as rhKGF in stimulating proliferation of MCF-7 epithelial cells. Our results identified A51E as an equally potent, more stable, and less aggregation-prone analog of rhKGF which could be a promising alternative drug candidate for the commercially available rhKGF (Palifermin).

Equine chorionic gonadotropin (eCG) is a glycoprotein hormone widely used in timed artificial ovulation (TAI) and superovulation protocols to improve the reproductive performance in livestock. Until recently, the only eCG products available in the market for veterinary use consisted in partially purified preparations of pregnant mare serum gonadotropin (PMSG). Here, a bioactive recombinant eCG (reCG) produced in suspension CHO-K1 cells was purified employing different chromatographic methods (hydrophobic interaction chromatography and reverse-phase (RP)-HPLC) and compared with a RP-HPLC-purified PMSG. To gain insight into the structural and functional characteristics of reCG, a bioinformatics analysis was performed. An exhaustive characterization comprising the determination of the purity degree, aggregates and nicked forms through SDS-PAGE, RP-HPLC and SEC-HPLC was performed. Higher order structures were studied by fluorescence spectroscopy and SEC-HPLC. Isoforms profile were analyzed by isoelectric focusing. Glycosylation analysis was performed through pulsed amperometric detection and PNGase F treatment following SDS-PAGE and weak anion exchange-HPLC. Slight differences between the purified recombinant hormones were found. However, recombinant molecules and PMSG exhibited variations in the glycosylation pattern. In fact, differences in sialic acid content between two commercial preparations of PMSG were also obtained, which could lead to differences in their biological potency. These results show the importance of having a standardized production process, as occurs in a recombinant protein bioprocess. Besides, our results reflect the importance of the glycan moieties on eCG conformation and hence in its biological activity, preventing denaturing processes such as aggregation.

Cholinesterase enzyme family consists of acetylcholinesterase (AChE, 3.1.1.7), the major enzyme responsible for hydrolysis of acetylcholine at cholinergic synapses, and butyrylcholinesterase (BChE, 3.1.1.8) a detoxification enzyme of plasma. Statins are cholesterol-lowering medications utilized as protective medicaments in stroke and Alzheimer's disease, which cholinesterases are associated with. Thus, in this study, we characterized the inhibitory effects and mechanisms of common statins, rosuvastatin, atorvastatin, simvastatin and lovastatin, on human erythrocyte AChE and purified serum BChE using in vitro and in silico methods. Kinetic assays identified statins as selective non-competitive inhibitors of human serum BChE. The IC50 and Km values were found as 194.7 ± 55.2 µM and 1.03 ± 0.2 µM for rosuvastatin, 492.5 ± 55.1 µM and 7.2 ± 0.3 µM for atorvastatin, 14.2 ± 0.3 µM and 202.7 ± 23.2 µM for lovastatin, and 17.6 ± 0.1 µM and 207.2 ± 13.2 µM for simvastatin, respectively. The compounds did not display considerable inhibition against AChE. Molecular docking predicted good affinity and strong interactions with the BChE active site for atorvastatin and rosuvastatin. Current study identifies rosuvastatin as the most specific and selective inhibitor of human BChE among the tested statins. As selective inhibitors of BChE statins have the potential to be re-evaluated as medicaments due to their pleiotropic effects.

Glycosylation is an important post translational modification in plants. First analysis of N-linked glycosylated proteins of Dioscorea alata using Concanavalin A lectin affinity chromatography enrichment coupled with label free quantification is presented. In total, 114 enriched glycoproteins were detected. Signal P and sub-cellular localization showed 42.2% of proteins to be secretory. These included peroxidases, endochitinases, calreticulin, calnexin, thaumatins and lipid transfer proteins. Gene Ontology and MapMan analysis predicted the enriched glycoproteins to be involved in processes essential for tuber maturation namely: signal transduction, lignification, protein trafficking, endoplasmic reticulum quality control and cell wall remodeling. This was supported by biochemical validation of the essential glycoproteins. Interestingly, out of the two dioscorin isoforms, Dio B was the only N-glycosylated form. In silico analysis showed O-glycosylation sites in the other form, Dio A suggesting its similarity with sporamin, the storage protein of sweet potato. Absence of signal peptide in Dio B and the presence of non-canonical motif hints towards its atypical glycosylation. The analysis revealed that N-glycosylation of Dio B isoform maintains the activities associated with Dioscorin at maturity and provides an overview of protein N-glycosylation, enriching the glycoproteome database of plants especially tubers.

The study of interactions between proteins and surfactants is of relevance in a diverse range of applications including food, enzymatic detergent formulation, and drug delivery. In spite of sodium dodecyl sulfate (SDS)-induced unfolding has been studied in detail at the protein level, deciphering the conformation-activity relationship of a recombinant γ-glutamyltranspeptidase (BlrGGT) from Bacillus licheniformis remains important to understand how the transpeptidase activity is related to its conformation. In this study, we examined the enzyme catalysis and conformational transition of BlrGGT in the presence of SDS. Enzymatic assays showed that the transpeptidase activity of BlrGGT was greatly affected by SDS in a concentration-dependent manner with approximately 90% inactivation at 6 mM. Native polyacrylamide gel electrophoresis of SDS-treated samples clearly revealed that the heterodimeric enzyme was apparently dissociated into two different subunits at concentrations above 2 mM. The study of enzyme kinetics showed that SDS can act as a mixed-type inhibitor to reduce the catalytic efficiency of BlrGGT. Moreover, the t1/2 value of the enzyme at 55 °C was greatly reduced from 495.1 min to 7.4 min in the presence of 1 mM SDS. The I3/I1 ratio of pyrene excimer fluorescence emission changed around 3.7 mM SDS in the absence of BlrGGT and the inflection point of enzyme samples was reduced to less than 2.7 mM. The Far-UV CD spectrum of the native enzyme had two negative peaks at 208 and 222 nm, respectively; however, both negative peaks increased in magnitude with increasing SDS concentration and reached maximal values at above 4.0 mM. The intrinsic fluorescence spectra of tryptophan further demonstrated that the SDS-induced enzyme conformational transition occurred at approximately 5.1 mM. Tween 20 significantly suppressed the interaction of BlrGGT with SDS by forming mixed micelles at a molar ratio of 1.0. Taken together, this study definitely promotes our better understanding of the relationship between the conformation and catalysis of BlrGGT.

Paired box 4 (PAX4) is a pivotal transcription factor involved in pancreatogenesis during embryogenesis, and in adults, it is key for β-cell proliferation and survival. Additionally, PAX4 also functions as a tumor suppressor protein in human melanomas. The present study demonstrates the production of bioactive recombinant human PAX4 transcription factor. At first, the inserts (PAX4 protein-coding sequence having tags at either ends) were cloned in an expression vector to give rise to pET28a(+)-HTN-PAX4 and pET28a(+)-PAX4-NTH genetic constructs, and these were then transformed into Escherichia coli (E. coli) for their expression. The HTN-PAX4 and PAX4-NTH fusion proteins produced were purified with a yield of ~ 3.15 mg and ~ 0.83 mg, respectively, from 1.2 L E. coli culture. Further, the secondary structure retention of the PAX4 fusion proteins and their potential to internalize the mammalian cell and its nucleus was demonstrated. The bioactivity of these fusion proteins was investigated using various assays (cell migration, cell proliferation and cell cycle assays), demonstrating it to function as a tumor suppressor protein. Thus, this macromolecule can prospectively help understand the function of human PAX4 in cellular processes, disease-specific investigations and direct cellular reprogramming.

INTRODUCTION: L-asparaginase (also known as L-ASNase) is a crucial therapeutic enzyme that is widely used in treatment of ALL (acute lymphoblastic leukemia) as a chemotherapeutic drug. Besides, this enzyme is used in the food industry as a food processing reagent to reduce the content of acrylamide in addition to the clinical industry. The improvement of activity and kinetic parameters of the L-ASNase enzyme may lead to higher efficiency resulting in practical achievement. In order to achieve this goal, we chosen glycine residue in position 88 as a potential mutation with advantageous outcomes. METHOD: In this study, firstly to find the appropriate mutation on glycine 88, various in silico analyses, such as MD simulation and molecular docking, were carried out. Then, the rational design was adopted as the best strategy for molecular modifications of the enzyme to improve its enzymatic properties. RESULT: Our in silico findings show that the four mutations G88Q, G88L, G88K, and G88A may be able to increase L-ASNase's asparaginase activity. The catalytic efficiency of each enzyme (kcat/Km) is the most important feature for comparing the catalytic activity of the mutants with the wild type form. The laboratory experiments showed that the kcat/Km for the G88Q mutant is 36.32% higher than the Escherichia coli K12 ASNase II (wild type), which suggests that L-ASNase activity is improved at lower concentration of L-ASN. Kinetic characterization of the mutants L-ASNase activity confirmed the high turnover rate (kcat) with ASN as substrate relative to the wild type enzyme. CONCLUSION: In silico analyses and laboratory experiments demonstrated that the G88Q mutation rather than other mutation (G88L, G88K, and G88A) could improve the kinetics of L-ASNase.

Molecular dynamics (MD) simulations are routinely performed of biomolecules in solution, because this is their native environment. However, the structures used in such simulations are often obtained with X-ray crystallography, which provides the atomic coordinates of the biomolecule in a crystal environment. With the advent of free electron lasers and time-resolved techniques, X-ray crystallography can now also access metastable states that are intermediates in a biochemical process. Such experiments provide additional data, which can be used, for example, to optimize MD force fields. Doing so requires that the simulation of the biomolecule is also performed in the crystal environment. However, in contrast to simulations of biomolecules in solution, setting up a crystal is challenging. In particular, because not all solvent molecules are resolved in X-ray crystallography, adding a suitable number of solvent molecules, such that the properties of the crystallographic unit cell are preserved in the simulation, can be difficult and typically is a trial-and-error based procedure requiring manual interventions. Such interventions preclude high throughput applications. To overcome this bottleneck, we introduce gmXtal, a tool for setting up crystal simulations for MD simulations with GROMACS. With the information from the protein data bank (rcsb.org) gmXtal automatically (i) builds the crystallographic unit cell; (ii) sets the protonation of titratable residues; (iii) builds missing residues that were not resolved experimentally; and (iv) adds an appropriate number of solvent molecules to the system. gmXtal is available as a standalone tool https://gitlab.com/pbuslaev/gmxtal .

γδ T cells, especially Vγ9Vδ2 T cells, play an important role in mycobacterial infection. We have identified some Vγ9Vδ2 T cells that recognize protein/peptide antigens derived from mycobacteria, which may induce protective immune responses to mycobacterial infection. To clarify the structural basis of the molecular recognition mechanism, we tried many methods to express the Vγ9Vδ2 T-cell receptor (TCR). The Vγ9Vδ2 TCR was not expressed well in a prokaryotic expression system or a baculovirus expression system, even after extensive optimization. In a mammalian cell expression system, the Vγ9Vδ2 TCR was expressed in the form of a soluble heterodimer, which was suitable for crystal screening. Reduced-temperature cultivation (cold shock) increased the yield of the recombinant TCR. The recombinant purified TCR was used for crystal trials, and crystals that could be used for X-ray diffraction were obtained. Although we have not yet determined the crystal structure of the Vγ9Vδ2 TCR, we have established a procedure for Vγ9Vδ2 TCR expression and purification, which is useful for basic research and potentially for clinical application.

Monkeypox, a viral zoonotic disease resembling smallpox, has emerged as a significant national epidemic primarily in Africa. Nevertheless, the recent global dissemination of this pathogen has engendered apprehension regarding its capacity to metamorphose into a sweeping pandemic. To effectively combat this menace, a multi-epitope vaccine has been meticulously engineered with the specific aim of targeting the cell envelope protein of Monkeypox virus (MPXV), thereby stimulating a potent immunological response while mitigating untoward effects. This new vaccine uses T-cell and B-cell epitopes from a highly antigenic, non-allergenic, non-toxic, conserved, and non-homologous A30L protein to provide protection against the virus. In order to ascertain the vaccine design with the utmost efficacy, protein-protein docking methodologies were employed to anticipate the intricate interactions with Toll-like receptors (TLR) 2, 3, 4, 6, and 8. This meticulous approach led the researchers to discern an optimal vaccine architecture, bolstered by affirmative prognostications derived from both molecular dynamics (MD) simulations and immune simulations. The current research findings indicate that the peptides ATHAAFEYSK, FFIVVATAAV, and MNSLSIFFV exhibited antigenic properties and were determined to be non-allergenic and non-toxic. Through the utilization of codon optimization and in-silico cloning techniques, our investigation revealed that the prospective vaccine exhibited a remarkable expression level within Escherichia coli. Moreover, upon conducting immune simulations, we observed the induction of a robust immune response characterized by elevated levels of both B-cell and T-cell mediated immunity. Moreover, as the initial prediction with in-silico techniques has yielded promising results these epitope-based vaccines can be recommended to in vitro and in silico studies to validate their immunogenic properties.

Bovine lactoferrin peptide (LFcinB), as an antimicrobial peptide, is expected to be an alternative of antibiotics owing to its broad-spectrum antimicrobial activity and specific mechanism. However, the weak antimicrobial activity, high hemolysis, and poor stability of LFcinB limited its applications in the field of biomedicine, food and agriculture. In order to improve the antimicrobial activity of LFcinB, five mutants were designed rationally, of which mutant LF4 (M10W/P16R/A24L) showed highest antimicrobial activity. The bioinformatics analysis indicated that the improved antimicrobial activity of LF4 was related to its increased cations, higher amphiphilicity and the extension of the β-sheet in the structure. These studies will highlight the important role of bioinformatic tools in designing ideal biopeptides and lay a foundation for further development of antimicrobial peptides.

Acetylcholinesterase (AChE, E.C. 3.1.1.7) termed as the true cholinesterase functions to end cholinergic transmission at synapses. Due to its diverse expression in non-neural tissues such as erythrocytes and bones along with its various molecular forms, researchers seek a non-classical role for this protein. Here, the inhibitory action of unsaturated 18 carbon fatty acids linoleic acid and alpha-linolenic acid and 20 carbon fatty acid arachidonic acid on AChE were investigated. Enzyme activity was measured in kinetic assay method according to Ellman assay utilizing acetylthiocholine. Analysis of the activity data revealed that among the fatty acids examined the IC50 values differed according to the length of the fatty acid and the number of the double bonds. Arachidonic acid, a 20-carbon fatty acid with 4 unsaturated bonds (20:4 n-6, cis 5,8,11,14) displayed an IC50 value of 2.78 µM and Ki value of 396.35 µM. Linoleic acid, an essential 18-carbon fatty acid (18:2 n-6, cis 9,12) had an IC50 value of 7.95 µM and Ki value of 8027.55 µM. The IC50 value of alpha-linolenic acid, 18-carbon fatty acid (18:3 n-3, cis-9,12,15) was found as 179.11 µM. Analysis of the data fit the inhibition mechanism for linoleic, alpha-linolenic and arachidonic acid as mixed-type; non-competitive. Molecular docking complied with these results yielding the best score for arachidonic acid. The alkenyl chain of the fatty acids predictably reached to the catalytic site while the carboxylate strongly interacted with the peripheric anionic site.

Malate is an important material to various industrials and clinical applications. Bacillus subtilis is a widely used biocatalyst tool for chemical production. However, the specific enzymatic properties of malate dehydrogenase from Bacillus subtilis (BsMDH) remain largely unknown. In the present study, BsMDH was cloned, recombinantly expressed and purified to test its enzymatic properties. The molecular weight of single unit of BsMDH was 34,869.7 Da. Matrix-Assisted Laser-Desorption Ionization-Time-of-Flight Mass Spectrometry and gel filtration analysis indicated that the recombinant BsMDH could form dimers. The kcat/Km values of oxaloacetate and NADH were higher than those of malate and NAD+, respectively, indicating a better catalysis in the direction of malate synthesis than the reverse. Furthermore, six BsMDH mutants were constructed with the substitution of amino acids at the coenzyme binding site. Among them, BsMDH-T7 showed a greatly higher affinity and catalysis efficiency to NADPH than NADH with the degree of alteration of 2039, suggesting the shift of the coenzyme dependence from NADH to NADPH. In addition, BsMDH-T7 showed a relatively lower Km value, but a higher kcat and kcat/Km than NADPH-dependent MDHs from Thermus flavus and Corynebacterium glutamicum. Overall, these results indicated that BsMDH and BsMDH-T7 mutant might be promising enzymes for malate production.

Ginkgo seed is an important Chinese medicine and food resource in China, but the toxicity of ginkgo acid in it limits its application. Previous studies have found that salicylic acid decarboxylase (Sdc) has a decarboxylation degradation effect on ginkgo acid. In order to improve the decarboxylation ability of Sdc to Ginkgo acid, 11 residues of the Sdc around the substrate (salicylic acid) were determined as mutation targets according to the analysis of crystal structure of Sdc (PDB ID:6JQX), from Trichosporon moniliiforme WU-0401, and a total of 30 single point mutant enzymes and one compound mutant enzyme were obtained. With Ginkgo acid C15:1 as the substrate, it was found from activity assay that Sdc-Y64T and Sdc-P191A had higher decarboxylation activity, which increased by 105.18% and 116.74% compared with that of wild type Sdc, respectively. The optimal pH for Sdc Y64T and Sdc-P191A to decarboxylate Ginkgo acid C15:1 was 5.5, which is the same as the wild type Sdc. The optimal temperature of Sdc-P191A was 50 °C, which was consistent with that of the wild type Sdc, but the optimal temperature of the mutant Sdc-Y64T was 40 °C, which was 10 °C lower than that of wild type Sdc.

In this contribution I describe my personal memories of my work with Herman Berendsen and the Biophysical chemistry group in Groningen in the period of years from 1967 to 1975.

My memories of professor Herman Berendsen cover roughly two periods during which I had many contacts with him. Between 1966 and 1973 I was his MSc student and later his PhD student in the Department of Biophysical Chemistry at the University of Groningen. The second period started in 1991 when I returned to the University of Groningen as a professor of environmental sciences.