Structural determinants of antimicrobial and antiplasmodial activity and selectivity in histidine-rich amphipathic cationic peptides

Designed histidine-rich amphipathic cationic peptides, such as LAH4, have enhanced membrane disruption and antibiotic properties when the peptide adopts an alignment parallel to the membrane surface. Although this was previously achieved by lowering the pH, here we have designed a new generation of histidine-rich peptides that adopt a surface alignment at neutral pH. In vitro, this new generation of peptides are powerful antibiotics in terms of the concentrations required for antibiotic activity; the spectrum of target bacteria, fungi, and parasites; and the speed with which they kill. Further modifications to the peptides, including the addition of more hydrophobic residues at the N terminus, the inclusion of a helix-breaking proline residue or using D-amino acids as building blocks, modulated the biophysical properties of the peptides and led to substantial changes in toxicity to human and parasite cells but had only a minimal effect on the antibacterial and antifungal activity. Using a range of biophysical methods, in particular solid-state NMR, we show that the peptides are highly efficient at disrupting the anionic lipid component of model membranes. However, we also show that effective pore formation in such model membranes may be related to, but is not essential for, high antimicrobial activity by cationic amphipathic helical peptides. The information in this study comprises a new layer of detail in the understanding of the action of cationic helical antimicrobial peptides and shows that rational design is capable of producing potentially therapeutic membrane active peptides with properties tailored to their function.     

Structural characterization of the large soluble oligomers of the GTPase effector domain of dynamin

Dynamin, a protein playing crucial roles in endocytosis, oligomerizes to form spirals around the necks of incipient vesicles and helps their scission from membranes. This oligomerization is known to be mediated by the GTPase effector domain (GED). Here we have characterized the structural features of recombinant GED using a variety of biophysical methods. Gel filtration and dynamic light scattering experiments indicate that in solution, the GED has an intrinsic tendency to oligomerize. It forms large soluble oligomers (molecular mass > 600 kDa). Interestingly, they exist in equilibrium with the monomer, the equilibrium being largely in favour of the oligomers. This equilibrium, observed for the first time for GED, may have regulatory implications for dynamin function. From the circular dichroism measurements the multimers are seen to have a high helical content. From multidimensional NMR analysis we have determined that about 30 residues in the monomeric units constituting the oligomers are flexible, and these include a 17 residue stretch near the N-terminal. This contains two short segments with helical propensities in an otherwise dynamic structure. Negatively charged SDS micelles cause dissociation of the oligomers into monomers, and interestingly, the helical characteristics of the oligomer are completely retained in the individual monomers. The segments along the chain that are likely to form helices have been predicted from five different algorithms, all of which identify two long stretches. Surface electrostatic potential calculation for these helices reveals that there is a distribution of neutral, positive and negative potentials, suggesting that both electrostatic and hydrophobic interactions could be playing important roles in the oligomer core formation. A single point mutation, I697A, in one of the helices inhibited oligomerization quite substantially, indicating firstly, a special role of this residue, and secondly, a decisive, though localized, contribution of hydrophobic interaction in the association process.     

Characterization of Molecular Determinants of the Conformational Stability of Macrophage Migration Inhibitory Factor: Leucine 46 Hydrophobic Pocket

Macrophage Migration Inhibitory Factor (MIF) is a key mediator of inflammatory responses and innate immunity and has been implicated in the pathogenesis of several inflammatory and autoimmune diseases. The oligomerization of MIF, more specifically trimer formation, is essential for its keto-enol tautomerase activity and probably mediates several of its interactions and biological activities, including its binding to its receptor CD74 and activation of certain signaling pathways. Therefore, understanding the molecular factors governing the oligomerization of MIF and the role of quaternary structure in modulating its structural stability and multifunctional properties is crucial for understanding the function of MIF in health and disease. Herein, we describe highly conserved intersubunit interactions involving the hydrophobic packing of the side chain of Leu46 onto the β-strand β3 of one monomer within a hydrophobic pocket from the adjacent monomer constituted by residues Arg11, Val14, Phe18, Leu19, Val39, His40, Val41, Val42, and Pro43. To elucidate the structural significance of these intersubunit interactions and their relative contribution to MIF’s trimerization, structural stability and catalytic activity, we generated three point mutations where Leu46 was replaced by glycine (L46G), alanine (L46A) and phenylalanine (L46F), and their structural properties, stability, oligomerization state, and catalytic activity were characterized using a battery of biophysical methods and X-ray crystallography. Our findings provide new insights into the role of the Leu46 hydrophobic pocket in stabilizing the conformational state of MIF in solution. Disrupting the Leu46 hydrophobic interaction perturbs the secondary and tertiary structure of the protein but has no effect on its oligomerization state.     

Molecular dynamics investigation of the influence of anionic and zwitterionic interfaces on antimicrobial peptides' structure: implications for peptide toxicity and activity

Molecular dynamics simulations of three related helical antimicrobial peptides have been carried out in zwitterionic diphosphocholine (DPC) micelles and anionic sodiumdodecylsulfate (SDS) micelles. These systems can be considered as model mammalian and bacterial membrane interfaces, respectively. The goal of this study is to dissect the differences in peptide composition which make the mutant peptides (novispirin-G10 and novispirin-T7) less toxic than the parent peptide ovispirin (OVIS), although all three peptides have highly antibacterial properties. Compared to G10 and T7, OVIS inserts deepest into the DPC micelle. This correlates well with the lesser toxicity of G10 and T7. There is strong evidence which suggests that synergistic binding of hydrophobic residues drives binding of OVIS to the micelle. The helical content of G10 and T7 is reduced in the presence of DPC, and this leads to less amphipathic peptide structures, which bind weakly to the micelle. Simulations in SDS were carried out to compare the influence of membrane electrostatics on peptide structure. All three peptides bound strongly to SDS, and retained helical form. This corresponds well with their equally potent antibacterial properties. Based on the simulations, we argue that secondary structure stability often leads to toxic properties. We also propose that G10 and T7 operate by the carpet mechanism of cell lysis. Toxicity of peptides operating by the carpet mechanism can be attenuated by reducing the peptide helical content. The simulations successfully capture experimental binding states, and the different depths of binding of the three peptides to the two micelles correlate with their antibacterial and toxic properties.     

Amyloid oligomer formation probed by water proton magnetic resonance spectroscopy

Formation of amyloid oligomers, the most toxic species of amyloids in degenerative diseases, is critically coupled to the interplay with surrounding water. The hydrophobic force driving the oligomerization causes water removal from interfaces, changing the surface-hydration properties. Here, we show that such effects alter the magnetic relaxation response of local water in ways that may enable oligomer detection. By using water proton magnetic resonance spectroscopy, we measured significantly longer transverse magnetic relaxation (T₂) times in mixtures of serum and amyloidogenic Aβ_1-42 peptides versus similar concentration solutions of serum and nonamyloidogenic scrambled Aβ_42-1 peptides. Immunochemistry with oligomer-specific antibodies, electron microscopy and computer simulations demonstrated that the hyperintense magnetic signal correlates with Aβ_1-42 oligomerization. Finding early biophysical markers of the oligomerization process is crucial for guiding the development of new noninvasive imaging techniques, enabling timely diagnosis of amyloid-related diseases and pharmacological intervention.     

Relative spatial positions of tryptophan and cationic residues in helical membrane-active peptides determine their cytotoxicity

The cytotoxic activity of 10 analogs of the idealized amphipathic helical 21-mer peptide (KAAKKAA)3, where three of the Ala residues at different positions have been replaced with Trp residues, has been investigated. The peptide's cytotoxic activity was found to be markedly dependent upon the position of the Trp residues within the hydrophobic sector of an idealized α-helix. The peptides with Trp residues located opposite the cationic sector displayed no antitumor activity, whereas those peptides with two or three Trp residues located adjacent to the cationic sector exhibited high cytotoxic activity when tested against three different cancer cell lines. Dye release experiments revealed that in contrast to the peptides with Trp residues located opposite the cationic sector, the peptides with Trp residues located adjacent to the cationic sector induced a strong permeabilizing activity from liposomes composed of a mixture of zwitterionic phosphatidylcholine and negatively charged phosphatidylserine (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)) (2:1) but not from liposomes composed of zwitterionic phosphatidylcholine, POPC. Fluorescence blue shift and quenching experiments revealed that Trp residues inserted deeper into the hydrophobic environment of POPC/POPS liposomes for peptides with high cytotoxic activity. Through circular dichroism studies, a correlation between the cytotoxic activity and the α-helical propensity was established. Structural studies of one inactive and two active peptides in the presence of micelles using NMR spectroscopy showed that only the active peptides adopted highly coiled to helical structures when bound to a membrane surface.     

Role of recurrent hydrophobic residues in catalysis of helix formation by T cell-presented peptides in the presence of lipid vesicles

We tested the hypothesis that the recurrence of hydrophobic amino acids in a polypeptide at positions falling in an axial, hydrophobic strip if the sequence were coiled as an alpha helix, can lead to helical nucleation on a hydrophobic surface. The hydrophobic surface could anchor such residues, whereas the peptide sequence grows in a helical configuration that is stabilized by hydrogen bonds among carbonyl and amido NH groups along the peptidyl backbone of the helix, and by other intercycle interactions among amino acid side chains. Such bound, helical structures might protect peptides from proteases and/or facilitate transport to a MHC-containing compartment and thus be reflected in the selection of T cell-presented segments. Helical structure in a series of HPLC-purified peptides was estimated from circular dichroism measurements in: 1) 0.01 M phosphate buffer, pH 7.0, 2) that buffer with 45% trifluoroethanol (TFE), and 3) that buffer with di-O-hexadecyl phosphatidylcholine vesicles. By decreasing the dielectric constant of the buffer, TFE enhances intrapeptide interactions generally, whereas the lipid vesicles only provide a surface for hydrophobic interactions. The peptides varied in their strip-of-helix hydrophobicity indices (SOHHI; the mean Kyte-Doolittle hydrophobicities of residues in an axial strip of an alpha helix) and in proline content. Structural order for peptides with helical circular dichroism spectra was estimated as percentage helicity from circular dichroism theta 222 nm values and peptide concentration. A prototypic alpha helical peptide with three cycles plus two amino acids and an axial hydrophobic strip of four leucyl residues (SOHHI = 3.8) was disordered in phosphate buffer, 58% helical in that buffer with 48% TFE, and 36% helical in that buffer with vesicles. Percentage helicity in the presence of vesicles of the subset of peptides without proline followed their SOHHI values. Peptides with multiple prolyl residues had circular dichroism spectra with strong signals, but since they did not have altered spectra in the presence of vesicles relative to phosphate buffer alone, the hydrophobic surface of the vesicle did not appear to stabilize those structures.     

Characterization of the Oligomerization and Aggregation of Human Serum Amyloid A

The fibrillation of Serum Amyloid A (SAA) – a major acute phase protein – is believed to play a role in the disease Amyloid A (AA) Amyloidosis. To better understand the amyloid formation pathway of SAA, we characterized the oligomerization, misfolding, and aggregation of a disease-associated isoform of human SAA – human SAA1.1 (hSAA1.1) – using techniques ranging from circular dichroism spectroscopy to atomic force microscopy, fluorescence spectroscopy, immunoblot studies, solubility measurements, and seeding experiments. We found that hSAA1.1 formed alpha helix-rich, marginally stable oligomers in vitro on refolding and cross-beta-rich aggregates following incubation at 37°C. Strikingly, while hSAA1.1 was not highly amyloidogenic in vitro, the addition of a single N-terminal methionine residue significantly enhanced the fibrillation propensity of hSAA1.1 and modulated its fibrillation pathway. A deeper understanding of the oligomerization and fibrillation pathway of hSAA1.1 may help elucidate its pathological role.     

Removing an interhelical salt bridge abolishes coiled-coil formation in a de novo designed peptide

Alpha-helical coiled coils represent a common protein oligomerization motif that are mainly stabilized by hydrophobic interactions occurring along their coiled-coil interface, the so-called hydrophobic seam. We have recently de novo designed and optimized a series of two-heptad repeat long coiled-coil peptides which are further stabilized by a complex network of inter- and intrahelical salt bridges. Here we have extended the de novo design of such two heptad-repeat long peptides by removing the central and most important g-e' Arg to Glu (g-e'RE) ionic interhelical interaction and replacing these residues by alanine residues. The effect of the missing interhelical ionic interaction on coiled-coil formation and stability has been analyzed by CD spectroscopy, analytical ultracentrifugation, and X-ray crystallography. We show that the peptide, while being highly alpha-helical, is no longer able to form a parallel coiled-coil structure but rather assumes an octameric globular helical assembly devoid of any coiled-coil interactions.     

α-Synuclein Oligomers Induced by Docosahexaenoic Acid Affect Membrane Integrity

A key feature of Parkinson disease is the aggregation of α-synuclein and its intracellular deposition in fibrillar form. Increasing evidence suggests that the pathogenicity of α-synuclein is correlated with the activity of oligomers formed in the early stages of its aggregation process. Oligomers toxicity seems to be associated with both their ability to bind and affect the integrity of lipid membranes. Previously, we demonstrated that α-synuclein forms oligomeric species in the presence of docosahexaenoic acid and that these species are toxic to cells. Here we studied how interaction of these oligomers with membranes results in cell toxicity, using cellular membrane-mimetic and cell model systems. We found that α-synuclein oligomers are able to interact with large and small unilamellar negatively charged vesicles acquiring an increased amount of α-helical structure, which induces small molecules release. We explored the possibility that oligomers effects on membranes could be due to pore formation, to a detergent-like effect or to fibril growth on the membrane. Our biophysical and cellular findings are consistent with a model where α-synuclein oligomers are embedded into the lipid bilayer causing transient alteration of membrane permeability.     

Antimicrobial activity profiles of Amphiphilic Xanthone derivatives are a function of their molecular Oligomerization

Currently, membrane-targeting small antimicrobial peptidomimetics (SAP) are important in antibiotic development because bacteria appear to develop resistance to these surface-active compounds less readily. However, the molecular membrane-targeting action of SAPs has received little attention. In this study, we investigated the effect of oligomerization of amphiphilic xanthone, a model SAP, on its antimicrobial properties against both Gram-positive and Gram-negative bacteria. First, oligomer formation by an amphiphilic xanthone, compound 2 (also coded as AM052), was investigated via solution-state nuclear magnetic resonance (NMR) spectroscopy. Then, the effects of oligomerization on membrane disruption were further studied via biophysical approaches. The results showed that the antimicrobial activities of SAPs develop in several stages: oligomer formation in aqueous solution, initial binding of oligomers to the membrane-water interface followed by insertion into the membrane bilayer, aggregation of antimicrobial oligomers in the membrane, and induced membrane leakage. Ultimately, the presence of the oligomers in the bacterial membrane leads to decreased membrane fluidity and bacterial cell death. Interestingly, the early formation of large oligomers leads to stronger membrane disruption and more rapid bacterial killing. However, reduced antimicrobial activities against Gram-negative bacteria were observed for compounds that formed larger oligomers because the LPS layer acts as a barrier to large complexes. Taken together, our results suggest that oligomerization of SAPs has a strong impact on their antimicrobial properties.     

Self-assembly and structural characterization of Echinococcus granulosus antigen B recombinant subunit oligomers

Echinococcus granulosus antigen B is an oligomeric protein of 120-160 kDa composed by 8-kDa (AgB8) subunits. Here, we demonstrated that the AgB8 recombinant subunits AgB8/1, AgB8/2 and AgB8/3 are able to self-associate into high order homo-oligomers, showing similar properties to that of parasite-produced AgB, making them valuable tools to study AgB structure. Dynamic light scattering, size exclusion chromatography and cross-linking assays revealed approximately 120- to 160-kDa recombinant oligomers, with a tendency to form populations with different aggregation states. Recombinant oligomers showed helical circular dichroism spectra and thermostability similar to those of purified AgB. Cross-linking and limited proteolysis experiments indicated different degrees of stability and compactness between the recombinant oligomers, with the AgB8/3 one showing a more stable and compact structure. We have also built AgB8 subunit structural models in order to predict the surfaces possibly involved in electrostatic and hydrophobic interactions during oligomerization.     

Origin of low mammalian cell toxicity in a class of highly active antimicrobial amphipathic helical peptides

We recently described a novel antimicrobial peptide, RTA3, derived from the commensal organism Streptococcus mitis, with strong anti-Gram-negative activity, low salt sensitivity, and minimal mammalian cell toxicity in vitro and in vivo. This peptide conforms to the positively charged, amphipathic helical peptide motif, but has a positively charged amino acid (Arg-5) on the nonpolar face of the helical structure that is induced upon membrane binding. We surmised that disruption of the hydrophobic face with a positively charged residue plays a role in minimizing eukaryotic cell toxicity, and we tested this using a mutant with an R5L substitution. The greatly enhanced toxicity in the mutant peptide correlated with its ability to bind and adopt helical conformations upon interacting with neutral membranes; the wild type peptide RTA3 did not bind to neutral membranes (binding constant reduced by at least 1000-fold). Spectroscopic analysis indicates that disruption of the hydrophobic face of the parent peptide is accommodated in negatively charged membranes without partial peptide unfolding. These observations apply generally to amphipathic helical peptides of this class as we obtained similar results with a peptide and mutant pair (Chen, Y., Mant, C. T., Farmer, S. W., Hancock, R. E., Vasil, M. L., and Hodges, R. S. (2005) J. Biol. Chem. 280, 12316-12329) having similar structural properties. In contrast to previous interpretations, we demonstrate that these peptides simply do not bind well to membranes (like those of eukaryotes) with exclusively neutral lipids in their external bilayer leaflet. We highlight a significant role for tryptophan in promoting binding of amphipathic helical peptides to neutral bilayers, augmenting the arsenal of strategies to reduce mammalian toxicity in antimicrobial peptides.     

Zinc Induces Temperature-Dependent Reversible Self-Assembly of Tau

Tau is an intrinsically disordered microtubule-associated protein that is implicated in several neurodegenerative disorders called tauopathies. In these diseases, Tau is found in the form of intracellular inclusions that consist of aggregated paired helical filaments (PHFs) in neurons. Given the importance of this irreversible PHF formation in neurodegenerative disease, Tau aggregation has been extensively studied. Several different factors, such as mutations or post translational modifications, have been shown to influence the formation of late-stage non-reversible Tau aggregates. It was recently shown that zinc ions accelerated heparin-induced oligomerization of Tau constructs. Indeed, in vitro studies of PHFs have usually been performed in the presence of additional co-factors, such as heparin, in order to accelerate their formation. Using turbidimetry, we investigated the impact of zinc ions on Tau in the absence of heparin and found that zinc is able to induce a temperature-dependent reversible oligomerization of Tau. The obtained oligomers were not amyloid-like and dissociated instantly following zinc chelation or a temperature decrease. Finally, a combination of isothermal titration calorimetry and dynamic light scattering experiments showed zinc binding to a high-affinity binding site and three low-affinity sites on Tau, accompanied by a change in Tau folding. Altogether, our findings stress the importance of zinc in Tau oligomerization. This newly identified Zn-induced oligomerization mechanism may be a part of a pathway different of and concurrent to Tau aggregation cascade leading to PHF formation.     

Differences in the hydrophobic properties of discrete alpha-helical domains of rat and human apolipoprotein A-IV

The primary structure of human apolipoprotein A-IV has been analyzed to identify amphipathic alpha-helical regions and compare their hydrophobic properties with analogous segments of rat apolipoprotein A-IV. Significant differences were found in the mean hydrophobicity and/or mean helical hydrophobic moment of several discrete alpha-helical domains of rat and human apolipoprotein A-IV, which may determine the differences in their biophysical and biological behavior.     

Dissecting the assembly of Abeta16-22 amyloid peptides into antiparallel beta sheets

Multiple long molecular dynamics simulations are used to probe the oligomerization mechanism of Abeta(16-22) (KLVFFAE) peptides. The peptides, in the monomeric form, adopt either compact random-coil or extended beta strand-like structures. The assembly of the low-energy oligomers, in which the peptides form antiparallel beta sheets, occurs by multiple pathways with the formation of an obligatory alpha-helical intermediate. This observation and the experimental results on fibrillogenesis of Abeta(1-40) and Abeta(1-42) peptides suggest that the assembly mechanism (random coil --> alpha helix --> beta strand) is universal for this class of peptides. In Abeta(16-22) oligomers both interpeptide hydrophobic and electrostatic interactions are critical in the formation of the antiparallel beta sheet structure. Mutations of either hydrophobic or charged residues destabilize the oligomer, which implies that the 16-22 fragments of Arctic (E22G), Dutch (E22Q), and Italian (E22K) mutants are unlikely to form ordered fibrils.     

The crystal structure of the calcium-bound con-G[Q6A] peptide reveals a novel metal-dependent helical trimer

The ability to form and control both secondary structure and oligomerization in short peptides has proven to be challenging owing to the structural instability of such peptides. The conantokin peptides are a family of γ-carboxyglutamic acid containing peptides produced in the venoms of predatory sea snails of the Conus family. They are examples of short peptides that form stable helical structures, especially in the presence of divalent cations. Both monomeric and dimeric conantokin peptides have been identified and represent a new mechanism of helix association, “the metallozipper motif” that is devoid of a hydrophobic interface between monomers. In the present study, a parallel/antiparallel three-helix bundle was identified and its crystal structure determined at high resolution. The three helices are almost perfectly parallel and represent a novel helix–helix association. The trimer interface is dominated by metal chelation between the three helices, and contains no interfacial hydrophobic interactions. It is now possible to produce stable monomeric, dimeric, or trimeric metallozippers depending on the peptide sequence and metal ion. Such structures have important applications in protein design.     

Binding of amphipathic α-helical antimicrobial peptides to lipid membranes: Lessons from temporins B and L

Temporins constitute a family of amphipathic α-helical antimicrobial peptides (AMP) and contain some of the shortest cytotoxic peptides, comprised of only 10-14 residues. General characteristics of temporins parallel those of other AMP, both in terms of structural features and biophysical properties relating to their interactions with membrane lipids, with selective lipid-binding properties believed to underlie the discrimination between target vs host cells. Lipid-binding properties also contribute to the cytotoxicity AMP, causing permeabilization of their target cell membranes. The latter functional property of AMP involves highly interdependent acidic phospholipid-induced conformational changes, aggregation, and formation of toxic oligomers in the membrane. These oligomers are subsequently converted to amyloid-type fibers, as demonstrated for e.g. temporins B and L in our laboratory, and more recently for dermaseptins by Auvynet et al. Amyloid state represents the generic minimum in the folding/aggregation free energy landscape, and for AMP its formation most likely serves to detoxify the peptides, in keeping with the current consensus on mature amyloid being inert and non-toxic. The above scenario is supported by sequence analyses of temporins as well as other amphipathic α-helical AMP belonging to diverse families. Accordingly, sequence comparison identifies ‘conformational switches’, domains with equal probabilities for adopting random coil, α-helical and β-sheet structures. These regions were further predicted also to aggregate and assemble into amyloid β-sheets. Taken together, the lipid-binding properties and structural characterization lend support to the notion that the mechanism of membrane permeabilization by temporins B and L and perhaps of most AMP could be very similar, if not identical, to that of the paradigm amyloid forming cytotoxic peptides, responsible for degenerative cell loss in e.g. prion, Alzheimer's and Parkinson's disease, and type 2 diabetes.     

Secondary structure of cell‐penetrating peptides during interaction with fungal cells

Cell‐penetrating peptides (CPPs) are peptides that cross cell membranes, either alone or while carrying molecular cargo. Although their interactions with mammalian cells have been widely studied, much less is known about their interactions with fungal cells, particularly at the biophysical level. We analyzed the interactions of seven CPPs (penetratin, Pep‐1, MPG, pVEC, TP‐10, MAP, and cecropin B) with the fungal pathogen Candida albicans using experiments and molecular simulations. Circular dichroism (CD) of the peptides revealed a structural transition from a random coil or weak helix to an α‐helix occurs for all peptides when the solvent is changed from aqueous to hydrophobic. However, CD performed in the presence of C. albicans cells showed that proximity to the cell membrane is not necessarily sufficient to induce this structural transition, as penetratin, Pep‐1, and MPG did not display a structural shift in the presence of cells. Monte Carlo simulations were performed to further probe the molecular‐level interaction with the cell membrane, and these simulations suggested that pVEC, TP‐10, MAP, and cecropin B strongly penetrate into the hydrophobic domain of the membrane lipid bilayer, inducing a transition to an α‐helical conformation. In contrast, penetratin, Pep‐1 and MPG remained in the hydrophilic region without a shift in conformation. The experimental data and MC simulations combine to explain how peptide structure affects their interaction with cells and their mechanism of translocation into cells (direct translocation vs. endocytosis). Our work also highlights the utility of combining biophysical experiments, biological experiments, and molecular modeling to understand biological phenomena.