rowan-python 3.2.1__py3-none-any.whl → 3.2.3__py3-none-any.whl

rowan/__init__.py

@@ -15,6 +15,7 @@ from stjames import (
     Constraint,
     ConstraintType,
     ConformerClusteringSettings,
+    ConformerProperties,
     ConformerGenSettingsUnion,
     Correction,
     DNASequence,
@@ -67,6 +68,7 @@ from stjames.engine_compatibility import (
     METHOD_ENGINES,
     get_supported_corrections,
 )
+from stjames.pdb import ResidueRef
 from stjames.optimization.band_method import NEBSettings
 from stjames.optimization.interpolation import Interpolation
 from stjames.optimization.string_method import FSMSettings, StringMethodSettings
@@ -82,5 +84,6 @@ from .types import ProteinUUID, RdkitMol, SMILES, StJamesMolecule, StructureInpu
 from .workflows import *
 from .project import *
 from .protein import *
+from .spectra import HARTREE_TO_EV, broaden_spectrum
 from .user import *
 from .utils import *

rowan/project.py

@@ -197,8 +197,15 @@ def default_project() -> Project:
 
     Raises:
         httpx.HTTPStatusError: API request fails
+        ValueError: user has no default project in the API key's organization
     """
     with api_client() as client:
         response = client.get("/user/me/default_project")
         response.raise_for_status()
-        return Project(**response.json())
+        project = response.json()
+    if project is None:
+        raise ValueError(
+            "You have no default project in this organization. "
+            "Set a default project, or pass a folder or project."
+        )
+    return Project(**project)

rowan/spectra.py

@@ -0,0 +1,68 @@
+"""Broaden stick spectra into continuous line shapes."""
+
+import math
+from collections.abc import Sequence
+from typing import Literal
+
+HARTREE_TO_EV = 27.211386245988
+
+_GAUSSIAN_FWHM_TO_SIGMA = 1 / (2 * math.sqrt(2 * math.log(2)))
+
+
+def broaden_spectrum(
+    positions: Sequence[float],
+    intensities: Sequence[float],
+    grid: Sequence[float],
+    fwhm: float,
+    lineshape: Literal["gaussian", "lorentzian"] = "gaussian",
+) -> list[float]:
+    """Sum line shapes centered on each stick, each with area equal to its intensity.
+
+    Use it on the stick spectra of `ConformerSearchResult.ensemble_properties`, for
+    example excitation energies with oscillator or rotatory strengths for UV-vis or ECD,
+    or vibrational frequencies with IR intensities. `grid` and `fwhm` share the units of
+    `positions`; convert excitation energies from Hartree with `HARTREE_TO_EV` first to
+    broaden in eV.
+
+    Args:
+        positions: line positions
+        intensities: line intensities, parallel to `positions`
+        grid: points at which to evaluate the spectrum
+        fwhm: full width at half maximum of each line
+        lineshape: line shape of each stick
+
+    Returns:
+        spectrum at each point of `grid`
+
+    Raises:
+        ValueError: `fwhm` is not positive or `intensities` differs in length from `positions`
+
+    >>> [round(y, 4) for y in broaden_spectrum([0.0], [1.0], [0.0, 0.5], fwhm=1.0)]
+    [0.9394, 0.4697]
+    >>> [round(y, 4) for y in broaden_spectrum([0.0], [2.0], [0.0], 1.0, "lorentzian")]
+    [1.2732]
+    """
+    if fwhm <= 0:
+        raise ValueError(f"fwhm must be positive, got {fwhm}")
+    if len(positions) != len(intensities):
+        raise ValueError(
+            f"positions ({len(positions)}) and intensities ({len(intensities)}) "
+            "must have the same length"
+        )
+
+    if lineshape == "gaussian":
+        sigma = fwhm * _GAUSSIAN_FWHM_TO_SIGMA
+        norm = 1 / (sigma * math.sqrt(2 * math.pi))
+
+        def shape(dx: float) -> float:
+            return norm * math.exp(-0.5 * (dx / sigma) ** 2)
+
+    else:
+        gamma = fwhm / 2
+
+        def shape(dx: float) -> float:
+            return gamma / (math.pi * (dx**2 + gamma**2))
+
+    return [
+        sum(i * shape(x - p) for p, i in zip(positions, intensities, strict=True)) for x in grid
+    ]

rowan/utils.py

@@ -4,7 +4,7 @@ from contextvars import ContextVar
 from dataclasses import dataclass
 from pathlib import Path
 from tempfile import NamedTemporaryFile
-from typing import Any, Generator, Literal
+from typing import Any, Callable, ContextManager, Generator, Literal
 
 import httpx
 import stjames
@@ -167,24 +167,14 @@ def api_client() -> Generator[httpx.Client, None, None]:
         yield client
 
 
-def download_file(
-    path: Path,
-    method: Literal["GET", "POST"],
-    url: str,
-    *,
-    params: dict[str, Any] | None = None,
-    json: Any = None,
-) -> Path:
-    """Stream an authenticated API response into one local file atomically."""
+def _stream_to_file(path: Path, stream: Callable[[], ContextManager[httpx.Response]]) -> Path:
+    """Write a streamed response into one local file atomically."""
     path.parent.mkdir(parents=True, exist_ok=True)
     temporary_path: Path | None = None
     try:
         with NamedTemporaryFile(dir=path.parent, prefix=f".{path.name}.", delete=False) as output:
             temporary_path = Path(output.name)
-            with (
-                api_client() as client,
-                client.stream(method, url, params=params, json=json) as response,
-            ):
+            with stream() as response:
                 response.raise_for_status()
                 for chunk in response.iter_bytes():
                     output.write(chunk)
@@ -196,6 +186,65 @@ def download_file(
     return path
 
 
+def download_file(
+    path: Path,
+    method: Literal["GET", "POST"],
+    url: str,
+    *,
+    params: dict[str, Any] | None = None,
+    json: Any = None,
+) -> Path:
+    """Stream an authenticated API response into one local file atomically."""
+
+    @contextmanager
+    def stream() -> Generator[httpx.Response, None, None]:
+        with (
+            api_client() as client,
+            client.stream(method, url, params=params, json=json) as response,
+        ):
+            yield response
+
+    return _stream_to_file(path, stream)
+
+
+def download_presigned_files(
+    directory: Path,
+    url: str,
+    *,
+    params: dict[str, Any] | None = None,
+    json: Any = None,
+) -> list[Path]:
+    """Fetch presigned download links from the API and stream each file into `directory`.
+
+    Links are short-lived, so an expired one is refreshed once before giving up.
+    """
+
+    def fetch_links() -> list[dict[str, str]]:
+        with api_client() as client:
+            return client.post(url, params=params, json=json).json()
+
+    def download(link: dict[str, str]) -> Path:
+        # Plain client: the API key must never be sent to S3.
+        @contextmanager
+        def stream() -> Generator[httpx.Response, None, None]:
+            with httpx.Client(timeout=httpx.Timeout(120, read=None)) as client:
+                with client.stream("GET", link["url"]) as response:
+                    yield response
+
+        return _stream_to_file(directory / Path(link["filename"]).name, stream)
+
+    links = fetch_links()
+    paths = []
+    for i, link in enumerate(links):
+        try:
+            paths.append(download(link))
+        except httpx.HTTPStatusError as error:
+            if error.response.status_code != 403:
+                raise
+            paths.append(download(fetch_links()[i]))
+    return paths
+
+
 ATOMIC_NUMBER_TO_ATOMIC_SYMBOL = {
     "1": "H",
     "2": "He",

rowan/workflows/_molecular_dynamics.py

@@ -8,7 +8,7 @@ from typing import Literal, Protocol
 from stjames.pdb import PDB
 
 from rowan.protein import Protein, retrieve_protein
-from rowan.utils import api_client, download_file
+from rowan.utils import api_client, download_file, download_presigned_files
 
 from .base import WorkflowResult
 
@@ -222,7 +222,9 @@ class _MolecularDynamicsResult(WorkflowResult):
         name: str | None = None,
         path: Path | str | None = None,
     ) -> Path:
-        """Download DCD trajectory files for selected replicates.
+        """Download DCD trajectory files for selected replicates as one `.tar.gz` archive.
+
+        Long trajectories can time out; use `download_trajectory_files` for those.
 
         Args:
             replicates: zero-based replicate indices
@@ -244,6 +246,29 @@ class _MolecularDynamicsResult(WorkflowResult):
             json=replicates,
         )
 
+    def download_trajectory_files(
+        self,
+        replicates: list[int],
+        path: Path | str | None = None,
+    ) -> list[Path]:
+        """Download each selected replicate's DCD file directly from storage.
+
+        Args:
+            replicates: zero-based replicate indices
+            path: output directory; defaults to the current directory
+
+        Returns:
+            paths to the downloaded `.dcd` files, one per replicate
+
+        Raises:
+            httpx.HTTPStatusError: API or storage request fails
+        """
+        return download_presigned_files(
+            Path(path) if path is not None else Path.cwd(),
+            f"/trajectory/{self.workflow_uuid}/trajectory_dcd_urls",
+            json=replicates,
+        )
+
     def _trajectory(self, replicate: int) -> _MolecularDynamicsTrajectory:
         """Return one trajectory with a consistent bounds error."""
         trajectories = self.trajectories

rowan/workflows/binding_affinity.py

@@ -153,9 +153,15 @@ def submit_binding_affinity_workflow(
         binding_affinity_settings=binding_affinity_settings or SinglePointEnergySettings(),
     )
 
+    workflow_data = workflow.model_dump(serialize_as_any=True, mode="json")
+    if isinstance(workflow.binding_affinity_settings, SinglePointEnergySettings):
+        workflow_data["binding_affinity_settings"]["truncation_radius"] = (
+            workflow.binding_affinity_settings.truncation_radius
+        )
+
     data = {
         "workflow_type": "binding_affinity",
-        "workflow_data": workflow.model_dump(serialize_as_any=True, mode="json"),
+        "workflow_data": workflow_data,
         "name": name,
         "folder_uuid": folder_uuid,
         "max_credits": max_credits,

rowan/workflows/conformer_search.py

@@ -31,6 +31,8 @@ from .base import (
 )
 from .constants import to_relative_kcal
 
+ConformerProperties = stjames.ConformerProperties
+
 
 @register_result("conformer_search")
 class ConformerSearchResult(WorkflowResult):
@@ -72,22 +74,59 @@ class ConformerSearchResult(WorkflowResult):
         return to_relative_kcal(energies) if relative else energies
 
     @property
-    def radii_of_gyration(self) -> list[float]:
-        """Radius of gyration for each conformer (A)."""
+    def radii_of_gyration(self) -> list[float | None]:
+        """Radius of gyration for each conformer (A), None where it failed."""
         return [p.radius_of_gyration for p in self._workflow.conformer_properties]
 
     @property
-    def sasa(self) -> list[float]:
-        """Solvent accessible surface area for each conformer (A^2)."""
+    def sasa(self) -> list[float | None]:
+        """Solvent accessible surface area for each conformer (A^2), None where it failed."""
         return [p.solvent_accessible_surface_area for p in self._workflow.conformer_properties]
 
     @property
-    def polar_sasa(self) -> list[float]:
-        """Polar solvent accessible surface area for each conformer (A^2)."""
+    def polar_sasa(self) -> list[float | None]:
+        """Polar solvent accessible surface area for each conformer (A^2), None where it failed."""
         return [
             p.polar_solvent_accessible_surface_area for p in self._workflow.conformer_properties
         ]
 
+    @property
+    def dipoles(self) -> list[float | None]:
+        """Dipole moment magnitude about the center of mass for each conformer (Debye).
+
+        None where the method reports no dipole.
+        """
+        return [p.dipole for p in self._workflow.conformer_properties]
+
+    @property
+    def ensemble_properties(self) -> ConformerProperties | None:
+        """Ensemble properties weighted by `boltzmann_weights`.
+
+        Each property is averaged over the conformers that have it: SASA, polar SASA, Rg,
+        the root-mean-square dipole (Debye), and per-atom Mulliken charges. When the run
+        computed them, it also holds stick spectra with every conformer's lines, each
+        intensity scaled by its conformer's population: `excitation_energies` (Hartree) with
+        `oscillator_strengths` (UV-vis) and `rotatory_strengths` (ECD), and
+        `vibrational_frequencies` (cm^-1) with `ir_intensities` (km/mol). Pass them to
+        `rowan.broaden_spectrum` for the spectrum. Rotatory strengths are omitted when the
+        conformers differ in stereochemistry, since mixing enantiomers cancels the ECD signal.
+        """
+        return self._workflow.ensemble_properties
+
+    @property
+    def boltzmann_weights(self) -> list[float]:
+        """Population of each conformer at the weighting temperature (sum to 1.0).
+
+        From free energies when every conformer has frequencies, else from ground-state
+        energies. Empty for workflows run before weights were reported.
+        """
+        return list(self._workflow.boltzmann_weights)
+
+    @property
+    def weighting_temperature(self) -> float:
+        """Temperature of `boltzmann_weights` and `ensemble_properties` (K)."""
+        return self._workflow.weighting_temperature
+
     def get_conformers(self, n: int | None = None) -> list[Molecule]:
         """Fetch conformer molecules.
 
@@ -156,6 +195,7 @@ def _mso_for_final_method(
     final_method: Method,
     solvent: SolventInput = None,
     transition_state: bool = False,
+    frequencies: bool = False,
 ) -> MultiStageOptSettings:
     """Build optimization settings for the final method.
 
@@ -163,9 +203,9 @@ def _mso_for_final_method(
     without a named preset.
 
     `solvent` and `transition_state` are merged into each stage's `Settings`
-    when supplied; the named preset itself is independent of them (mirrors how
-    tinbergen's form keeps the preset dropdown separate from the solvent / TS
-    toggles).
+    when supplied, and `frequencies` runs on the optimization stage; the named
+    preset itself is independent of them (mirrors how tinbergen's form keeps the
+    preset dropdown separate from the solvent / TS toggles).
     """
     opt_method, sp_method = _FINAL_METHOD_TO_MSO_RECIPE.get(final_method, (final_method, None))
 
@@ -196,6 +236,7 @@ def _mso_for_final_method(
     return MultiStageOptSettings(
         optimization_settings=[opt_stage],
         singlepoint_settings=sp_stage,
+        frequencies=frequencies,
     )
 
 
@@ -224,6 +265,8 @@ def submit_conformer_search_workflow(
     final_method: Method | str = "aimnet2_wb97md3",
     solvent: SolventInput = None,
     transition_state: bool = False,
+    frequencies: bool = False,
+    weighting_temperature: float = 298.15,
     multistage_opt_settings: MultiStageOptSettings | None = None,
     conformer_clustering_settings: ConformerClusteringSettings | None = None,
     initial_conformers: Sequence[StructureInput] | None = None,
@@ -264,6 +307,11 @@ def submit_conformer_search_workflow(
             `multistage_opt_settings` is provided
         transition_state: whether to optimize the transition state. Ignored
             if `multistage_opt_settings` is provided
+        frequencies: compute vibrational frequencies of each conformer at the optimization
+            level, for the ensemble IR spectrum. Ignored if `multistage_opt_settings` is
+            provided; set its `frequencies` instead
+        weighting_temperature: temperature for the conformer populations, in Kelvin. With
+            frequencies, populations use free energies at this temperature
         multistage_opt_settings: optimization stages and singlepoint settings
             for ranking conformers. When provided, takes precedence over
             `final_method` / `solvent` / `transition_state`. When omitted, an MSO is
@@ -359,7 +407,10 @@ def submit_conformer_search_workflow(
         if isinstance(final_method, str):
             final_method = Method(final_method)
         multistage_opt_settings = _mso_for_final_method(
-            final_method, solvent=solvent, transition_state=transition_state
+            final_method,
+            solvent=solvent,
+            transition_state=transition_state,
+            frequencies=frequencies,
         )
     multistage_opt_settings = _release_ts_optimization_constraints(multistage_opt_settings)
 
@@ -370,6 +421,7 @@ def submit_conformer_search_workflow(
         multistage_opt_settings=multistage_opt_settings,
         conf_gen_settings=conf_gen_settings,
         conformer_clustering_settings=conformer_clustering_settings,
+        weighting_temperature=weighting_temperature,
     )
 
     workflow_data = workflow.model_dump(serialize_as_any=True, mode="json")

rowan/workflows/protein_md.py

@@ -149,7 +149,7 @@ def submit_protein_md_workflow(
     binder: Binder | None = None,
     protein_restraint_cutoff: float | None = None,
     protein_restraint_constant: float = 100,
-    analysis_interval_ps: float | None = None,
+    analysis_interval_ps: float | None = 100,
     clustering: KMeansClusteringSettings | GreedyClusteringSettings | None = None,
     validate_forcefield: bool = True,
     name: str = "Protein MD Workflow",
@@ -194,8 +194,8 @@ def submit_protein_md_workflow(
         protein_restraint_cutoff: cutoff distance from the binder past which Cα atoms are
             harmonically restrained, in Å; None disables restraints
         protein_restraint_constant: force constant for Cα backbone restraints, in kcal/mol/Ų
-        analysis_interval_ps: interval at which to compute per-frame SASA and polar SASA, in ps.
-            None disables those analyses
+        analysis_interval_ps: interval at which to compute per-frame SASA and polar SASA
+            (default 100), in ps. None disables those analyses
         clustering: how to cluster trajectory frames. None disables clustering; pass a
             KMeansClusteringSettings (num_clusters) or GreedyClusteringSettings (cutoff_angstrom)
         validate_forcefield: validate the protein forcefield

rowan/workflows/relative_binding_free_energy_perturbation.py

@@ -12,7 +12,7 @@ from rowan.folder import Folder
 from rowan.molecule import Molecule
 from rowan.protein import Protein
 from rowan.types import ProteinUUID
-from rowan.utils import api_client, download_file
+from rowan.utils import api_client, download_file, download_presigned_files
 
 from .base import Workflow, WorkflowResult, molecule_to_dict, register_result
 from .rbfe_graph import RelativeBindingFreeEnergyGraphEdge, RelativeBindingFreeEnergyGraphResult
@@ -143,7 +143,9 @@ class RelativeBindingFreeEnergyPerturbationResult(WorkflowResult):
         *,
         leg: Literal["complex", "solvent"] = "complex",
     ) -> Path:
-        """Download DCD trajectory files for a specific perturbation edge.
+        """Download DCD trajectory files for a specific perturbation edge as one `.tar.gz` archive.
+
+        Long trajectories can time out; use `download_edge_trajectory_files` for those.
 
         Args:
             edge_index: index of the edge (0-based, matching `edges` order)
@@ -183,6 +185,40 @@ class RelativeBindingFreeEnergyPerturbationResult(WorkflowResult):
             json=lambda_vals,
         )
 
+    def download_edge_trajectory_files(
+        self,
+        edge_index: int,
+        lambda_vals: list[float] | None = None,
+        path: Path | str | None = None,
+        *,
+        leg: Literal["complex", "solvent"] = "complex",
+    ) -> list[Path]:
+        """Download each lambda window's DCD file for one edge directly from storage.
+
+        Args:
+            edge_index: index of the edge (0-based, matching `edges` order)
+            lambda_vals: lambda values to download. Defaults to all windows
+            path: directory to save the files to. Defaults to current directory
+            leg: thermodynamic leg whose trajectories should be downloaded
+
+        Returns:
+            paths to the downloaded `.dcd` files, one per lambda window
+
+        Raises:
+            IndexError: edge_index is out of range
+            httpx.HTTPStatusError: API or storage request fails
+        """
+        edges = self.edges
+        if edge_index < 0 or edge_index >= len(edges):
+            raise IndexError(f"Edge index {edge_index} out of range (0-{len(edges) - 1})")
+
+        return download_presigned_files(
+            Path(path) if path is not None else Path.cwd(),
+            f"/trajectory/{self.workflow_uuid}/rbfe_trajectory_dcd_urls",
+            params={"edge_index": edge_index, "leg": leg},
+            json=lambda_vals,
+        )
+
     def download_all_trajectories(
         self,
         path: Path | str | None = None,
@@ -205,6 +241,35 @@ class RelativeBindingFreeEnergyPerturbationResult(WorkflowResult):
             self.download_edge_trajectories(i, path=path, leg=leg) for i in range(len(self.edges))
         ]
 
+    def download_all_trajectory_files(
+        self,
+        path: Path | str | None = None,
+        *,
+        leg: Literal["complex", "solvent"] = "complex",
+    ) -> list[Path]:
+        """Download every edge's DCD files directly from storage.
+
+        Files for each edge go in an `edge_<index>` subdirectory of `path`.
+
+        Args:
+            path: directory to save the files to. Defaults to current directory
+            leg: thermodynamic leg whose trajectories should be downloaded
+
+        Returns:
+            paths to all downloaded `.dcd` files, in edge order
+
+        Raises:
+            httpx.HTTPStatusError: any API or storage request fails
+        """
+        directory = Path(path) if path is not None else Path.cwd()
+        return [
+            file
+            for i in range(len(self.edges))
+            for file in self.download_edge_trajectory_files(
+                i, path=directory / f"edge_{i}", leg=leg
+            )
+        ]
+
     @property
     def ligand_dg_results(self) -> dict[str, RelativeBindingFreeEnergyResult] | None:
         """Per-ligand binding free energy results, or None if not yet computed."""

{rowan_python-3.2.1.dist-info → rowan_python-3.2.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.5
 Name: rowan-python
-Version: 3.2.1
+Version: 3.2.3
 Summary: Rowan Python Library
 Project-URL: Homepage, https://github.com/rowansci/rowan-client
 Project-URL: Bug Tracker, https://github.com/rowansci/rowan-client/issues
@@ -12,7 +12,7 @@ Requires-Dist: httpx
 Requires-Dist: nest-asyncio
 Requires-Dist: rdkit
 Requires-Dist: setuptools
-Requires-Dist: stjames>=0.0.275
+Requires-Dist: stjames>=0.0.289
 Description-Content-Type: text/markdown
 
 # Rowan Python Library

{rowan_python-3.2.1.dist-info → rowan_python-3.2.3.dist-info}/RECORD

@@ -1,26 +1,27 @@
-rowan/__init__.py,sha256=inUxeeRp1yF7Jwn-_ZRjBBgmaZGlxY-9kmz_a0ejECE,2111
+rowan/__init__.py,sha256=pXnVxfLHs7s7AZZJmcaMk8CjTYz-g-IhMxjIV0B7i8Y,2224
 rowan/api_keys.py,sha256=ms8NAkLreFQcaEUPCf8lSf7T0x2OmrWCAT44N1cgocE,6955
 rowan/calculation.py,sha256=E1dAox4ScE3NHJXiKFZ3Lpmw2h7DKCVzcwRvGZ2H2cQ,3802
 rowan/config.py,sha256=ay6UxKLE9bbTxHFXxDfUXGRVJWJ1T-eS2_L4HF8mb_E,20927
 rowan/constants.py,sha256=emCH4m9OL2Hm5E-6mJGM_FgzrK_JrZT-FiKJ6pMNQ4Y,84
 rowan/folder.py,sha256=8Gb60oes78G7QTGm0EX-ofkmQEggv5cwseA7dRBIhok,13883
 rowan/molecule.py,sha256=1ueC_M2Yaubkt5eQlONCagShWPnoJfNNXgd0s7sslPM,15826
-rowan/project.py,sha256=N8DdNt7xCYHunIAqTvPmBGkH8thzwmguATo4Thfa13o,5479
+rowan/project.py,sha256=R5uwpzg7GOaGQ_ebCvTysYvkUHE_hHxWAhujp7CBqU0,5770
 rowan/protein.py,sha256=AGHarvWfciSwWXfHnVjdUW47FML7vI4Kk4dYTeFfKAQ,19322
 rowan/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+rowan/spectra.py,sha256=dNmE5th6wpTVAfvZZNDkg2jlBRAmPtFT6yDtxex8APw,2307
 rowan/types.py,sha256=D0d6MW6MvLMs6UtNgGKlapFGOnPlzqFTk9EAgB9dwks,636
 rowan/user.py,sha256=wGwyp0uzifzJeC88yCC7pgxMhxuZ-LLWHfCcuPXcx44,6199
-rowan/utils.py,sha256=d5s7UzjgJNtZ_NDLjvXSeSVnEBrKOqQiyVy-7C13rm0,8020
+rowan/utils.py,sha256=n82jhtMrHwIPXjEX4i2JFuQTLsanwixk5HfOfzPpKko,9656
 rowan/workflows/__init__.py,sha256=1JBuqs3KlizJEFhidvGCG7KtTdLOQKI2kspn0SK_KHw,5161
-rowan/workflows/_molecular_dynamics.py,sha256=ZkddNEo-aRzwFRopDh6LcPPRDL6Rrx4lwW-_roMW3Kc,9323
+rowan/workflows/_molecular_dynamics.py,sha256=fDQVAZam6Ic_Au1yz41bVWg2_NBfkPX2ZFJNDIbw55U,10204
 rowan/workflows/admet.py,sha256=92jEbZr9OuN0FzJtWmE4EdARzs6r0uiRwq3-kSPV7nE,2531
 rowan/workflows/analogue_docking.py,sha256=NZa-2rNzuyfVE_5WZm4pQspdtWFy8tEhxbqwOrXWGnU,10881
 rowan/workflows/base.py,sha256=DR2k5VFMmv2sa9v4hYU_Fvzf8k2unp5nhPqvOd7NX64,39520
 rowan/workflows/basic_calculation.py,sha256=7o9jxEpDMiong2k3Ke1WCLwaSWZvKY5tqhbixnbQWbo,12945
 rowan/workflows/batch_docking.py,sha256=AtCp2-D_LOVxzSPmQ144zYIz5kft9NeRljELdRDU6vY,4941
 rowan/workflows/bde.py,sha256=4Yw27U2RPsuUDf9tw-7wPbw_4gaUcvDlQP8je7VJktU,7099
-rowan/workflows/binding_affinity.py,sha256=O0fL_oQtMpYlmdaCvHTR7WSNfOENMPAN0OUBxO4640k,6553
-rowan/workflows/conformer_search.py,sha256=enx9HbFZVVYgHJP-HUFam63E0qCDq6jjGeXNuBXXlys,16311
+rowan/workflows/binding_affinity.py,sha256=CP5rsWNMhsdLGzLdPBIs89wZDW6sa8YwGEkqep0ls1Q,6821
+rowan/workflows/conformer_search.py,sha256=PSufIThwh2LoqgP6x5gjFxuNW5DP5ClJ5caBCuJoDKU,18910
 rowan/workflows/constants.py,sha256=iM9-6DyIe0N2ydyq-lRYP7VDufFmMrN5tDjnpUSlxF8,588
 rowan/workflows/covalent_inhibitor_scan.py,sha256=VcZhwW3HRuFuGZNedMPZBWbzuhGvz_B5pNJSIPyi9Kw,6721
 rowan/workflows/descriptors.py,sha256=mLo_VnqVrnEbiMHT1-SrbKhy_Dk5jpIM4btNGrlmNio,3041
@@ -43,18 +44,18 @@ rowan/workflows/pocket_detection.py,sha256=dFeCLB59v5-Izf54blrRzJhwUtBQX2SJyVKT_
 rowan/workflows/pose_analysis_md.py,sha256=ndfGS8oP8oQ3EXe83Wgbx1H8ZFURlm8IJZN34aDpqeE,10599
 rowan/workflows/protein_binder_design.py,sha256=ee3qYGzTzQsca9otnRz9Qjwj1X66LTCKqtTQ5EOzBHk,9181
 rowan/workflows/protein_cofolding.py,sha256=s4bLxlix4kAO_NzS75tc2wiEYqUxpFofxziUnxI2iUY,17422
-rowan/workflows/protein_md.py,sha256=M4lbfEEvyxmhBWfz_Hf3whmNIYuHe3Q3HZzDcq-jYII,12264
+rowan/workflows/protein_md.py,sha256=iDeKW0ONY3ZPKPsjmZVQhbFQbo0vmOp5M8ko2Hfc16U,12277
 rowan/workflows/protein_preparation.py,sha256=8C04rhR4AiGiYJWZ8V06O4kpV443tGihcTk8rnNgscE,5325
 rowan/workflows/rbfe_graph.py,sha256=DBNxtmgw7BNSwOXFLpHDxJtfpheTpDPs8U3SYB2MXdE,8070
 rowan/workflows/redox_potential.py,sha256=0AV2nKTngEFw7VDDpDQvdeBo5u03Aklc2ak4bDWxUz0,5366
-rowan/workflows/relative_binding_free_energy_perturbation.py,sha256=dBQkM3wtbNVVkiXbB3RBBDUcme2OnWZkWFZdTUR8vlM,15872
+rowan/workflows/relative_binding_free_energy_perturbation.py,sha256=KDVSKBqHQcL4zsjad8abmcT1Y-viXdmHedhVW2Xg55M,18319
 rowan/workflows/scan.py,sha256=MqV0NcjfInwxGlegOyWvAt8pwN0BdpsKvgCgTRkzdzU,7440
 rowan/workflows/solubility.py,sha256=Lws0lk4ifmtj84W_jyeOuaTpt8KM8-fceOggHrEYnmc,11211
 rowan/workflows/solvent_dependent_conformers.py,sha256=1dYf8UBVMRFELdwHzv3MI0pe12qoCAfohkaMiW3ak50,7460
 rowan/workflows/spin_states.py,sha256=C3hJEY4hccE1fCW4UjI4nIqGZLw9l2C2VKWCLiopSlY,9044
 rowan/workflows/strain.py,sha256=hq0dOlxnN72V4Sohr3qNFbOrjGmV01tc5zcF0HrL2-Y,6372
 rowan/workflows/tautomer_search.py,sha256=o5iFXHfcb6DvvWH7itA05hvWIrzxrc3a_9IDikYftiM,6351
-rowan_python-3.2.1.dist-info/METADATA,sha256=g0STX2NMLAvsvtvLRpJ8gaR6qBV-9SckC_SQzSaJzE0,3934
-rowan_python-3.2.1.dist-info/WHEEL,sha256=W3fkpkm7-wf9vBI5Z-7s0eWkeM-spu78I8Neb98DeEg,87
-rowan_python-3.2.1.dist-info/licenses/LICENSE,sha256=i05z7xEhyrg6f8j0lR3XYjShnF-MJGFQ-DnpsZ8yiVI,1084
-rowan_python-3.2.1.dist-info/RECORD,,
+rowan_python-3.2.3.dist-info/METADATA,sha256=5zYnTTmtF9VpgtLBN_VPLY2RCU_v42bNunjmoWaqr_U,3934
+rowan_python-3.2.3.dist-info/WHEEL,sha256=W3fkpkm7-wf9vBI5Z-7s0eWkeM-spu78I8Neb98DeEg,87
+rowan_python-3.2.3.dist-info/licenses/LICENSE,sha256=i05z7xEhyrg6f8j0lR3XYjShnF-MJGFQ-DnpsZ8yiVI,1084
+rowan_python-3.2.3.dist-info/RECORD,,