High purity and recovery of native filamentous hemagglutinin (FHA) from Bordetella pertussis using affinity chromatography.

Filamentous hemagglutinin (FHA) is a critical adhesion molecule produced by Bordetella pertussis (BP), the causative agent of highly contagious respiratory infection known as whooping cough. FHA plays a pivotal role in the pathogenesis of whooping cough and is a key component of acellular pertussis vaccines (aPV). However, conventional purification methods for FHA often involve labor-intensive processes and result in low purity and recovery rates. Therefore, this study explores the use of monoclonal and polyclonal antibodies as specific tools to achieve highly pure and efficient FHA purification. To generate FHA-specific antibodies, polyclonal antibodies were produced by immunizing sheep and monoclonal antibodies (MAbs) were generated by immunizing mice with recombinant and native FHA. The MAbs were selected based on affinity, isotypes, and specificity, which were assessed through ELISA and Western blot assays. Two immunoaffinity columns, one monoclonal and one polyclonal, were prepared for FHA antigen purification. The purity and recovery rates of these purifications were determined using ELISA, SDS-PAGE, and immunoblotting. Furthermore, the MAbs were employed to develop an ELISA assay for FHA antigen concentration determination. The study's findings revealed that immunoaffinity column-based purification of FHA resulted in a highly pure antigen with recovery rates of approximately 57% ± 6.5% and 59% ± 7.9% for monoclonal and polyclonal columns, respectively. Additionally, the developed ELISA exhibited appropriate reactivity for determining FHA antigen concentration. This research demonstrates that affinity chromatography is a viable and advantageous method for purifying FHA, offering superior purity and recovery rates compared to traditional techniques. This approach provides a practical alternative for FHA purification in the context of aPV development.

Secretion of Pertussis Toxin from Bordetella pertussis.

Production and secretion of pertussis toxin (PT) is essential for the virulence of Bordetella pertussis. Due to the large oligomeric structure of PT, transport of the toxin across bacterial membrane barriers represents a significant hurdle that the bacteria must overcome in order to maintain pathogenicity. During the secretion process, PT undergoes a two-step transport process. The first step involves transport of the individual polypeptide chains of PT across the inner membrane utilizing a generalized secretion pathway, most likely the bacterial Sec system. The second step involves the use of a specialized apparatus to transport the toxin across the outer membrane of the bacterial cell. This apparatus, which has been termed the Ptl transporter and which is unique to the PT secretion pathway, is a member of the type IV family of bacterial transporters. Here, the current understanding of the PT secretion process is reviewed including a description of the Ptl proteins that assemble to form the transporter, the general structure of type IV transporters, the known similarities and differences between canonical type IV substrate transport and Ptl-mediated transport of PT, as well as the known sequence of events in the assembly and secretion of PT.

Molecular Basis for Bordetella pertussis Interference with Complement, Coagulation, Fibrinolytic, and Contact Activation Systems: the Cryo-EM Structure of the Vag8-C1 Inhibitor Complex.

Complement, contact activation, coagulation, and fibrinolysis are serum protein cascades that need strict regulation to maintain human health. Serum glycoprotein, a C1 inhibitor (C1-INH), is a key regulator (inhibitor) of serine proteases of all the above-mentioned pathways. Recently, an autotransporter protein, virulence-associated gene 8 (Vag8), produced by the whooping cough pathogen, Bordetella pertussis, was shown to bind to C1-INH and interfere with its function. Here, we present the structure of the Vag8-C1-INH complex determined using cryo-electron microscopy at a 3.6-Å resolution. The structure shows a unique mechanism of C1-INH inhibition not employed by other pathogens, where Vag8 sequesters the reactive center loop of C1-INH, preventing its interaction with the target proteases.IMPORTANCE The structure of a 10-kDa protein complex is one of the smallest to be determined using cryo-electron microscopy at high resolution. The structure reveals that C1-INH is sequestered in an inactivated state by burial of the reactive center loop in Vag8. By so doing, the bacterium is able to simultaneously perturb the many pathways regulated by C1-INH. Virulence mechanisms such as the one described here assume more importance given the emerging evidence about dysregulation of contact activation, coagulation, and fibrinolysis leading to COVID-19 pneumonia.

Development of novel isatin-nicotinohydrazide hybrids with potent activity against susceptible/resistant Mycobacterium tuberculosis and bronchitis causing-bacteria.

Joining the global fight against Tuberculosis, the world's most deadly infectious disease, herein we present the design and synthesis of novel isatin-nicotinohydrazide hybrids (5a-m and 9a-c) as promising anti-tubercular and antibacterial agents. The anti-tubercular activity of the target hybrids was evaluated against drug-susceptible M. tuberculosis strain (ATCC 27294) where hybrids 5d, 5g and 5h were found to be as potent as INH with MIC = 0.24 µg/mL, also the activity was evaluated against Isoniazid/Streptomycin resistant M. tuberculosis (ATCC 35823) where compounds 5g and 5h showed excellent activity (MIC = 3.9 µg/mL). Moreover, the target hybrids were examined against six bronchitis causing-bacteria. Most derivatives exhibited excellent antibacterial activity. K. pneumonia emerged as the most sensitive strain with MIC range: 0.49-7.81 µg/mL. Furthermore, a molecular docking study has proposed DprE1 as a probable enzymatic target for herein reported isatin-nicotinohydrazide hybrids, and explored the binding interactions within the vicinity of DprE1 active site.

Properties of protein unfolded states suggest broad selection for expanded conformational ensembles.

Much attention is being paid to conformational biases in the ensembles of intrinsically disordered proteins. However, it is currently unknown whether or how conformational biases within the disordered ensembles of foldable proteins affect function in vivo. Recently, we demonstrated that water can be a good solvent for unfolded polypeptide chains, even those with a hydrophobic and charged sequence composition typical of folded proteins. These results run counter to the generally accepted model that protein folding begins with hydrophobicity-driven chain collapse. Here we investigate what other features, beyond amino acid composition, govern chain collapse. We found that local clustering of hydrophobic and/or charged residues leads to significant collapse of the unfolded ensemble of pertactin, a secreted autotransporter virulence protein from Bordetella pertussis, as measured by small angle X-ray scattering (SAXS). Sequence patterns that lead to collapse also correlate with increased intermolecular polypeptide chain association and aggregation. Crucially, sequence patterns that support an expanded conformational ensemble enhance pertactin secretion to the bacterial cell surface. Similar sequence pattern features are enriched across the large and diverse family of autotransporter virulence proteins, suggesting sequence patterns that favor an expanded conformational ensemble are under selection for efficient autotransporter protein secretion, a necessary prerequisite for virulence. More broadly, we found that sequence patterns that lead to more expanded conformational ensembles are enriched across water-soluble proteins in general, suggesting protein sequences are under selection to regulate collapse and minimize protein aggregation, in addition to their roles in stabilizing folded protein structures.

Negative charge of the AC-to-Hly linking segment modulates calcium-dependent membrane activities of Bordetella adenylate cyclase toxin.

Two distinct conformers of the adenylate cyclase toxin (CyaA) appear to accomplish its two parallel activities within target cell membrane. The translocating conformer would deliver the N-terminal adenylyl cyclase (AC) enzyme domain across plasma membrane into cytosol of cells, while the pore precursor conformer would assemble into oligomeric cation-selective pores and permeabilize cellular membrane. Both toxin activities then involve a membrane-interacting 'AC-to-Hly-linking segment' (residues 400 to 500). Here, we report the NMR structure of the corresponding CyaA411-490 polypeptide in dodecylphosphocholine micelles and show that it consists of two α-helices linked by an unrestrained loop. The N-terminal α-helix (Gly418 to His439) remained solvent accessible, while the C-terminal α-helix (His457 to Phe485) was fully enclosed within detergent micelles. CyaA411-490 weakly bound Ca2+ ions (apparent KD 2.6 mM) and permeabilized negatively charged lipid vesicles. At high concentrations (10 µM) the CyaA411-490 polypeptide formed stable conductance units in artificial lipid bilayers with applied voltage, suggesting its possible transmembrane orientation in the membrane-inserted toxin. Mutagenesis revealed that two clusters of negatively charged residues within the 'AC-to-Hly-linking segment' (Glu419 to Glu432 and Asp445 to Glu448) regulate the balance between the AC domain translocating and pore-forming capacities of CyaA in function of calcium concentration.

Structural and Functional Characterization of the Globin-Coupled Sensors of Azotobacter vinelandii and Bordetella pertussis.

Aims: Structural and functional characterization of the globin-coupled sensors (GCSs) from Azotobacter vinelandii (AvGReg) and Bordetella pertussis (BpeGReg). Results: Ultraviolet/visible and resonance Raman spectroscopies confirm the presence in AvGReg and BpeGReg of a globin domain capable of reversible gaseous ligand binding. In AvGReg, an influence of the transmitter domain on the heme proximal region of the globin domain can be seen, and k'CO is higher than for other GCSs. The O2 binding kinetics suggests the presence of an open and a closed conformation. As for BpeGReg, the fully oxygenated AvGReg show a very high diguanylate cyclase activity. The carbon monoxide rebinding to BpeGReg indicates that intra- and intermolecular interactions influence the ligand binding. The globin domains of both proteins (AvGReg globin domain and BpeGRegGb with cysteines (Cys16, 45, 114, 154) mutated to serines [BpeGReg-Gb*]) share the same GCS fold, a similar proximal but a different distal side structure. They homodimerize through a G-H helical bundle as in other GCSs. However, BpeGReg-Gb* shows also a second dimerization mode. Innovation: This article extends our knowledge on the GCS proteins and contributes to a better understanding of the GCSs role in the formation of bacterial biofilms. Conclusions:AvGReg and BpeGReg conform to the GCS family, share a similar overall structure, but they have different properties in terms of the ligand binding. In particular, AvGReg shows an open and a closed conformation that in the latter form will very tightly bind oxygen. BpeGReg has only one closed conformation. In both proteins, it is the fully oxygenated GCS form that catalyzes the production of the second messenger.

Irreversible versus repairable membrane poration: differences in permeabilization elicited by Bordetella Adenylate Cyclase Toxin and its hemolysin domain in macrophages.

Rapid plasma membrane repair in response to pore-forming toxins is crucial for cell survival, but the molecular mechanisms employed by eukaryotic nucleated cells to maintain membrane integrity and the specificities of such pathways remain poorly understood. Here, we have explored the permeabilization elicited by the Bordetella pertussis adenylate cyclase toxin, a 200-kDa protein toxin with α-helical pore-forming domain that forms pores of tunable size, and evaluated the response of target macrophages to such toxin poration. We show here that the response and the fate of target macrophages depend on toxin pore width. We find that the toxin's hemolysin moiety induces a transient membrane permeabilization by forming wide enough pores allowing Ca2+ influx into the target cell cytosol. This activates a Ca2+ -dependent cellular response involving exocytosis and endocytosis steps eliminating toxin pores and restoring membrane integrity. In contrast, the full-length native toxin, at low concentrations, forms very small pores that cause insidious perturbation of cell ion homeostasis that escapes control by the macrophage membrane repair response, eventually leading to cell death. Our data reveal that permeability to Ca2+ and ATP are key elements in the membrane repair pathway for eliminating α-helical pores of bacterial origin.

Single Molecule Force Spectroscopy Reveals the Mechanical Design Governing the Efficient Translocation of the Bacterial Toxin Protein RTX.

The efficient translocation of the bacterial toxin adenylate cyclase toxin (CyaA) from the bacterial cytosol to the extracellular environment by the type 1 secretion system (T1SS) is essential for the toxin to function. To understand the molecular features that are responsible for the efficient translocation of CyaA, here we used optical tweezers to investigate the mechanical properties and conformational dynamics of the RTX domain of CyaA at the single molecule level. Our results revealed that apo-RTX behaves like an ideal random coil. This property allows the T1SS to translocate RTX without overcoming the enthalpic resistance. In contrast, the folded holo-RTX is mechancially stable, and its folding occurs in a vectorial, cotranslocational fashion starting from its C-terminus. Moreover, our results showed that the folding of holo-RTX generates a stretching force, which can further facilitate the translocation of RTX. Our results highlight the important role played by the Ca2+-triggered folding of RTX in the translocation of RTX and provide mechanistic insights into the mechanical design that governs the efficient translocation of RTX.

Bordetella Adenylate Cyclase Toxin Inhibits Monocyte-to-Macrophage Transition and Dedifferentiates Human Alveolar Macrophages into Monocyte-like Cells.

Monocytes arriving at the site of infection differentiate into functional effector macrophages to replenish the resident sentinel cells. Bordetella pertussis, the pertussis agent, secretes an adenylate cyclase toxin-hemolysin (CyaA) that binds myeloid phagocytes through complement receptor 3 (CD11b/CD18) and swiftly delivers its adenylyl cyclase enzyme domain into phagocytes. This ablates the bactericidal capacities of phagocytes through massive and unregulated conversion of cytosolic ATP into the key signaling molecule cAMP. We show that exposure of primary human monocytes to as low a concentration as 22.5 pM CyaA, or a low (2:1) multiplicity of infection by CyaA-producing B. pertussis bacteria, blocks macrophage colony-stimulating factor (M-CSF)-driven differentiation of monocytes. CyaA-induced cAMP signaling mediated through the activity of protein kinase A (PKA) efficiently blocked expression of macrophage markers, and the monocytes exposed to 22.5 pM CyaA failed to acquire the characteristic intracellular complexity of mature macrophage cells. Neither M-CSF-induced endoplasmic reticulum (ER) expansion nor accumulation of Golgi bodies, mitochondria, or lysosomes was observed in toxin-exposed monocytes, which remained small and poorly phagocytic and lacked pseudopodia. Exposure to 22.5 pM CyaA toxin provoked loss of macrophage marker expression on in vitro differentiated macrophages, as well as on primary human alveolar macrophages, which appeared to dedifferentiate into monocyte-like cells with upregulated CD14 levels. This is the first report that terminally differentiated tissue-resident macrophage cells can be dedifferentiated in vitro The results suggest that blocking of monocyte-to-macrophage transition and/or dedifferentiation of the sentinel cells of innate immunity through cAMP-elevating toxin action may represent a novel immune evasion strategy of bacterial pathogens.IMPORTANCE Macrophages are key sentinel cells of the immune system, and, as such, they are targeted by the toxins produced by the pertussis agent Bordetella pertussis The adenylate cyclase toxin (CyaA) mediates immune evasion of B. pertussis by suspending the bactericidal activities of myeloid phagocytes. We reveal a novel mechanism of potential subversion of host immunity, where CyaA at very low (22 pM) concentrations could inhibit maturation of human monocyte precursors into the more phagocytic macrophage cells. Furthermore, exposure to low CyaA amounts has been shown to trigger dedifferentiation of mature primary human alveolar macrophages back into monocyte-like cells. This unprecedented capacity is likely to promote survival of the pathogen in the airways, both by preventing maturation of monocytes attracted to the site of infection into phagocytic macrophages and by dedifferentiation of the already airway-resident sentinel cells.

An intrinsic fluorescence method for the determination of protein concentration in vaccines containing aluminum salt adjuvants.

Determination of protein concentration in vaccines containing aluminum salt adjuvant typically necessitates desorption of the protein prior to analysis. Here we describe a method based on the intrinsic fluorescence of tyrosine and tryptophan that requires no desorption of proteins. Adjuvanted formulations of three model Bordetella pertussis antigens were excited at 280 nm and their emission spectra collected from 290 to 400 nm. Emission spectra of protein antigens in the presence of aluminum salt adjuvants were able to be detected, the effects of adjuvants on the spectra were analyzed, and linear regressions were calculated. The fluorescence method proved to be very sensitive with a limit of quantification between 0.4 and 4.4 µg/mL and limit of linearity between 100 and 200 µg/mL, across the formulations tested. The fluorescence method was found to be influenced by adjuvant presence, type of adjuvant, adjuvant concentration, buffer and pH conditions. The method also demonstrated ability to monitor the percent adsorption of antigens to the adjuvants. Furthermore, intrinsic fluorescence showed good correlation with micro-Kjeldahl elemental assay in quantifying protein concentration. Being a non-invasive, quick and sensitive method, intrinsic fluorescence has the potential to be utilized as a high throughput tool for vaccine development and conceivably implemented in-line, using in-line fluorimeters, to monitor antigen concentration during formulation processing.

Development of a thermostable spray dried outer membrane vesicle pertussis vaccine for pulmonary immunization.

Worldwide resurgence of whooping cough calls for improved, next-generation pertussis vaccines that induce broad and long-lasting immunity. A mucosal pertussis vaccine based on outer membrane vesicles (omvPV) is a promising candidate. Further, a vaccine that is stable outside the cold chain would be of substantial advantage for worldwide distribution and application. A vaccine formulated as a powder could both stabilize the vaccine as well as make it suitable for pulmonary vaccination. To that end, we developed a spray dried omvPV with improved stability compared to the liquid omvPV formulation. Spray drying did not affect the structural integrity of the omvPV. The antigenicity of Vag8, a major antigen in omvPV was diminished slightly and an altered tryptophan fluorescence indicated some changes in protein structure. However, when administered via the pulmonary route in mice after reconstitution, spray dried omvPV showed comparable immune responses and protection against challenge with live B. pertussis as liquid omvPV. Mucosal IgA and Th17 responses were established in addition to broad systemic IgG and Th1/Th17 responses, indicating the induction of an effective immunity profile. Overall, a spray dried omvPV was developed that maintained effective immunogenic properties and has an improved storage stability.

Cork-in-Bottle Occlusion of Fluoride Ion Channels by Crystallization Chaperones.

Crystallization of dual-topology fluoride (Fluc) channels requires small, soluble crystallization chaperones known as monobodies, which act as primary crystal lattice contacts. Previous structures of Flucs have been solved in the presence of monobodies that inhibit fluoride currents in single-channel electrophysiological recordings. These structures have revealed two-fold symmetric, doubly bound arrangements, with one monobody on each side of the membrane. The combined electrophysiological and structural observations raise the possibility that chaperone binding allosterically closes the channel, altering the structure from its conducting form. To address this, we identify and solve the structure with a different monobody that only partially blocks fluoride currents. The structure of the channel-monobody complex is asymmetric, with monobody bound to one side of the channel only. The channel conformation is nearly identical on the bound and uncomplexed sides, and to all previously solved structures, providing direct structural evidence that monobody binding does not induce local structural changes.

SEC-SAXS and HDX-MS: A powerful combination. The case of the calcium-binding domain of a bacterial toxin.

Small-angle X-ray scattering (SAXS) is a relatively simple experimental technique that provides information on the global conformation of macromolecules in solution, be they fully structured, partially, or extensively unfolded. Size exclusion chromatography in line with a SAXS measuring cell considerably improves the monodispersity and ideality of solutions, the two main requirements of a "good" SAXS sample. Hydrogen/deuterium exchange monitored by mass spectrometry (HDX-MS) offers a wealth of information regarding the solvent accessibility at the local (peptide) level. It constitutes a sensitive probe of local flexibility and, more generally, of structural dynamics. The combination of both approaches presented here is very powerful, as illustrated by the case of RD, a calcium-binding protein that is part of a bacterial virulence factor.

Calmodulin fishing with a structurally disordered bait triggers CyaA catalysis.

Once translocated into the cytosol of target cells, the catalytic domain (AC) of the adenylate cyclase toxin (CyaA), a major virulence factor of Bordetella pertussis, is potently activated by binding calmodulin (CaM) to produce supraphysiological levels of cAMP, inducing cell death. Using a combination of small-angle X-ray scattering (SAXS), hydrogen/deuterium exchange mass spectrometry (HDX-MS), and synchrotron radiation circular dichroism (SR-CD), we show that, in the absence of CaM, AC exhibits significant structural disorder, and a 75-residue-long stretch within AC undergoes a disorder-to-order transition upon CaM binding. Beyond this local folding, CaM binding induces long-range allosteric effects that stabilize the distant catalytic site, whilst preserving catalytic loop flexibility. We propose that the high enzymatic activity of AC is due to a tight balance between the CaM-induced decrease of structural flexibility around the catalytic site and the preservation of catalytic loop flexibility, allowing for fast substrate binding and product release. The CaM-induced dampening of AC conformational disorder is likely relevant to other CaM-activated enzymes.

Identification and characterization of a brilliant yellow pigment produced by Bordetella pertussis.

Culture supernatants of Bordetella pertussis are a brilliant yellow; however, the structure and biological role of the responsible pigment have not been investigated. In this study, a brilliant yellow-colored fraction was extracted from culture supernatants of B. pertussis and analyzed by HPLC. UV-visible spectral analysis and mass spectrometry identified the brilliant yellow pigment as riboflavin. Riboflavin production was high in lag and early log phases and riboflavin was found to enhance growth of B. pertussis in low-density cultures. Riboflavin production is not regulated by the BvgAS system. In addition, it was found that other Bordetella species, such as B. parapertussis, B. holmesii and B. bronchiseptica, also release riboflavin into their culture supernatants. This is the first report that B. pertussis secrets riboflavin to the extracellular space and that riboflavin may promote its growth. The mechanism may be associated with pathogenesis of B. pertussis.

Bordetella holmesii: Lipid A Structures and Corresponding Genomic Sequences Comparison in Three Clinical Isolates and the Reference Strain ATCC 51541.

Bordetella holmesii can cause invasive infections but can also be isolated from the respiratory tract of patients with whooping-cough like symptoms. For the first time, we describe the lipid A structure of B. holmesii reference strain ATCC 51541 (alias NCTC12912 or CIP104394) and those of three French B. holmesii clinical isolates originating from blood (Bho1) or from respiratory samples (FR4020 and FR4101). They were investigated using chemical analyses, gas chromatography-mass spectrometry (GC-MS), and matrix-assisted laser desorption ionization-mass spectrometry (MALDI-MS). The analyses revealed a common bisphosphorylated ß-(1→6)-linked d-glucosamine disaccharide with hydroxytetradecanoic acid in amide linkages. Similar to B. avium, B. hinzii and B. trematum lipids A, the hydroxytetradecanoic acid at the C-2' position are carrying in secondary linkage a 2-hydroxytetradecanoic acid residue resulting of post-traductional biosynthesis modifications. The three clinical isolates displayed characteristic structural traits compared to the ATCC 51541 reference strain: the lipid A phosphate groups are more or less modified with glucosamine in the isolates and reference strain, but the presence of 10:0(3-OH) is only observed in the isolates. This trait was only described in B. pertussis and B. parapertussis strains, as well as in B. petrii isolates by the past. The genetic bases for most of the key structural elements of lipid A were analyzed and supported the structural data.

Negatively charged residues of the segment linking the enzyme and cytolysin moieties restrict the membrane-permeabilizing capacity of adenylate cyclase toxin.

The whooping cough agent, Bordetella pertussis, secretes an adenylate cyclase toxin-hemolysin (CyaA) that plays a crucial role in host respiratory tract colonization. CyaA targets CR3-expressing cells and disrupts their bactericidal functions by delivering into their cytosol an adenylate cyclase enzyme that converts intracellular ATP to cAMP. In parallel, the hydrophobic domain of CyaA forms cation-selective pores that permeabilize cell membrane. The invasive AC and pore-forming domains of CyaA are linked by a segment that is unique in the RTX cytolysin family. We used mass spectrometry and circular dichroism to show that the linker segment forms α-helical structures that penetrate into lipid bilayer. Replacement of the positively charged arginine residues, proposed to be involved in target membrane destabilization by the linker segment, reduced the capacity of the toxin to translocate the AC domain across cell membrane. Substitutions of negatively charged residues then revealed that two clusters of negative charges within the linker segment control the size and the propensity of CyaA pore formation, thereby restricting the cell-permeabilizing capacity of CyaA. The 'AC to Hly-linking segment' thus appears to account for the smaller size and modest cell-permeabilizing capacity of CyaA pores, as compared to typical RTX hemolysins.

Calcium-Driven Folding of RTX Domain ß-Rolls Ratchets Translocation of RTX Proteins through Type I Secretion Ducts.

Calcium-binding RTX proteins are equipped with C-terminal secretion signals and translocate from the Ca(2+)-depleted cytosol of Gram-negative bacteria directly into the Ca(2+)-rich external milieu, passing through the "channel-tunnel" ducts of type I secretion systems (T1SSs). Using Bordetella pertussis adenylate cyclase toxin, we solved the structure of an essential C-terminal assembly that caps the RTX domains of RTX family leukotoxins. This is shown to scaffold directional Ca(2+)-dependent folding of the carboxy-proximal RTX repeat blocks into ß-rolls. The resulting intramolecular Brownian ratchets then prevent backsliding of translocating RTX proteins in the T1SS conduits and thereby accelerate excretion of very large RTX leukotoxins from bacterial cells by a vectorial "push-ratchet" mechanism. Successive Ca(2+)-dependent and cosecretional acquisition of a functional RTX toxin structure in the course of T1SS-mediated translocation, through RTX domain folding from the C-terminal cap toward the N terminus, sets a paradigm that opens for design of virulence inhibitors of major pathogens.

Implicit Time Integration for Multiscale Molecular Dynamics Using Transcendental Padé Approximants.

Molecular dynamics systems evolve through the interplay of collective and localized disturbances. As a practical consequence, there is a restriction on the time step imposed by the broad spectrum of time scales involved. To resolve this restriction, multiscale factorization was introduced for molecular dynamics as a method that exploits the separation of time scales by coevolving the coarse-grained and atom-resolved states via Trotter factorization. Developing a stable time-marching scheme for this coevolution, however, is challenging because the coarse-grained dynamical equations depend on the microstate; therefore, these equations cannot be expressed in closed form. The objective of this paper is to develop an implicit time integration scheme for multiscale simulation of large systems over long periods of time and with high accuracy. The scheme uses Padé approximants to account for both the stochastic and deterministic features of the coarse-grained dynamics. The method is demonstrated for a protein either undergoing a conformational change or migrating under the influence of an external force. The method shows promise in accelerating multiscale molecular dynamics without a loss of atomic precision or the need to conjecture the form of coarse-grained governing equations.