A molecular dynamics study of loop fluctuation in human papillomavirus type 16 virus-like particles: a possible indicator of immunogenicity

Experimental and molecular dynamics simulation studies on the physical properties of three HBc-VLP derivatives as nanoparticle protein vaccine candidates

Self-assembling virus-like particles (VLPs) are promising platforms for vaccine development. However, the unpredictability of the physical properties, such as self-assembly capability, hydrophobicity, and overall stability in engineered protein particles fused with antigens, presents substantial challenges in their downstream processing. We envision that these challenges can be addressed by combining more precise computer-aided molecular dynamics (MD) simulations with experimental studies on the modified products, with more to-date forcefield descriptions and larger models closely resembling real assemblies, realized by rapid advancement in computing technology. In this study, three chimeric designs based on the hepatitis B core (HBc) protein as model vaccine candidates were constructed to study and compare the influence of inserted epitopes as well as insertion strategy on HBc modifications. Large partial VLP models containing 17 chains for the HBc chimeric model vaccines were constructed based on the wild-type (wt) HBc assembly template. The findings from our simulation analysis have demonstrated good consistency with experimental results, pertaining to the surface hydrophobicity and overall stability of the chimeric vaccine candidates. Furthermore, the different impact of foreign antigen insertions on the HBc scaffold was investigated through simulations. It was found that separately inserting two epitopes into the HBc platform at the N-terminal and the major immunogenic regions (MIR) yields better results compared to a serial insertion at MIR in terms of protein structural stability. This study substantiates that an MD-guided design approach can facilitate vaccine development and improve its manufacturing efficiency by predicting products with extreme surface hydrophobicity or structural instability.

Integrating molecular dynamics simulation with small- and wide-angle X-ray scattering to unravel the flexibility, antigen-blocking, and protease-restoring functions in a hindrance-based pro-antibody

Spatial hindrance-based pro-antibodies (pro-Abs) are engineered antibodies to reduce monoclonal antibodies' (mAbs) on-target toxicity using universal designed blocking segments that mask mAb antigen-binding sites through spatial hindrance. By linking through protease substrates and linkers, these blocking segments can be removed site-specifically. Although many types of blocking segments have been developed, such as coiled-coil and hinge-based Ab locks, the molecular structure of the pro-Ab, particularly the region showing how the blocking fragment blocks the mAb, has not been elucidated by X-ray crystallography or cryo-EM. To achieve maximal effect, a pro-Ab must have high antigen-blocking and protease-restoring efficiencies, but the unclear structure limits its further optimization. Here, we utilized molecular dynamics (MD) simulations to study the dynamic structures of a hinge-based Ab lock pro-Ab, pro-Nivolumab, and validated the simulated structures with small- and wide-angle X-ray scattering (SWAXS). The MD results were closely consistent with SWAXS data (χ2 best-fit = 1.845, χ2 allMD = 3.080). The further analysis shows a pronounced flexibility of the Ab lock (root-mean-square deviation = 10.90 Å), yet it still masks the important antigen-binding residues by 57.3%-88.4%, explaining its 250-folded antigen-blocking efficiency. The introduced protease accessible surface area method affirmed better protease efficiency for light chain (33.03 Å2) over heavy chain (5.06 Å2), which aligns with the experiments. Overall, we developed MD-SWAXS validation method to study the dynamics of flexible blocking segments and introduced methodologies to estimate their antigen-blocking and protease-restoring efficiencies, which would potentially be advancing the clinical applications of any spatial hindrance-based pro-Ab.

Computational method for designing vaccines applied to virus-like particles (VLPs) as epitope carriers

Nonenveloped virus-like particles (VLPs) are self-assembled oligomeric structures composed of one or more proteins that originate from diverse viruses. Because these VLPs have similar antigenicity to the parental virus, they are successfully used as vaccines against cognate virus infection. Furthermore, after foreign antigenic sequences are inserted in their protein components (chimVLPs), some VLPs are also amenable to producing vaccines against pathogens other than the virus it originates from (these VLPs are named platform or epitope carrier). Designing chimVLP vaccines is challenging because the immunogenic response must be oriented against a given antigen without altering stimulant properties inherent to the VLP. An important step in this process is choosing the location of the sequence modifications because this must be performed without compromising the assembly and stability of the original VLP. Currently, many immunogenic data and computational tools can help guide the design of chimVLPs, thus reducing experimental costs and work. In this study, we analyze the structure of a novel VLP that originate from an insect virus and describe the putative regions of its three structural proteins amenable to insertion. For this purpose, we employed molecular dynamics (MD) simulations to assess chimVLP stability by comparing mutated and wild-type (WT) VLP protein trajectories. We applied this procedure to design a chimVLP that can serve as a prophylactic vaccine against the SARS-CoV-2 virus. The methodology described in this work is generally applicable for VLP-based vaccine development.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

B-cell epitope discovery: The first protein flexibility-based algorithm-Zika virus conserved epitope demonstration

Antibody-antigen interaction-at antigenic local environments called B-cell epitopes-is a prominent mechanism for neutralization of infection. Effective mimicry, and display, of B-cell epitopes is key to vaccine design. Here, a physical approach is evaluated for the discovery of epitopes which evolve slowly over closely related pathogens (conserved epitopes). The approach is 1) protein flexibility-based and 2) demonstrated with clinically relevant enveloped viruses, simulated via molecular dynamics. The approach is validated against 1) seven structurally characterized enveloped virus epitopes which evolved the least (out of thirty-nine enveloped virus-antibody structures), 2) two structurally characterized non-enveloped virus epitopes which evolved slowly (out of eight non-enveloped virus-antibody structures), and 3) eight preexisting epitope and peptide discovery algorithms. Rationale for a new benchmarking scheme is presented. A data-driven epitope clustering algorithm is introduced. The prediction of five Zika virus epitopes (for future exploration on recombinant vaccine technologies) is demonstrated. For the first time, protein flexibility is shown to outperform solvent accessible surface area as an epitope discovery metric.

Computational Design of a Multi-Epitope Vaccine Against Porphyromonas gingivalis

Porphyromonas gingivalis is a Gram-negative pathogenic bacterium associated with chronic periodontitis. The development of a chimeric peptide-based vaccine targeting this pathogen could be highly beneficial in preventing oral bone loss as well as other severe gum diseases. We applied a computational framework to design a multi-epitope-based vaccine candidate against P. gingivalis. The vaccine comprises epitopes from subunit proteins prioritized from the P. gingivalis reference strain (P. gingivalis ATCC 33277) using several reported vaccine properties. Protein-based subunit vaccines were prioritized through genomics techniques. Epitope prediction was performed using immunoinformatic servers and tools. Molecular modeling approaches were used to build a putative three-dimensional structure of the vaccine to understand its interactions with host immune cells through biophysical techniques such as molecular docking simulation studies and binding free energy methods. Genome subtraction identified 18 vaccine targets: six outer-membrane, nine cytoplasmic membrane-, one periplasmic, and two extracellular proteins. These proteins passed different vaccine checks required for the successful development of a vaccine candidate. The shortlisted proteins were subjected to immunoinformatic analysis to map B-cell derived T-cell epitopes, and antigenic, water-soluble, non-toxic, and good binders of DRB1*0101 were selected. The epitopes were then modeled into a multi-epitope peptide vaccine construct (linked epitopes plus adjuvant) to enhance immunogenicity and effectively engage both innate and adaptive immunity. Further, the molecular docking approach was used to determine the binding conformation of the vaccine to TLR2 innate immune receptor. Molecular dynamics simulations and binding free energy calculations of the vaccine–TLR2 complex were performed to highlight key intermolecular binding energies. Findings of this study will be useful for vaccine developers to design an effective vaccine for chronic periodontitis pathogens, specifically P. gingivalis.

Fighting viruses with computers, right now

The synergistic conjunction of various technological revolutions with the accumulated knowledge and workflows is rapidly transforming several scientific fields. Particularly, Virology can now feed from accurate physical models, polished computational tools, and massive computational power to readily integrate high-resolution structures into biological representations of unprecedented detail. That preparedness allows for the first time to get crucial information for vaccine and drug design from in-silico experiments against emerging pathogens of worldwide concern at relevant action windows. The present work reviews some of the main milestones leading to these breakthroughs in Computational Virology, providing an outlook for future developments in capacity building and accessibility to computational resources.

High-throughput computational pipeline for 3-D structure preparation and in silico protein surface property screening: A case study on HBcAg dimer structures

Knowledge-based experimental design can aid biopharmaceutical high-throughput screening (HTS) experiments needed to identify critical manufacturability parameters. Prior knowledge can be obtained via computational methods such as protein property extraction from 3-D protein structures. This study presents a high-throughput 3-D structure preparation and refinement pipeline that supports structure screenings with an automated and data-dependent workflow. As a case study, three chimeric virus-like particle (VLP) building blocks, hepatitis B core antigen (HBcAg) dimers, were constructed. Molecular dynamics (MD) refinement quality, speed, stability, and correlation to zeta potential data was evaluated using different MD simulation settings. Settings included 2 force fields (YASARA2 and AMBER03) and 2 pKa computation methods (YASARA and H++). MD simulations contained a data-dependent termination via identification of a 2 ns Window of Stability, which was also used for robust descriptor extraction. MD simulation with YASARA2, independent of pKa computation method, was found to be most stable and computationally efficient. These settings resulted in a fast refinement (6.6-37.5 h), a good structure quality (-1.17--1.13) and a strong linear dependence between dimer surface charge and complete chimeric HBcAg VLP zeta potential. These results indicate the computational pipeline's applicability for early-stage candidate assessment and design optimization of HTS manufacturability or formulability experiments.

Native State Organization of Outer Membrane Porins Unraveled by HDx-MS

Hydrogen-deuterium exchange (HDx) associated with mass spectrometry (MS) is emerging as a powerful tool to provide conformational information about membrane proteins. Unfortunately, as for X-ray diffraction and NMR, HDx performed on reconstituted in vitro systems might not always reflect the in vivo environment. Outer-membrane vesicles naturally released by Escherichia coli were used to carry out analysis of native OmpF through HDx-MS. A new protocol compatible with HDx analysis that avoids hindrance from the lipid contents was setup. The extent of deuterium incorporation was in good agreement with the X-ray diffraction data of OmpF as the buried β-barrels incorporated a low amount of deuterium, whereas the internal loop L3 and the external loops incorporated a higher amount of deuterium. Moreover, the kinetics of incorporation clearly highlights that peptides segregate well in two distinct groups based exclusively on a trimeric organization of OmpF in the membrane: peptides presenting fast kinetics of labeling are facing the complex surrounding environment, whereas those presenting slow kinetics are located in the buried core of the trimer. The data show that HDx-MS applied to a complex biological system is able to reveal solvent accessibility and spatial arrangement of an integral outer-membrane protein complex.

Multiscale Modeling in the Clinic: Drug Design and Development

A wide range of length and time scales are relevant to pharmacology, especially in drug development, drug design and drug delivery. Therefore, multiscale computational modeling and simulation methods and paradigms that advance the linkage of phenomena occurring at these multiple scales have become increasingly important. Multiscale approaches present in silico opportunities to advance laboratory research to bedside clinical applications in pharmaceuticals research. This is achievable through the capability of modeling to reveal phenomena occurring across multiple spatial and temporal scales, which are not otherwise readily accessible to experimentation. The resultant models, when validated, are capable of making testable predictions to guide drug design and delivery. In this review we describe the goals, methods, and opportunities of multiscale modeling in drug design and development. We demonstrate the impact of multiple scales of modeling in this field. We indicate the common mathematical and computational techniques employed for multiscale modeling approaches used in pharmacometric and systems pharmacology models in drug development and present several examples illustrating the current state-of-the-art models for (1) excitable systems and applications in cardiac disease; (2) stem cell driven complex biosystems; (3) nanoparticle delivery, with applications to angiogenesis and cancer therapy; (4) host-pathogen interactions and their use in metabolic disorders, inflammation and sepsis; and (5) computer-aided design of nanomedical systems. We conclude with a focus on barriers to successful clinical translation of drug development, drug design and drug delivery multiscale models.

Engineering virus-like particles as vaccine platforms

Virus-like particles (VLPs) have been utilized as vaccine platforms to increase the immunogenicity of heterologous antigens. A variety of diverse VLP types can serve as vaccine platforms, and research has focused on engineering VLPs to improve their efficacy as vaccines, enhance their stability, and allow for more versatile display of antigens. Here, we review selected VLP vaccine platforms, highlight efforts to improve these platforms through structure-informed rational design, and point to areas of future research that will assist in these efforts.

Lessons learned from successful human vaccines: Delineating key epitopes by dissecting the capsid proteins

Recombinant VLP-based vaccines have been successfully used against 3 diseases caused by viral infections: Hepatitis B, cervical cancer and hepatitis E. The VLP approach is attracting increasing attention in vaccine design and development for human and veterinary use. This review summarizes the clinically relevant epitopes on the VLP antigens in successful human vaccines. These virion-like epitopes, which can be delineated with molecular biology, cryo-electron microscopy and x-ray crystallographic methods, are the prerequisites for these efficacious vaccines to elicit functional antibodies. The critical epitopes and key factors influencing these epitopes are discussed for the HEV, HPV and HBV vaccines. A pentamer (for HPV) or a dimer (for HEV and HBV), rather than a monomer, is the basic building block harboring critical epitopes for the assembly of VLP antigen. The processing and formulation of VLP-based vaccines need to be developed to promote the formation and stabilization of these epitopes in the recombinant antigens. Delineating the critical epitopes is essential for antigen design in the early phase of vaccine development and for critical quality attribute analysis in the commercial phase of vaccine manufacturing.

Broad spectrum assessment of the epitope fluctuation--Immunogenicity hypothesis

Prediction of immunogenicity is a substantial barrier in vaccine design. Here, a molecular dynamics approach to assessing the immunogenicity of nanoparticles based on structure is presented. Molecular properties of epitopes on nonenveloped viral particles are quantified via a set of metrics. One such metric, epitope fluctuation (and implied flexibility), is shown to be inversely correlated with immunogenicity for each of a broad spectrum of nonenveloped viruses. The molecular metrics and experimentally determined immunogenicities for these viruses are archived in the open-source vaccine computer-aided design database. Results indicate the promise of computer-aided vaccine design to bring greater efficiency to traditional lab-based vaccine discovery approaches.

Modeling Viral Capsid Assembly

I present a review of the theoretical and computational methodologies that have been used to model the assembly of viral capsids. I discuss the capabilities and limitations of approaches ranging from equilibrium continuum theories to molecular dynamics simulations, and I give an overview of some of the important conclusions about virus assembly that have resulted from these modeling efforts. Topics include the assembly of empty viral shells, assembly around single-stranded nucleic acids to form viral particles, and assembly around synthetic polymers or charged nanoparticles for nanotechnology or biomedical applications. I present some examples in which modeling efforts have promoted experimental breakthroughs, as well as directions in which the connection between modeling and experiment can be strengthened.

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Differential structural dynamics and antigenicity of two similar influenza H5N1 virus HA-specific HLA-A*0201-restricted CLT epitopes

The free RI-10 and KI-10 peptides showed bent rather than extended conformations as binding in the cleft of the HLA-A*0201 molecule, and KI-10–HLA-A*0201 processed a much higher flexibility than RI-10–HLA-A*0201.

Bioengineering virus-like particles as vaccines

Virus-like particle (VLP) technology seeks to harness the optimally tuned immunostimulatory properties of natural viruses while omitting the infectious trait. VLPs that assemble from a single protein have been shown to be safe and highly efficacious in humans, and highly profitable. VLPs emerging from basic research possess varying levels of complexity and comprise single or multiple proteins, with or without a lipid membrane. Complex VLP assembly is traditionally orchestrated within cells using black-box approaches, which are appropriate when knowledge and control over assembly are limited. Recovery challenges including those of adherent and intracellular contaminants must then be addressed. Recent commercial VLPs variously incorporate steps that include VLP in vitro assembly to address these problems robustly, but at the expense of process complexity. Increasing research activity and translation opportunity necessitate bioengineering advances and new bioprocessing modalities for efficient and cost-effective production of VLPs. Emerging approaches are necessarily multi-scale and multi-disciplinary, encompassing diverse fields from computational design of molecules to new macro-scale purification materials. In this review, we highlight historical and emerging VLP vaccine approaches. We overview approaches that seek to specifically engineer a desirable immune response through modular VLP design, and those that seek to improve bioprocess efficiency through inhibition of intracellular assembly to allow optimal use of existing purification technologies prior to cell-free VLP assembly. Greater understanding of VLP assembly and increased interdisciplinary activity will see enormous progress in VLP technology over the coming decade, driven by clear translational opportunity.

Epitope fluctuations in the human papillomavirus are under dynamic allosteric control: a computational evaluation of a new vaccine design strategy

The dynamic properties of the capsid of the human papillomavirus (HPV) type 16 were examined using classical molecular dynamics simulations. By systematically comparing the structural fluctuations of the capsid protein, a strong dynamic allosteric connection between the epitope containing loops and the h4 helix located more than 50 Å away is identified, which was not recognized thus far. Computer simulations show that restricting the structural fluctuations of the h4 helix is key to rigidifying the epitopes, which is thought to be required for eliciting a proper immune response. The allostery identified in the components of the HPV is nonclassical because the mean structure of the epitope carrying loops remains unchanged, but as a result of allosteric effect the structural fluctuations are altered significantly, which in turn changes the biochemical reactivity profile of the epitopes. Exploiting this novel insight, a new vaccine design strategy is proposed wherein a relatively small virus capsid fragment is deposited on a silica nanoparticle in such a way that the fluctuations of the h4 helix are suppressed. The structural and dynamic properties of the epitope carrying loops on this hybrid nanoparticle match the characteristics of epitopes found on the full virus-like particle precisely, suggesting that these nanoparticles may serve as potent, cost-effective, and safe alternatives to traditionally developed vaccines. The structural and dynamic properties of the hybrid nanoparticle are examined in detail to establish the general concepts of the proposed new design.

Epitope engineering and molecular metrics of immunogenicity: a computational approach to VLP-based vaccine design

Developing antiviral vaccines is increasingly challenging due to associated time and cost of production as well as emerging drug-resistant strains. A computer-aided vaccine design strategy is presented that could greatly accelerate the discovery process and yield vaccines with high immunogenicity and thermal stability. Our strategy is based on foreign viral epitopes engineered onto well-established virus-like particles (VLPs) and demonstrates that such constructs present similar affinity for antibodies as does a native virus. This binding affinity serves as one molecular metric of immunogenicity. As a demonstration, we engineered a preS1 epitope of hepatitis B virus (HBV) onto the EF loop of human papillomavirus VLP (HPV-VLP). HBV-associated HzKR127 antibody displayed binding affinity for this structure at distances and strengths similar to those for the complex of the antibody with the full HBV (PDBID: 2EH8). This antibody binding affinity assessment, along with other molecular immunogenicity metrics, could be a key component of a computer-aided vaccine design strategy.

Variational methods for time-dependent classical many-particle systems

A variational method for the classical Liouville equation is introduced that facilitates the development of theories for non-equilibrium classical systems. The method is based on the introduction of a complex-valued auxiliary quantity Ψ that is related to the classical position-momentum probability density ρ via ρ = Ψ*Ψ. A functional of Ψ is developed whose extrema imply that ρ satisfies the Liouville equation. Multiscale methods are used to develop trial functions to be optimized by the variational principle. The present variational principle with multiscale trial functions can capture both the microscopic and the coarse-grained descriptions, thereby yielding theories that account for the two way exchange of information across multiple scales in space and time. Equations of the Smoluchowski form for the coarse-grained state probability density are obtained. Constraints on the initial state of the N-particle probability density for which the aforementioned equation is closed and conserves probability are presented. The methodology has applicability to a wide range of systems including macromolecular assemblies, ionic liquids, and nanoparticles.

Patent Highlights

A snapshot of recent key developments in the patent literature of relevance to the advancement of pharmaceutical and medical R&D